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Magnetically tunable Feshbach resonances were employed to associate cold diatomic molecules in a series of 
experiments involving both atomic Bose as well as two spin component Fermi gases. This review illustrates 
theoretical concepts of both the particular nature of the highly excited Feshbach molecules produced and the 
techniques for their association from unbound atom pairs. Coupled channels theory provides the rigorous for- 
mulation of the microscopic physics of Feshbach resonances in cold gases. Concepts of dressed versus bare 
energy states, universal properties of Feshbach molecules, as well as the classification in terms of entrance- and 
closed-channel dominated resonances are introduced on the basis of practical two-channel approaches. Their sig- 
nificance is illustrated for several experimental observations, such as binding energies and lifetimes with respect 
to collisional relaxation. Molecular association and dissociation are discussed in the context of techniques involv- 
ing linear magnetic field sweeps in cold Bose and Fermi gases as well as pulse sequences leading to Ramsey-type 
interference fringes. Their descriptions in terms of Landau-Zener, two-level mean field as well as beyond mean 
field approaches are reviewed in detail, including the associated ranges of validity. 
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IS I. INTRODUCTION 



Resonances in general refer to the energy dependent en- 
hancement of inter-particle collision cross sections due to the 
existence of a meta-stable state. Since the beginning of quan- 
tum mechanics such phen omena have been the subject of nu- 
merous studies in nuclear (Blatt and Weisskopl 19521) as well 
as atomic and m olecular physics dBransden an d Joachain, 
12003b IChildL Il974h . The meta-stable state may be described 
in terms of tunnelling across a barrier of the potential energy 
or coupling of a bound level of a subsystem to its envir onment 
dFanoll 19351 fT%ll iFeshbachl 1 1 95 8ll 1 9621 iRicel 1 1 93 3l) . These 
scenarios are respectively referred to as shape and Feshbach 
resonances. 
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In the context of cold atomic gases, collision phenom- 
ena associated with Feshbach resonances were predicted 
first for systems of spin polarised hydrog en and deuterium 
(Stwalley, 1976; Ua ng and Stwallevill980l) as well as lithium 
(lUang et all 1 1 98 lb in the presence of magnetic fields. The 
associated resonance energies depend on the field strength 
via the Zeeman effect in the hyperfine levels. This re- 
search has gained substantially new experimental perspec- 
tives since the achievement of Bo se-Einstein condensation 
dBosel Il924t lEinsteinl 1 1924 1 19251) of dilute vapours of al- 
kali atoms J Anderson et adll995HBradlev et al.lll997l.ll995t 



16000 



et at 



1995). Contrary to conventional gases, such 



iDavis 

atomic clouds with densities five orders of magnitude less 
than air and sub yuK temperatures give rise to binary colli- 
sion energies close to the threshold between scattering and 
molecular binding. In this extraordinary regime, magneti- 
cally tunable Feshbach resonances can be employed to ma- 
nipulat e the inter-atomic for ces determined by the scattering 
length dTiesinga et all Il993[). as well as for the produc tion of 
diatomic molecules at rest (iTimmermans et all [l999b). 



— van Abeelen & Verhaar 
• Inouyeefa/. 




895 900 905 910 

magnetic field strength B [G] 



915 



FIG. 1 (Colour in online edition) Order of magnitude variation of the 
scattering length, a, as a function of the magnetic field strength, B, 
in the vicinity of the 907 G (90.7 mT) zero energy resonance of 23 Na. 
The circl es refer to measurem ents using sodium Bose-Einstein con- 
densates dlnouve ega/.LIl998l ) , while the solid curve indicates asso- 
ciated theoretical predictions ( Ivan Abeel en and Verhaarl [19981) . The 
scattering length is normalised to its asymptotic value a\, s away from 
the singularity at Bo = 907 G. Adapte d by permission fro m Macmil- 
lan Publishers Ltd: Nature (London), dlnouve et a/.LIl998h . copyright 
(1998). 

Both techniques were demonstrated in several series of 
experiments with widespread applications throughout the 
physics of cold gases. Figure [TJ illustrates the manipulation 
of the scattering length, a, in a 23 Na Bose-Einstein conden- 
sate exposed to a s patially homogeneo us magnetic field of 
variable strength, B dlnouve et all [1998!). The pole in Fig. [TJ 
at about 907 G 1 is due to the near degeneracy of the en- 




FIG. 2 (Colour in online edition) Ramsey fringes between atomic 
and molecular components produced from a Bo se-Einstein conden- 
sate of 17100 85 Rb atoms dPonlev et all 12002) via a sequence of 
magnetic field pulses in the vicinity of the 155 G zero energy reso- 
nance. Circles and squares refer to the number of particles in the rem- 
nant condensate and in an atomic burst consisting of correlated pairs 
with a comparatively high relative velocity, respectively. Diamonds 
indicate the total amount of particles in both of these components. Its 
difference with respect to the total number of atoms (dotted line) indi- 
cates the production of undetected Feshbach molecules. The fringes 
are shown as a function of the delay time of the interferometer, fevoive, 
in which the atomic and molecular states acquire their phase dif- 
ference. The frin ge frequency provid es an accurate measure of the 
binding energy (Claussen et ah, 2003). Adapte d by permission from 
Macmillan Publishers Ltd: Nature (London), dPonlev et all 1 20021) . 
copyright (2002). 



ergy associated with a Feshbach resonance and the thresh- 
old. Such a singularity, u sually referred to as a zero energy 
resonance dTavlorl Il972l) 2 . allows the scattering length, in 
principle, to assume all values between -cm and oo. Cold 
gases with such widely tunable interactions were subse- 
quentl y realised for several s pecies of alkali atom s, such as 

~ 3 Cs 



quentiy reatise d tor several s pecies oi alkali atom s, si 
85 Rb dCourteille et all [19981 (Roberts et all 1 1998b and 



dChin et all 120031: 1 Vuletic et aZ.Ul998l) . Their applications in- 
volve the Bose-Ei nstein condensation of 85 Rb dCornish et all 
l2000t) and 133 Cs dWeber et all 120031) as well as studies of 
the collapse of condensates with negative scattering lengths 
dPonlev et aZ.ll200UlRoberts et a/.l,l2000al) . 

Quite generally, the zero energy resonance position exactly 
coincides with the field strength at which the energy of a di- 
atomic vibrational bound state becomes degenerate with the 
threshold for dissociation into an atom pair at rest. On the 
side of positive scattering lengths, this Feshbach molecular 
state describes a stable molecule in the absence of background 
gas collisions which ceases to exist at the position of the sin- 
gularity. In the context of two-body systems involving the 



1 The Standard International unit for magnetic field is Tesla, whereas most 



of the papers quoted in this review use Gauss as the unit. Consequently, we 
use G here, where 1 G=10 T. 

In the context of cold gases, a singularity of the scattering length is also 
often referred to simply as a Feshbach resonance. 
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same species of atoms such a bound state is usually referred 
to as a dimer. The relation between collision resonances above 
threshold and bound states below it links the manipulation 
of interactions to the molecular conversion of separated atom 
pairs. 

Production of cold dimers w as de monstrated first 

dFioretti et all T l998; Ni kolov et all 1 19991: iTakekoshi etail 
[19 



1999; Wvnar et al. 



(Weineref al., 199 



2000) via photo-association of atoms 
. This achievement was followed by 
studies of condensed gases of 85 Rb exposed to time dependent 
magnetic field variations consisting of pairs o f pulses in the 
vicini ty of the 155 G zero energy resonance (iDonlev et all 
2002). These experiments temporarily probed the regime of 
strong interactions where the size of the scattering length 
was comparable to the average inter-atomic distance. Such 
perturbations led to the three distinct components of the 
gas illustrated in Fig. [2] The oscillatory behaviour of their 
occupations as a function of the time delay between the pulses 
implied an interpretation in terms of Ramsey interference 
fringes due to a superposition state of separated atoms and 
Feshbach molecules. 




FIG. 3 (Colour in online edition) A molecular component of about 
3000 dimers is out-coupled from a dilute cloud of 25000 ground state 
cesium atoms using Stern-Gerlach separation in an inhomogeneous 
magnetic field (Herbig et al., 2003). The left and middle images re- 
fer to situations in which the magnetic field is calibrated in such 
a way that it exactly compensates for the gravitational force acting 
on the atoms in gases without and with a molecular component, re- 
spectively. Due to their magnetic moment difference of -0.57 //Boh, 
(//Bohr = 9.27400949X 10~ 24 J/T is the Bohr magneton) with respect to 
separated atoms, the Feshbach molecules in the middle image drop 
down from the atomic cloud which is levitated and centred at the 
same position as in the left reference image. Conversely, the right 
image shows levitation of molecules and upward acceleration of the 
separated atoms using a suitably adjusted inhomogeneous magnetic 
field. 

According to Fig. |2j the pulse sequence allowed a con- 
version of up to about 16% of the 85 Rb atoms into Fesh- 
bach molecules. Subsequent experiments improved the pro- 
duction efficiency by using magnetic field sweeps from neg- 
ative to positive scattering lengths across a zero energy res- 



onance. This technique was applied to cold gases consist- 
ing of incoher ent two spin component mixtures of ei t her 4 K 
or 



6 Li atoms dCubizolles et al. ,_ 2003 ;_ Jochim et all l2003at 
"ID 12003a 



Regal et all 2003a; IStrecker et all 120031) as we ll as Bose- 
Einste i n condensate s of 133 Cs, 87 Rb and 2 3 Na dDiirr et all 
l2004at iHerbig et all 12003b IXu et all l2003h . In the context 
of these experiments, new schemes for the detection of Fesh- 
bach molecules were developed. These techniques involve ra- 
dio frequency (rf) photo-dissociation, atom loss and recovery, 
as well as the spatial separation of molecules from the rem- 
nant atomic cloud followed by their dissociation using mag- 
netic field sweeps. Separation of Feshbach molecules from 
an atomic gas, for instance, may be achi e ved v i a the Stern- 
Gerlach approach of Fig. [3] dChin et all 12005b iDurr et all 
l2004iJ: IHerbig et all 120031) . probing the magnetic moments 
of dimers at magnetic fields away from the zero energy reso- 
nance. 




FIG. 4 (Colour in online edition) Density distributions of dilute ther- 
mal (left image) and partially Bose-Ei nstein condensed (rig ht im- 
age) gases of 40 K2 Feshbach molecules iGreiner et al ., 2003). Both 
molecular vapours were produced using the technique of linear mag- 
netic field sweeps across the 202 G zero energy resonance. The initial 
two spin component clouds of 470000 and 250000 Fermi atoms were 
prepared above (250 nK) and below (79 nK) the associated critical 
temperatures for condensation, respectively. The condensate fraction 
in the right image is 12%. Adapte d by permission from Macmillan 
Publishers Ltd: Nature (London), jGreiner et all 120031) . copyright 
(2003). 

Near reso nance, their l a rge spa tial extent associated with a 
loose bond dKohler et all 12003 a b) can lead to a remarkable 
stability of Feshbach molecules with respect to inelastic colli- 
sions in the environment of a two spin component Fermi gas 
dPetrov et all 12004 l2005bl) . The lifetimes of s uch dimers , 
ranging from about 100 ms in the case of 4 K. dReeal et all 
2004al) to several seconds in 6 Li2 gases dCubizolles et al. 



20031; iJochim et all l2003ah . were sufficient for the observa- 
tion of their Bose-Einstein condensation. This achievement 
was based on eit her the magne t ic fiel d sweep technique illus- 
trated in Fig. [4] dGreiner et all 20031) or evaporative cooling 
of a m olecular cloud (Jochim et al, 2003b; IZwierlein et all 
120031) . Such pionee ring experiments gave rise to an ongo 



2004 1 [Chin et al.l 
2002[|O , Haral7al. 
2004bP 



2004a 



ing series of stu dies (Bartenstein et al., 2004 



Bourdel et al. 
~r — 



Kinastefa/., 2004; Loftu s et al. 
20021: iPartridge et al.1 12005: Regal et al. 



Zwierlein et all 120051 [2004) probing the crossover 



from Bardeen-Cooper-Schrieffer (BCS) pairing at negative 
scattering leng ths to molecular Bose-Einstein condensation 
(Eagles], 1 



Il969t iLeggettl Il980t iNozieres and Schmitt-Rinkl 
1 1 985I: iRanderial 1 1 995l) . 
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This article reviews the theoretical background of the ex- 
otic species of highly excited Feshbach molecules produced 
as well as their formation in the environment of a cold gas. 

Section HT1 briefly introduces the concept of universality of 
weakly bound dimers, their relevant length and energy scales 
as well as the general form of their wave functions. This dis- 
cussion integrates near resonant Feshbach molecules into the 
general class of quantum halo systems whose classic exam- 
ples are the deuteron of nuclear physics and the helium dimer 
4 He2 molecule. More details about such exotic two-particle 
states as well as their extensions to few-b ody systems are 
given in a recent review dJensen et al 

Section |nl] discusses those concepts of diatomic scattering 
and molecular physics that are particular to the applications 
of Feshbach resonances in cold gases. Subsection [OLA] in- 
troduces the microscopic origin of the resonance enhance- 
ment of interactions illustrated in Fig. Q] scattering chan- 
nels and rigorous coupled channels theory (IMies et all \ 199 A 
Stoof et al. ., 1988), as well as the relation between zero en- 
ergy resonances and molecular energy spectra. Several as- 
pects of the rigorous method are well recovered in terms of 
two-channel approaches whose general concepts, such as the 
two-channel Hamiltonian and the meta-stable Feshbach reso- 
nance state, are discussed in Subsection IIII.BI On this basis, 
Subsection IIII.CI introduces the bare and dressed bound and 
continuum energy levels. The universal properties of Fesh- 
bach molecules near resonance are strictly derived in Subsec- 
tion llll.Dl Their physical significance is illustrated in Subsec- 
tion [TILE] for several experimentally relevant examples, such 
as molecular binding energies and the lifetimes of dimers in 
cold Fermi and Bose gases. The size of the magnetic field 
range associated with universality implies a distinction be- 
tween entrance- and closed-channel dominated zero energy 
resonances, whose physical origin is discussed in Subsec- 
tion IIII.FI Implementations of two-channel approaches are 
given in detail in Subsection IIII.GI describing properties of 
Feshbach molecules close to as well as away from zero en- 
ergy resonances. Their applications involve the Stern-Gerlach 
separation of dimers shown in Fig. [3] Characteristic parame- 
ters relevant to two-channel approaches, such as, for instance, 
the magnetic moments associated with Feshbach resonances 
of several atomic species, are summarised in Tables Hyl andlVl 

Section [TV] reviews dynamical approaches describing the 
production and dissociation of cold Feshbach molecules via 
linear magnetic field sweeps across zero energy resonances. 
Its introductory paragraphs outline, on the basis of the dis- 
crete energy spectrum of a trapped atom pair, the adiabatic 
transfer from quasi-continuum to dimer states. Molecular as- 
sociation of an atom pair via linear magnetic field sweeps 
falls into the category of rare dynamical two-body prob- 
lems whose solutions can be r epresented in an analytic form 
(iDemkov and Osherovl Il968t iMacek and Cavagnerol Il998l) . 
A detailed derivation of the associated transition amplitudes 
in Subsection llV.AI provides the foundation fo r all subsequent 
applications of the Landau-Zener approach dLandaul Il932t 
IZenerl Il932l) . These involve predictions of final populations 
of the quasi-continuum to dimer state transfer in tight atom 
traps, as well as fast sweep limits of molecule production in 



cold Bose and two spin component Fermi gases. Such con- 
version efficiencies can be sensitive to the quantum statis- 
tics associated with identical atoms whose effects are dis- 
cussed in Subsection IIV.BI The opposite saturation limits of 
molecule production via magnetic field sweeps in cold gases 
are a subject of ongoing research requiring dynamical descrip- 
tions of many-particle systems. Among such methods, the 
two-level mean field appro ach to Bose-Einstein condensates 
dTim mermans et al ., 1999b) is outlined in detail, including its 
relati o n to the associated Landau-Zener theory dGoral et all 
120041: IMies et all [2000). Intuitive as well as quantitative 
methods in the context of molecular Bose-Einstein conden- 
sation illustrated in Fig. |4] are reviewed in Subsection IIV.CI 
outlining the concept of an adiabatic production of dimers via 
magnetic field sweeps in cold gases with a significant momen- 
tum spread. Subsection llV.DI discusses the theory of Feshbach 
molecular dissociation, demonstrating the accuracy of ana- 
lytic treatments of linear magnetic field sweeps across zero 
energy resonances. 

Section [V] addresses the production of dimers via non- 
linear magnetic field variations, illustrated, for instance, in 
Fig. [2] in the context of Ramsey interferometry with atoms 
and Feshbach molecules. The description of these experi- 
ments requires general theoretical concepts, such as a precise 
treatment of molecular populations in gases, which are ob- 
scured in many applications involving linear magnetic field 
sweeps. Subsection lV.AI gives an intuitive explanation for the 
observations of the Ramsey fringes of Fig. [2] using a simpli- 
fying two-body approach. The general observable describ- 
ing the number of m olecules in the environment of a gas 
(Kohleref al., 2003a) is introduced in Subsection IV.BI This 
introduction includes simple applications to the fast sweep 
limits of dimer production in Bose and two spin component 
Fermi gases. Subsection [VC] addresses the complete descrip- 
tion of the Ramsey interferometry experi ments in terms of dy- 
namical beyond mean fiel d approaches dHolland et all l200lt 
iKohler and Burnetii 120021) . 

Section [VI] concludes this review and provides an outlook 
on related research, such as p-wave and optical Feshbach res- 
onances, cold dipolar and ground state molecules as well as 
extensions to few-body physics and studies of Efimov's effect. 



II. WEAKLY BOUND DIATOMIC MOLECULES 

Weakly bound diato mic molecule s are s pecial cases of low 
energy halo systems dJensen et all |2004|) . These remark- 
able quantum states are characterised by a large mean sepa- 
ration between the constituent particles which very much ex- 
ceeds the outer classical turning point, r c i ass ; ca i, determined 
by their binding energy and their attractive interaction. Ha- 
los are therefore quite distinct from spatially extended quasi 
classical states, such as the Rydberg levels associated with 
the long range Coulomb potential. Due to their non-classical 
but comparatively simple nature, two-body halos have played 
a significant role in the present understanding of compos- 
ite few-particle syst ems. Classic exampl e s are the deuteron 
in nuclear physics (Blat t and Weisskopn Il952l) as well as 
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the weakly bound helium dimer 4 He? molecule dLuo et a El 
1 19931 ISchollkopf and Toenniesl 1 19941) . whose typical halo 
wave function is illustrated in Fig. [5] We note that the modulus 
of the bound state energy of this halo molecule with respect to 
the dissociation threshold, i.e. the zero of energy in Fig. [5] is 
negligibly small compared to the well depth of the pair inter- 
action. The 4 He atoms are separated by distances larger than 
the classical turning point, r i ass i ca i, with a probability of about 
80 %. Such distances very much exceed the intuitive force 
range associated with the potential well. 
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FIG. 5 (Colour in o nline edition) Hel ium dimer interaction poten- 
tial (dashed curve) dTang et oil 1 19951) and radial probability den- 
sity of the 4 He2 molecule (solid curve). The dotted and dot-dashed 
lines indicate the bound state energy E b = -4.2 x 10~ 9 a.u. (1 a.u.= 
4.35974417 x 1(T 18 J is the atomic unit of energy) and the classical 
turning point r c i assical = 27 a Bobr (a Bobl = 0.052917 nm is the Bohr 
radius), respectively. The inter-atomic distance r is given on a loga- 
rithmic scale. 

The characteristic long range of isotropic diatomic halo 
molecules, such as the helium dimer, is determined mainly by 
a single parameter of the interaction potential, the s- wave scat- 
tering length a. Whenever the inter-atomic interaction sup- 
ports a weakly bound halo state the scattering l ength is posi- 
tive. In accordance with effective range theory dBethel 1 19491 
ISchwingeiifT947albl) . the length scale a may be interpreted in 
terms of the radius of a hypothetical hard sphere which mim- 
ics the scattering properties of the real microscopic potential 
in the limit of zero collision energy. This radius is closely re- 
lated to the bound state energy of a weakly bound diatomic 
halo molecule which is well approximated by the formula: 



E b 



-H 2 /(ma 2 ). 



(1) 



Here m is twice the reduced mass of the atom pair which in 
the case of identical particles coincides with the mass of a sin- 
gle atom. At inter-atomic distances, r, large compared to the 
classical turning point, ^classical, the associated isotropic bound 
state wave function assumes the following general form: 



<Ph(r) 



1 



,-r/a 



(2) 



J2na r 

This wave function gives the mean distance between the 



atoms, i.e. the bond length of the molecule, to be: 



<r> = 4n 



o 



r 2 drr\cf> h (r)\ 2 



a/2. 



(3) 



For the typical example of the helium dimer in Fig. [5]the so- 
lution of the stationary two -body Schr odinger equation with 
a realistic pair interaction dTang et al.l 1 1995b predicts (r) = 
5.1 nm, while the estimate of Eq. (f3]l yields (r) = 5.2 nm. For 
comparison, measurements based on the diffraction of a he- 
lium molecular b eam from a micro-fab ricated material trans- 
mission grating dGrisenti et al.l l2000h determined the bond 
length of 4 He 2 to be (r) = 5.2(4) nm. 

A dependence of physical quantities, such as those of 
Eqs. ([1) and (|2j, only on the scattering length rather than 
the details of the microscopic forces is usually referred to 
as universality. The relation between Eqs. ([]]) and fol- 
lows immediately from the stationary Schrodinger equation, 
Ho^tphif) - Eb<p\,(r), using the general two-body Hamiltonian 
associated with the relative motion of an atom pair, 



H 



2B 



n 2 , 

— V 2 + V(r). 
m 



(4) 



A typical realistic molecular potential V(r) dTang et q/.lll995l) 
is depicted in Fig. [5] for the example of the helium dimer. Its 
behaviour in the limit of large inter-atomic distances is domi- 
nated by the van der Waals interaction, 



V(r) 



-C 6 /r b 



(5) 



The constant C(, is known as the van der Waals dispersion co- 
efficient. In accordance with Eq. ©, the outer classical turn- 
ing point of the relative motion of an atom pair with the energy 
E\, of Eq. (Q} is given by the following formula: 



'"classical - [«(2/vdw) 2 ] 



(6) 



Here / v dw is a characteristic range associated with the van der 
Waals interaction between atoms, usually referred to as the 
van der Waals length, 



1 



^vdw = -(mC 6 /h ) 



2x1/4 



(7) 



The characteristic property of diatomic halo molecules that 
their spatial extent, determined by the scattering length a, 
very much exceeds r c i ass i ca i therefore implies the condition 
a » / v dw for their existence. In the range of such typ- 
ically large inter-atomic distances, r » ^ v dw, the station- 
ary Schrodinger equation for a halo molecule reduces to its 
interaction-free counterpart, i.e. 



E b (f> b (r) ~ -(ft 2 /m)V 2 4> b (r), 



(8) 



whose unit normalised solution associated with the bound 
state energy of Eq. (Q]i is given by Eq. (f2]). 

The universal properties of Eqs. (Q} and and the length 
scales, a and Z v dw, associated with diatomic halo systems char- 
acterise, in the same manner, the helium dimer as well as 
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highly excited long range Feshbach molecules. Any other 
details of their binary interactions are obscured by the large 
spatial extent of the bound states. The size of the scattering 
length, however, which determines the long range of these 
states depends sensitively on the microscopic potential whose 
detailed structure varies considerably among different species. 



III. FESHBACH RESONANCES 

Inter-atomic collisions in cold gases are characterised by 
de Broglie wavelengths much larger than the van der Waals 
length of the microscopic potential. Similarly to the universal 
properties of diatomic halo molecules, such a separation of 
length scales implies that the associated low energy interac- 
tions are determined mainly by the s-wave scattering length, 
a. Feshbach resonances provide an opportunity of manipulat- 
ing these inter-atomic forces by exposing a cold gas of alkali 
atoms to a spatially homogeneous magnetic field of strength 
B. Based on an introduction to the microscopic physics, this 
section describes the concept of magnetic tuning of the scat- 
tering length and the associated cross section (given, for in- 
stance, by 87Tfl 2 in the case of identical Bose atoms) as well as 
its relation to the properties of molecular states. 



A. Molecular physics of resonances 

Resonance enhancement of collision cross sections relies 
upon the existence of meta-stable states. The two-body 
molecular physics of such neutral atom Feshbach resonance 
states is the subject of this subsection. This involves the mi- 
croscopic origin of the inter-atomic interactions as well as the 
atomic and molecular symmetries that permit the classifica- 
tion of scattering channels and molecular bound and meta- 
stable energy levels. Such a basic understanding of low energy 
neutral atom scattering and bound states provides the ground- 
ing for practical two-channel approaches of the subsequent 
subsections. 



1. Inter-atomic interactions 

The helium atom of Section Ullhas a 'So configuration with 
no unpaired electron, and the interaction of two such atoms 
is represented by the single molecular Born-Oppenheimer po- 
tential illustrated in Fig. [5] Such a system is too simple to have 
a magnetically tunable resonance state. It is therefore nec- 
essary to consider here the experimentally relevant case of a 
pair of 2 Si/2 atoms of the same species. The unpaired electron 
spins Si and S2 from each atom can be coupled to a total spin 
S = Si + S2 with the associated quantum numbers S = or 1. 
States with S = or S = 1 are called singlet or triplet states, 
respectively. The electronic part of the inter-atomic i nterac- 
tion i s represented, as for the simplest molecule H2 JPaulingl 
Il939l) . by singlet and triplet molecular Born-Oppenheimer po- 
tentials of l Ti and symmetry, corresponding respectively 
to 25 + 1 =1 and 3. Here the notation "X" refers to zero 
projection of the electronic orbital angular momentum on the 



inter-atomic axis. The label "+" indicates that the electronic 
wave function is left unchanged upon reflection in a plane con- 
taining the nuclei. Finally, "g" and "u" are associated with the 
gerade (even) and ungerade (odd) symmetry upon inversion 
through the geometric centre of the molecule, respectively. 
Consequently, the latter symmetry is absent when the atoms 
are of different species. 
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FIG. 6 (Colour in online edition) The and 3 E+ Born- 

Oppenheimer potential energy curves for the Rbo dimer molecule 
correlating with two separated 2 Si/2 Rb atoms. An energy E asso- 
ciated with the wave number unit E/hc = 1 cm 1 corresponds to 
E/h = 29.979 GHz . The i nset s hows the long range adiabatic po- 
tentials dMies et al .1 H996L l200Cf) on a much smaller energy scale 
with the zero field (B = 0) hyperfine structure of the 85 Rb iso- 
tope with Ehf/h = 3.035 GHz. The arrows show the van der Waals 
length Z vc iw of Eq. (0 and the uncoupling distance, r u , where the 
hyperfine energy Em equals the difference between the 3 £+ and 
Born-Oppenheimer potential curves. The 87 Rb isotope has the 
same Born-Oppenheimer potential energy curves, but the long range 
curves would be different with atomic levels / = 1 and 2 and 
E M /h = 6.835 GHz. 

Figure [6] shows potential energy curves Vs (r) for 5=0 
and 1 for two Rb atoms as an illustrative case. Similar poten- 
tials exist for any pair of like alkali metal atoms, and if the 
gerade and ungerade symmetry labels g and u are dropped, 
for pairs of unlike alkali metal atoms. For the case of two 
S-state atoms, the long range form of the potential is the van 
der Waals dispersion with the lead term -Ce/r 6 of Eq. (0 that 
is identical for the and 3 £+ potentials. The splitting be- 
tween these two potentials comes from the difference in chem- 
ical bonding interactions when the charge clouds of the two 
atoms overlap at small inter-atomic distances, r < 1 nm. The 
long-range form of this splitting is associated with the elec- 
tron exchange interaction, which decreases exponentially as r 
increases. This interaction is attractive for the 'X* state and 
repulsive for the 3 X+ state. 

More diatomic states need to be accounted for when the 
atomic nuclear spin quantum number i is nonzero. At zero 
magnetic field, i.e. B — 0, the unpaired electrons in each atom 
with total electronic angular momentum j can interact with the 
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FIG. 7 (Colour in online edition) Zeeman levels of the 85 Rb 2 Si/2 
atom versus the magnetic field strength, B. Solid and dashed curves 
refer to the f = 2 and / = 3 multiplets, respectively. The projection 
quantum numbers m/ and the alphabetic labels a,b,c, ... are shown. 



nuclear spin i to produce the hyperfine levels of total angular 
momentum f = j + i. For alkali metal atoms in their electronic 
ground states the identity j = s gives the total angular momen- 
tum quantum number to be either / = i — 1/2 or i + 1/2. The 
difference in energy between these two hyperfine levels is the 
hyperfine splitting £hf- If / is an integer, the atom will be a 
composite boson, whereas if / is half an odd integer, the atom 
will be a composite fermion. The interaction with a magnetic 
field will further split the energy of each / state into its Zee- 
man components. FigureQillustrates the B-dependent energy 
levels of 85 Rb, a bosonic isotope. In the limit of high magnetic 
field strengths the Zeeman energy of the unpaired electron be- 
comes large compared to . Consequently, / will no longer 
be a good quantum number. Rotation about the magnetic field 
axis, however, is a symmetry transformation. This implies 
that nij, the quantum number associated with the projection 
of the total angular momentum onto this axis, is conserved. 
Figure [7] also shows that the Zeeman levels of a 85 Rb atom 
do not anti-cross in a range of magnetic fields up to a few 
hundred G. As for notation, it is therefore often convenient 
to label the Zeeman levels by fm / where / is the total angu- 
lar momentum quantum number at B — with which the level 
adiabatically correlates. We will use such an explicit reference 
to the different Zeeman states in the subsequent subsections . 
An al ternative, brief notation (Mi es et all 120001 IStoof et aH 
1988) indicated in Fig. |7] is to label the Zeeman components 
in order of increasing energy as a,b,c, . . . and shall be used 
throughout this subsection. 



2. Coupled channels scattering 

The hyperfine splitting £hf and the magnetic Zeeman en- 
ergy are typically large compared to the kinetic energy for 
cold collisions. Consequently, the scattering properties of 
two atoms depend strongly on the Zeeman levels in which 



they are prepared. A scattering channel is defined by spec- 
ifying the quantum numbers that describe each of the two 
initially separated atoms. In the following it is convenient 
to perform a partial wave expansion of the stationary scat- 
tering wave function of the relative motion of an atom pair. 
Its components thus depend on the inter-atomic distance, r, 
in addition to the quantum numbers I and mi associated with 
the angular momentum of the relative motion and its orienta- 
tion with respect to the magnetic field axis. Given this con- 
vention, a scattering channel is specified by the channel in- 
dex a = {f\mf l f2inf 2 £m(}. Here the braces signify that the 
wave function is (anti-)symmetric with respect to ex change of 
identical Bose (Fermi) particles (IStoof et all Il988l) . Conse- 
quently, the order of atomic indices 1 and 2 is unimportant. 
Only those channels with even (odd) £ exist for two identical 
bosons (fermions) in the same Zeeman level. All £ values are 
possible for identical bosons or fermions in different Zeeman 
levels or for nonidentical species. The partial waves associ- 
ated with £ = 0, 1, or 2 are designated s, p, or d waves. If 
the channel energy, E a = Ef mf + Ef %mf , is less than the total 
energy of the system, E, the channel a is said to be open. Con- 
versely, if E a > E, channel a is referred to as closed. When 
E is lower than the lowest channel energy, all channels are 
closed and E can refer only to one of the discrete molecular 
levels. If E is higher than the lowest channel energy, then at 
least one channel is open, and E is associated with a stationary 
energy level in the scattering continuum. 

In the case of the 'So 4 He atom of Section |ll]the two-body 
interaction Hamiltonian is specified by a single-channel po- 
tential describing the helium dimer bound state as well as 
cold collisions. Such a treatment is not adequate for alkali 
atoms with nuclear spin because of intrinsic coupling among 
the various channels. It is therefore convenient to represent 
the wave functions for scattering or bound states as a coupled 
channels expansion into their compon ents E) in the sep- 
arated atom spin basis l abelled by a dGaol 1 19961; iMies et all 
1 1996t IStoof et a/.l[l988l) . The associated radial wave function 
F a (r, E) of the relative motion of an atom pair with the energy 
E is determined by the following formula: 



W a (r,E) = F a (r,E)/r. 



(9) 



Substituting Eq. (0 into the stationary Schrodinger equa- 
tion gives the following coupled channels equations which, 
in combination with the boundary conditions imposed on 
F a (r, E), determine both the continuum and bound states: 



d 2 F a (r ,E) 
8r 2 



+ 5 2 [ E6 "P ~ Va P (r) ] F ^ E) = °" (10) 



We note that m should be identified in Eq. ( fTOb with twice 
the reduced mass of an atom pair. The potential matrix, V(r), 
takes on the following form in the separated atom spin basis: 



Ef\m ft + Ef 2 , nf2 + 



tl 1 



■+1) 



8aft + V%(j). (11) 



Here the atomic hyperfine and magnetic interaction terms are 
given by the experimentally known separated atom energies 
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Ef lM/i and Ef im , and the kinetic energy of axis rotation is 
given by the centrifugal energy term proportional to 1(1 + 1). 
These terms are diagonal in the asymptotic basis. The com- 
plicated part of the scattering due to the electronic Born- 
Oppenheimer potentials and electron spin-dependent interac- 
tions is contained in the interaction matrix Vi a t(r) of Eq. (fTTh . 
which is comprised of two parts: 



Vint(r) = Vd(r) + V M (r). 



(12) 



The contribution V e \{r) represents the strong electronic inter- 
action. It is diagonal in ( but non-diagonal in the atomic 
channel quantum numbers /im^/im/,. Its elements diago- 
nal in fimf l f2irif 2 depend on weighted sums of the two Born- 
Oppenheimer potentials Vs = o(r) and V s=1 (r), whereas the off- 
diagonal terms depend on the difference between Vs=o(r) and 
Vs=\(r). The strong electronic interaction, V e \(r), is responsi- 
ble for elastic scattering and inelastic spin-e xchange collisions 
dBenderl fl%3l: iDalgarno and Rudgel 1 1 965l) . and gives rise to 
the broadest scattering resonances. 

The term V ss (r) in Eq. ( fT2l re presents the weak relativistic 
spin-s pin potential energy (Moerdii k et a/.L[T995tlStoof et all 
1988). It is due to the anisotropic dipolar interaction between 
the two electron spins, and is non-diagonal in both fimf l f2inf 1 
and (, i.e., it couples different partial waves. In the limit of 
large inter-atomic distances V ss (r) is proportional to ar 2 /r 3 , 
where a = 1/137.0426 is the fine structure constant. As 
V ss (r) is a tensor operator of rank 2, only blocks that differ in I 
by zero or two units have non-vanishing matrix elements, ac- 
cordi ng to the Wigner-Eckart theorem dEckari fl 9 3 Ot Iwi gner . 
119271) . In addition, there are no s-wave diagonal potentials 
varying as 1/r 3 . 

At short range the spin-dipole interaction can be modi- 
fied by second-order spin-orbit co ntributions, which are im- 
portant for a hea vy atom like Cs (Kotochi gova et all l200lt 
iMies et all 1 1996b . In general, the potential energy contribu- 
tion V ss (r) is responsible for weak inelastic spin-dipolar relax- 
ation and gives rise to narrow scattering resonances. 

The low energy collision physics of alkali atoms in spe- 
cific hyperfine states is sensitive to the re-coupling of electron 
spins between the separated atoms and the short range zone of 
strong chemical interactions. At small inter-atomic distances 
the potential energy scale is orders of magnitude larger than 
£hf ■ Whereas the electron spin is coupled to the nuclear spin 
on the same atom when the atoms are far separated, the elec- 
tron spins become uncoupled from the nuclear spin and cou- 
ple strongly to one another at small r to make the S — and 
1 states of the Born-Oppenheimer potentials. The distance 
where this re-coupling occurs is near r u of Fig. [6] where the 
difference in Born-Oppenheimer potentials Vs=i(r u )—Vs=o(r u ) 
due to the exchange potential is equal to the atomic hyperfine 
energy £),f ■ This occurs typically in the distance range of 20 
to 25 flBohr for alkali atoms. 

The co u pled channels method dMies et all 1 19961; 
IStoofef all (1988) of Eq. ([TOj properly accounts for the 
dynamical changes in the couplings among the five angular 
momenta si, S2, ii, 12, and I as the atoms move through the 
region near r u . Basic symmetries of the coupling terms in 
Vint(r) of Eq. ( fT2b allow us to separate the interaction matrix 



into blocks, within which the coupling is strong and between 
which the coupling is intrinsically weak. Such a separation 
gives rise to classifications of the various stationary energy 
states in terms of their predominant symmetry properties. 
Projected levels that are bound just within a particular block 
are, in general, associated with scattering resonances when 
their energy is above the scattering threshold. We note, 
however, that in contrast to the complete stationary states de- 
termined by Eq. ( [Tol l, such Feshbach resonance levels depend 
on the separation of Vi n t(>") into blocks or, equivalently, on the 
choice of basis set. 

In this context, the separated atom spin basis referred to in 
Eq. ( fTTT i is convenient at long range, but leads to off-diagonal 
elements at short range. One alternative, short range basis 
would first couple Si and S2 to a resultant S and ii and 12 to 
a resultant I. Then S and I can be coupled to a resultant F, 
which in tu rn couples to I to give the total angular m omen- 
tumFtotai dMoerdiik et aZ.lll995l:lTTesinga et aall993l) . Then 
we could set up a molecular basis set with quantum num- 
bers (S I)F£F tot . d iM. Here M refers to the orientation quan- 
tum number associated with the projection of F tot ai onto the 
magnetic field axis. This short range basis takes advantage of 
the fact that V e \(r) is diagonal in S . An alternative separated 
atom basis set could couple f 1 and f2 to a resultant F, and give 
the basis {f\f2)F(F lota \M. This basis is useful at low B fields 
at which the Zeeman levels do not anti-cross, and where F 
may be viewed as a good quantum number. The unitary trans- 
formation between separated atom and molecular basis sets 
is cal l ed a frame transformation dBendeii [19631; iBurke et all 
119981; IDalgarno and Rudgel Il965t iGao et a/.l 120051) . Bound 
and meta-stable states of light element s like Li and Na are 
best classified by the m olecular basis (Moe rdiik et all [19951: 
ISimonucci et al. , 120051) whereas the sepa rated atom basis is 
better for heavy e lements like Rb or Cs dChin et all l2004bt 
iMarte et all 120021) . Similarly, the long-range part of a scat- 
tering wave function is best described in the separated atom 
basis set, whereas the molecular basis is more appropriate for 
the short range part. 



3. Threshold collisions 

The quantum numbers f\m^ and fim,f 2 associated with the 
Zeeman states of the separated atoms in which a dilute gas is 
prepared determine the entrance channel of a two-body col- 
lision. In the context of cold collisions, it is usually suffi- 
cient to consider just the s-wave (6 = 0) component of an ini- 
tial plane wave momentum state of the relative motion of an 
atom pair. To avoid inelastic processes known as spin relax- 
ation, most experimental applications of Feshbach resonances 
involve atom pairs in the lowest energetic Zeeman states for 
which i-wave scattering is allowed. In the following, we as- 
sume such a case and choose the zero of energy at the entrance 
channel energy, E a . This convention implies that the total en- 
ergy E of a colliding atom pair is identical to its positive ki- 
netic energy, H 2 k 2 /m. Here k is referred to as the wave number 
and p = hk is the momentum of the relative motion. On the 
other hand, bound state energies, E\,, are always negative. 
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The amplitudes associated with transitions between the 
initial and final states of a diatomic collision may be in- 
ferred from the asymptotic form of the scattering solutions 
to Eq. ( TTOl l in the limit of large distances, r — > oo. As the only 
open channel is assumed to be the entrance channel, the long 
range boundary condition imposed on the associated compo- 
nent of the radial wave function reads dTavloiill972l) : 



F a {r, h 2 k 2 lm) oc sin [kr + . 



(13) 



Here the absolute magnitude of F a (r,H 2 k 2 /m) is determined 
up to an overall energy dependent pre-factor whose value is 
a matter of convention. All information about a collision is 
contained in the scatt ering phase shift, $jhk). Accord ing to ef- 
fective range theory (Bethe, 1949t lSchwingerlll947ilbh . g(M) 
determines the scattering length, a, via the following low en- 
ergy asymptotic expansion: 



1 1 9 

k cot ij(%k) = — + -k r eff . 

a 2 



(14) 



Here r e a is known as the effective range. For most of the appli- 
cations in this review, the scattering length alone is sufficient 
for the description of cold diatomic collisions. 

In general, the scattering length is only weakly dependent 
on the magnetic field strength, B, unless B can be tuned in 
such a way that a closed channel Feshbach resonance level 
crosses the entrance channel scattering threshold. Such a sce- 
nario may occur due to a difference in magnetic moments, 
dE a /dB and dEp/dB, associated with the entrance and closed 
channels, respectively. The meta-stability of the resonance 
state leads to a time delay during a collision when the energies 
of the scattered atoms and of the resonance level are nearly 
matched. This results in an enhancement of the scattering 
cross section whose width in energy depends on the strength 
of the coupling between the entrance and closed channels via 
the lifetime of the resonance state. As the zero energy cross 
section is proportional to a 2 , such a resonance enhancement 
of collisions may be u sed to widely tune the scattering length 
dTiesinga et all 1993b as illustrated in Fig. [TJ We note, how- 
ever, that because of the inter-channel coupling the magnetic 
field strength Bo associated with the singularity of a differs 
from the crossing point between the resonance energy and the 
scattering threshold. The magnetic field width and shift of 
such zero energy resonances will be the subject of Subsec- 
tion nrc] 

Singularities of the s-wave scattering length are always ac- 
companied by the degenerac y of a bound vibrational level 
with the scattering threshold dTavlorl,[i972l) . In the context of 
magnetic Feshbach resonances, the properties of such a cou- 
pled channels stationary energy state, termed the Feshbach 
molecule, may be inferred from Eq. ( TTOb in the zero bound 
state energy limit, E\> — > 0. Similarly to the studies of Sec- 
tion [II] the derivation just relies upon the fact that the poten- 
tial matrix Vint(r) of Eq. (fTTT i vanishes at large separations, 
r — > oo. In such an asymptotic distance range the atoms cease 
to interact, and the solution of Eq. (TlOb associated with any 
s-wave (( = 0) channel is given by the following formula: 

F a (r, E b ) oc exp [- y/-m(E h - E a )r/h] . (15) 



Its pre-factor depends on the overall normalisation of the 
bound state and is thus determined by all components. In the 
limit E\, — > 0, however, only the entrance channel radial wave 
function acquires a long range and, therefore, predominates 
all the others. The relation between the binding energy and 
the scattering length in the vicinity of a zero energy resonance 
will be the subject of Subsection IIII.DI This discussion will 
show, on the basis of general arguments, that E\, is determined 
by Eq. dTJ, while the Feshbach molecular wave function re- 
duces to its entrance channel component given by Eq. ©. 

If the entrance channel is not the lowest in energy the above 
scenario of zero energy resonances and Feshbach molecules 
is only approximate. This implies that the scattering length 
is always finite even when a closed channel Feshbach reso- 
nance level is magnetically tuned to cross the entrance chan- 
nel energy. The associated Feshbach molecule can decay 
into the lower energetic open channels via spin relaxation 
dKohler et all 120051 ^Thompson et al.i l2005bl) . Such a rather 
general situation underlies, for instance, the experiments of 
Fig. |2] involving gases of 85 Rb, and serves as an example for 
the following explicit channel classification. 



4. Example of channel classification 

The most rigorous classification of scattering channels is 
by the total projection quantum number M = /«/,+/«/, + rti(. 
In the absence of external fields other than B, the symmetry 
with respect to rotation about the magnetic field axis implies 
that M is strictly conserved during the course of a collision, 
i.e., states with different M can not couple. The next useful 
classification is by the partial wave quantum number t. Only 
weak coupling is normally possible between states of different 
t, since it can only originate from the intrinsically small and 
anisotropic V ss (r) matrix elements. 



TABLE I Separated atom quantum numbers for the s-wave (I = 0) 
M = -4 block of the coupling matrix Va(r) for 85 Rb, for which / 
assumes the values 2 and 3. The separated atom energies relative to 
E ee = are shown for B = and 160G (16.0mT). 



Cfi h) 


m h m h 


a 


E a /h [GHz] 


E a /h [GHz] 








f or B = 


forB = 160 G 


(2 2) 


-2 -2 


ee 








(2 3) 


-1 -3 


df 


3.035732 


2.591623 


(2 3) 


-2 -2 


eg 


3.035732 


2.756754 


(3 3) 


-3-1 


fh 


6.071464 


5.508558 


(3 3) 


-2 -2 


gg 


6.071464 


5.513508 



Table U shows an example of the quantum numbers needed 
to describe the s-wave channels associated with the interac- 
tion of a pair of / = 2, m/ = -2 85 Rb atoms, which is the 
e state in the alphabetic notation. This state is one for which 
there is a broad Feshbach resonance close to the scattering 
threshold near 155 G. Both the binding energies of the asso- 
ciated Feshbach molecules and the scattering length may be 
inferred from solutions of Eq. ( fTUl ), using just the matrix el- 
ements of the potential V e \(r) between the channel states of 
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Table U Since the nuclear spin of a 85 Rb atom is i = 5/2, 
the two ground state / values are 2 and 3. The total pro- 
jection quantum number for s waves is M — -4 for any B. 
There are only five possible separated atom spin channels. As 
the hyperfine splitting is E^f/h = 3.035 Ghz for this species, 
cold collisions associated with temperatures T on the order of 
1 fjK (where k^T/h = 21 kHz given the Boltzmann constant 
£b = 1.3806505 x 10~ 23 J/K) have only a single open s-wave 
channel, the ee channel. All the other ( = channels, df, eg, 
fh, and gg are closed. 

TABLE II Separated atom quantum numbers for the open channel 
rf-wave (I = 2) M = -4 block of the coupling matrix V ss (r) for 85 Rb. 
The separated atom energies relative to E„ = are shown for B = 
and 160 G (16.0 mT). There are also 19 closed channels in the d- 
wave M = -4 block. 



(/1/2) 


m/j rrif 2 ni[ 


a 


E a /h [GHz] 


E a /h [GHz] 








f or B = 


forB = 160 G 


(2 2) 


-1 -1 -2 


dd 





-0.161496 


(2 2) 


0-2 -2 


ce 





-0.157321 


(2 2) 


-1 -2-1 


de 





-0.080748 


(2 2) 


-2-2 


ee 









The anisotropic interaction V ss (r) weakly couples the s- 
wave block to the af-wave block, and the c/-wave block to the 
g-wave block, etc. This has two consequences: extra spin- 
relaxation channels are possible, and projected energy states 
of c/-wave character (or even higher partial wave character) 
can give rise to scattering resonances for s-wave collisions. 
In the 85 Rb case there are a total of 23 M = -4 cZ-wave 
channels that couple to the s-wave block illustrated in Ta- 
ble 1] Only the four listed in Table HI] are open with respect 
to the ee channel energy. These channels are all degenerate 
with the ee s-wave channel at B = 0. Such a degeneracy 
at zero magnetic field implies a suppress ion of inelastic col- 
lisions due to the Wigner threshold law dWignerl fl948h out- 
lined in Subsection IIII.DI The of-wave channels of Table HI] 
become open, however, by a relatively large amount of energy 
on the cold /iK temperature scale as B increases. The associ- 
ated energy gap gives rise to inelastic decay by which a pair of 
e-state atoms can relax in a collision when B increases from 
zero dRoberts et a/.[|2000bl) . These open li-wave channels are 
also responsible for the observed spontaneous dissociation of 
Feshbach molec ules with binding energ ies near the ee scatter- 
ing threshold dThompson et all l2005bl) . Associated lifetimes 
as a function of the magnetic field strength will be discussed 
in Subsection llH.EI 

Table Hill illustrates the blocks of the matrix Vmt(r) accord- 
ing to \ f\f2\Fl quantum numbers. The example refers to Bose 
atoms with a nuclear spin quantum number of i = 5/2, such 
as 85 Rb. Basis states in the same vertical column of Table HTT1 
that is, with the same Ft quantum numbers, are coupled by 
the strong exchange interactions in V e i(r). Basis states from 
different vertical columns can only be coupled by the weak 
interactions in V ss (r). Classification using Ft blocks was used 
for Feshbach resonance states with I > obse rved for the 
comparatively heavy alkali metal atoms 87 Rb dMarte et all 



TABLE III Separated atom quantum numbers {fifyFi for the s- 
and rf-wave part of V mi (r) for the bosonic isotope 85 Rb, for which 
/ assumes the values 2 and 3. Odd values of F are missing when 
f\ equals f%, because of Bose symmetry. Spin-exchange interactions 
couple states in the same column, whereas states in different columns 
can only be coupled by spin-dipolar interactions. 



(/. h) 


Ft 


(3 3) 


Qs 


2s 


4s 


6s Od 


2d 


4d 6d 


(2 3) 




Is 2s 


3s 4s 5s 




Id 2d 


3d 4d 5d 


(22) 


0s 


2s 


4s 


Od 


2d 


4d 



120021) and 133 Cs dChin et aZ.ll2004bl) . In particular, the 133 Cs 2 
Feshbach molecules associated with the Stern-Gerlach sepa- 
ration experiments of Fig. [3] are mainly of g-wave character 
with F — 4 and 1-4. Such a separated atom classification 
shall be applied in the following to Feshbach resonance and 
molecular states of 85 Rb associated with the observations of 
Fig. [2] 



5. Near threshold bound states 

The illustration of the relation between resonance and 
bound states on the one hand and singularities of the scat- 
tering length on the other hand relies upon coupled channels 
calculations. Figure [8] shows the five separated atom channel 
energies E a (B) for the s-wave block described in Table U In 
accordance with the convention used throughout this review, 
the scattering threshold associated with the ee entrance chan- 
nel defines the zero of energy. As the five channel energies 
just reflect the sums of their associated atomic Zeeman ener- 
gies, they cluster into three groups. The lowest energy group 
consists just of the ee entrance channel. The next group with 
energies on the order of £hf includes the df and eg channels. 
Finally, the group with energies of about 2£),f refers to the 
fh and gg channels. Both channel energies associated with 
each one of the latter two groups are degenerate at B = 0. At 
the low magnetic field strengths shown in Fig. [8] the choice 
of (fifz)F quantum numbers well characterises the two-body 
physics. The (22) entrance channel has F — 4, the (23) group 
has F = 4 and 5 and the (33) group has F — 4 and 6. Only F 
values of 4 or more are possible because the projection quan- 
tum number is M = -4; the odd value F = 5 is ruled out for 
identical bosons with (/1/2) = (33). 

The solid curves in Figure |8] show the s-wave bound states 
of the 85 Rb2 molecular dimer with negative energies, £\, < 0. 
The curves virtually parallel to the E = threshold refer to 
the v = — 1, -2, and -3 states of the ee entrance channel, la- 
belled by their vibrational quantum numbers, v, starting with 
V = — 1 for the highest excited level. In addition to v, their low 
field (22)4 quantum numbers are indicated in Fig. [8] These 
levels all have nearly the same magnetic moment as a pair of 
separated atoms in the entrance channel Zeeman state config- 
uration. Each one of the other four, closed channels also has 
a vibrational series leading to each of the four closed-channel 
thresholds. The pair of (23)F , v = — 3 levels with F — 4 and 
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FIG. 8 (Colour in online edition) Magnetic field dependent M = -4 
.s-wave energy levels of the 85 Rb2 dimer. The dotted curves labelled 
by ee, df, eg, fh and gg show the energies of the five separated atom 
channels of Table|I] Solid curves indicate the calculated .s-wave cou- 
pled channels bound state energies labelled by the vibrational quan- 
tum number, v. Their symmetry labels refer to the set of quantum 
numbers {fifi)F. The ee limit gives rise to a single vibrational pro- 
gression with F = 4; the df and eg limits give rise to two (23)F series 
with F = 4 and 5; and the fh and gg limits give rise to two (33)F 
series with F = 4 and 6. In all cases, the F = 4 level has the lower 
energy of the (fifz) pairs. At B = the v = -1 and -2, (23)F levels 
and the v = -1, -2, and -3, (33)F levels represent meta-stable states 
with E > embedded in the (22)4 ee scattering continuum. These 
meta-stable levels are represented by dashed curves, which give the 
approximate positions of scattering resonances in the ee channel. The 
meta-stable levels with F = 4 are coupled to the ee s-wave scattering 
continuum through the exchange interaction. The F = 5 and 6 meta- 
stable levels are not exchange coupled to the ee entrance channel at 
B = 0, but become weakly coupled at higher fields where F is no 
longer a good quantum number. 



5 and the pair of (33)F , v = -4 levels with F = 4 and 6 are 
both bound with respect to the separated atoms. For both pairs 
the lowest level has F — 4. These levels have different mag- 
netic moments from a pair of ee separated atoms. This implies 
that their bound state energies, Eb(B), relative to E = vary 
significantly with the magnetic field strength B. 

While the 85 Rb2 molecular dimer states involve the com- 
plete strong electronic interaction, V e j(r), the Feshbach reso- 
nance levels of Fig. [8] refer to blocks of the potential matrix 
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FIG. 9 (Colour in online edition) Scattering length (upper panel) and 
bound state energy levels (lower panel) versus B for two 85 Rb atoms 
in the e state. The separated atom energy of two e-state atoms is 
taken to be zero. A singularity of the scattering length occurs where 
an energy level becomes degenerate with the scattering threshold at 
E = 0. The resolution of the figure is not sufficient to show the 
variation of the bound state energies with B just below threshold; 
this variation is discussed in Fig.[l2]of Subsection llll.DI The broad 
resonance is due to the threshold crossing of the strongly coupled 
(fifi)F£ = (33)4*, v = -3 level, whereas the narrow resonance is 
due to the crossing of the weakly coupled (fifz)F£ = (33)6.?, v = -3 
level. 



associated with their symmetry labels. At zero magnetic field 
the (23), v = — 1 and -2 levels and the (33), v = — 1, -2, and -3 
levels have positive energies so that they are meta-stable lev- 
els embedded in the ee scattering continuum. The two {33)F 
levels are especially interesting, in that they cross the ee E = 
threshold and turn into bound states when B is sufficiently 
large. Figure [9] shows an expanded view of the bound state 
energy levels in the near-threshold crossing region as well as 
the i-wave scattering length a of the ee entrance channel as a 
function of B. The scattering length has a singularity at the 
magnetic field value Bo where a new bound state appears with 
E\,(Bq) = 0. The F = 4 (33) level is coupled by the exchange 
interaction to the F — 4 ee entrance channel, and gives rise 
to a broad resonance with an associated singularity of a near 
155 G. On the other hand, the F — 6 (33) level is only weakly 
coupled to the F — 4 ee entrance channel, and gives rise to 
a narrow zero energy resonance near 220 G. In addition, the 
interaction between the (33)F levels and the highest excited 
entrance channel bound state at E-i/h = -0.22 GHz is evi- 
dent in Fig. [9] It results, for instance, in the avoided crossing 
of the (22)4, v = -1 and (33)4, v = -3 levels between 200 G 
and 250 G which is due to the same coupling that leads to the 
broad singularity of the scattering length near 155 G. 



B. Two-channel approach 

Coupled channels calculations based on realistic molec- 
ular potentials and known atomic properties are capa- 
ble of accounting for a variety of experimental colli- 
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(van Abeelen and Verhaar, 1998; Bartenstein et al, 


2005 


Chin et al. . 


2004b 


: Houbiers et al 


, 1998; Leo era/., 


2000 


Loftus et al. 


, 2002 


; Marcelis et al, 


2004; Marte et ai, 2002) 



Their accuracy, including predictions of new resonances, 
is possible once the actual potentials have been calibrated 
in such a way that they recover the correct scattering 
lengths for the Born-Oppenheimer and potentials 
dAbraham et q/.L Il997t I van Kempen et all 120021) in addition 
to the van der Waals coefficient, C^. Coupled channels ap- 
proaches have the drawback, however, that they are not read- 
ily accessible. Consequently, it is desirable to find simpler ap- 
proaches to Feshbach resonance and molecular levels close to 
the dissociation threshold energy. To this end, the physics of 
binary collisions as well as the properties of the highly excited 
molecular bound states can us ually be well described in terms 



of just two scattering channels (lclrildlll974l;lMies et a/.ll2000l 
iMoerdiik et all 1 19951; iTknmerman s et al. , Il999bl) . In the fol- 



lowing, we denote by the entrance channel the Zeeman state 
configuration of a pair of asymptotically separated atoms in 
which a dilute gas is initially prepared. Under the conditions 
of resonance enhancement, this spin configuration is strongly 
coupled, in general, to several energetically closed scattering 
channels. In idealised treatments, however, this coupling is 
usually due to the near degeneracy of the energies of a sin- 
gle meta-stable vibrational state, the Feshbach resonance state 
<p les (r), and the colliding atoms. The spin configuration asso- 
ciated with the Feshbach resonance state is referred to, in the 
following, simply as the closed channel. 



1 . Two-channel Hamiltonian 

The general Hamiltonian of the relative motio n of an atom 
pair, i.e., the basis of all two-channel approaches dChildl 1 1 974 ; 
Mies et all l2000t iMoerdijk et all 1 1995c iTimmermans et al. 
1999b), is given by the following matrix: 



H2B 



H bg W(r) 
W(r) H C \(B) ) 



(16) 



Its off diagonal elements, W(r), are the energies associated 
with the spin exchange (or dipole) interaction and provide the 
inter-channel coupling as a function of the distance r between 
the atoms. The diagonal elements, H\, g and H C \(B), can be 
interpreted in terms of entrance and closed channel Hamil- 
tonians in the hypothetical absence of coupling, respectively. 
These single channel Hamiltonians consist of kinetic and po- 
tential energy contributions given by the following formulae: 



H bg = --V 2 + V bg (r), 
m 

H d (B) = V 2 + V c i(B,r). 

m 



(17) 
(18) 



Here and in the following we choose the zero of energy as 
the threshold for dissociation of the entrance channel. Conse- 
quently, the background scattering potential Vb g (r) of Eq. (Tl\ 



vanishes in the limit of large inter-atomic separations, in ac- 
cordance with the general asymptotic behaviour of the van der 
Waals interaction described by Eq. ((5]). In this approach the 
complete two-body Hamiltonian of Eq. ( TToT ) depends on the 
magnetic field strength B simply via an overall shift of the 
closed channel potential V c i(B, r) with respect to the entrance 
channel dissociation threshold. The amount of this shift is 
determined by the sum of the single particle Zeeman ener- 
gies associated with the closed channel spin configuration of a 
pair of atoms. Typical diagonal potential energy contributions 
to the two-body Hamiltonian are illustrated schematically in 
Fig. [10] whereas special requirements for them will be dis- 
cussed in Subsection llll.GI 
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FIG. 10 (Colour in online edition) Scheme of the entrance (solid 
curve) and closed channel (dashed curve) potentials associated with 
a model representatio n of the 155 G zero energy resonance of 85 Rb 
<K6hlerera/.Ll2003bh , The horizontal dashed and vertical dot-dashed 
lines indicate the magnetic field dependent energy E mi {B) of the 
Feshbach resonance level <f> rcs (r) and the van der Waals length of 
'vdw = 82 a Bohl , respectively. In this model the outer classical turning 
point of the resonance level if> KS (r) associated with the energy E ICS (B) 
in the closed channel potential is located at about r dassil . a i = 72a Bohr . 
Its vibrational quantum number v = — 1 with respect to the disso- 
ciation threshold of V c i(S, r) is chosen arbitrarily. The horizontal 
solid lines indicate the energies £_i and £L 2 of the bare vibrational 
bound states associated with the entrance channel potential. This 
bare interaction is adjusted in such way that it mimics the highest 
excited (/i/ 2 )F = (22)4, v = -1 and -2 levels of Fig. [8] In this 
model the off-diagonal spin exchange coupling has been chosen to 
be W(r) = Bexp(-r/a), using a = 5«Bohi and B/k % = 38.5 mK to 
match the measured Feshbach resonance parameters (Clauss en et all 
120031) . 



2. The bare Feshbach resonance state 

The closed channel Hamiltonian supports the bare Fesh- 
bach resonance state in accordance with the following 
Schrodinger equation: 



H cl (B)cf, ves (r) = £ res (B)0 res (r). 



(19) 
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The associated resonance energy, E les (B), can be tuned with 
respect to the entrance channel dissociation threshold, i.e. the 
zero of energy in Fig. QJJ by varying the magnetic field 
strength. In the hypothetical absence of inter-channel cou- 
pling the Feshbach resonance level usually represents a com- 
paratively tightly bound state of H C \(B) with a classical ra- 
dius r c ] ass i ca i smaller than the van der Waals length of Eq. (7). 
Under realistic conditions, however, (f> KS (r) is only a meta- 
stable state with a decay width depending on the strength of 
the inter-channel coupling. Its precise form can not always be 
unambiguously identified from full coupled channels calcula- 
tions because the two-channel picture is an idealisation. The 
remarkable accuracy of such a simplifying approach is due 
to the fact that near resonance the physically relevant dressed 
stationary energy levels are largely insensitive to the detailed 
structure of the resonance state. 



C. Dressed energy states 

The dressed stationary energy states consist of two orthog- 
onal components associated with the spin configurations |bg) 
and |cl) of the entrance and closed channels, respectively. 
Their wave functions are ther efore of the following general 
form dDrummond et all Il998l) : |bg)^ bg (r) + |cl)0 cl (r). The 
components (f>b g (r) and <f> c \(r) depend on the relative position 
r of the atoms. In accordance with Eq. ( fToT l, these compo- 
nents are determined by the following set of coupled station- 
ary Schrodinger equations: 

#b g 0bg(r) + W(r)0 cl (r) = E<f> hg (r), (20) 
W(r)0 bg (r) + H cl (B)<f> cl (r) = E<p cl (r). (21) 

All applications in this review will involve either symmetric 
or anti-symmetric spin configurations of identical Bose and 
Fermi atom pairs, respectively. Their physical spatial wave 
functions will only consist of an s-wave component. The so- 
lutions of Eqs. (f20b and ( |2TI ) with a negative energy E be- 
low the entrance channel dissociation threshold are associated 
with molecular bound states. Solutions with a positive energy 
belong to the continuum spectrum of the Hamiltonian and de- 
scribe collisions of initially separated pairs of atoms. 



1 . Dressed continuum states 

The dressed continuum states can be labelled by the rela- 
tive momentum p of a pair of asymptotically separated atoms 
in the entrance channel spin configuration and are associated 
with the collision energies, E = p 2 /m. Due to the continuous 
range of angles of incidence between the atoms, any given col- 
lision energy is infinitely degenerate. In this context it is con- 
venient to postpone the introduction of any symmetry prop- 
erties of the spatial wave functions associated with a possible 
identical nature of the atoms. Continuum states are usually 
chosen in such a way that their entrance channel component, 
0p 8 (r), behaves at asymptotically large inter-atomic distances 
like a superposition of an inciden t plane wave and an outgo- 
ing spherical wave (iTavlori Il972h . This choice of the wave 



function leads to the following boundary condition: 



s (r) 



1 



(2nh) 3 > 2 



Jpr/h 



(22) 



The quantity f(§, p) is known as the scattering amplitude and 
depends on the modulus p of the relative momentum of the 
atoms and on the scattering angle determined by the relation: 
cos $ = p • r/(pr). Due to the off diagonal coupling, W(r), 
the stationary continuum states also have a closed channel 
component. The associated wave function, 0p(r), decays at 
asymptotically large inter-atomic distances because the cold 
collision energies of interest are below the closed channel dis- 
sociation threshold. This property reflects physical intuition 
as closed channel atom pairs are spatially confined by the po- 
tential V c i(r) of Fig. [10] 

Due to their long range of Eq. d22l . stationary continuum 
wave functions may be interpreted in terms of amplitu des for 
the de nsity of particle flux rather than physical states dTavlorl 
1 1972b . In the present context, the entrance channel com- 
ponent, 0p g (r), determines observable low energy scattering 
properties of colliding atom pairs, such as, for instance, the 
differential cross section of distinguishable atoms, |/(#, p)\ 2 . 
These properties can be represented in terms of bare energy 
states associated with the Hamiltonians H bg and H C \(B). To 
this end, it is instructive to reformulate the coupled set of sta- 
tionary Schrodinger equations (201 ) and (2T1 ). including the 
boundary condition of Eq. d22l . in terms of the associated 
Green's functions in addition to the entrance channel contin- 
uum states. The bare Green's functions depend on a complex 
variable z with the dimension of an energy in accordance with 
the following formulae: 



Gbg(z) = (z - H bg ) , 
G cl (B,z)=[z-H a (B)T 1 



(23) 
(24) 



The entrance channel continuum wave functions, p +) (r), also 
referred to as background scattering states, satisfy the station- 
ary Schrodinger equation associated with the Hamiltonian H bg 
ofEq. (Ell, i-e.: 



„2 

H bg <P ( P + \r) = !-<f> { ;\r). 
r m 



(25) 



Their long range behaviour is determined by boundary condi- 
tions analogous to Eq. (22l with f(§, p) replaced by the bare 
amplitude, f bg {§,p), associated with the background scatter- 
ing. 

Expressed in terms of the bare Green's functions and con- 
tinuum states, the Schrodinger equations (l20l i and $2l[ read: 

I0P E > = l<>> + G bg (E + *WI<>, (26) 
|<) = G cl (B,£)W|^ E ). (27) 

The argument "z = E + iO" of the entrance channel Green's 
function in Eq. (26l > indicates that the physical collision en- 
ergy E = p 2 /m is approached from the upper half of the 
complex plane. This choice of the energy argument ensures 
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that the scattering wave function 0p 8 (r) is compatible with 
Eq. d22b . in accordance with the following long range asymp- 
totic behaviour of Gb g (z) in spatial coordinates: 



Gbgfer.r') 



m(2ntif 12 e ipr/n 
Anh 2 r 



(28) 



Here p '(r') = [^^(r')] is the entrance channel continuum 
energy state with in coming spherical wave boundary condi- 
tions (Tavlorl [l972h . Its label, p = (mE) ll2 v I r, may be in- 
terpreted as the asymptotic momentum associated with the 
relative motion of the scattered atoms. The closed channel 
continuum wave function p (r) of Eq. ( |27] i decays at asymp- 
totically large inter-atomic distances because the bare Green's 
function on the right hand side is evaluated at the collision 
energy E = p 2 /m below the dissociation threshold of V c \(r). 

Direct application of the two-channel Hamiltonian (fT~6b to 
Eqs. d26l i and (l27T i verifies their equivalence to the coupled 
Schrodinger equations ( f20b and (f2Tb . The associated set of 
integral equations for the spatial wave functions </> p 8 (r) and 
<pp(r) may be used, for instance, to numerically determine 
the exact dressed ene rgy states in the two-channel approach 
(Kohl er et al .1 l2003bh . This formulation of the two-channel 
scattering problem is also particularly useful for an approxi- 
mate treatment based on the singularities of the closed channel 
Green's function. The single resonance approach, underly- 
ing this treatment, provides analytic formulae for the dressed 
continuum wave functions as well as their associated collision 
parameters in terms of the bare energy states. 



2. Single resonance approach 

In accordance with Eqs. ([T9T l and (1241 . the closed chan- 
nel Green's function has a singularity at the resonance energy 
E K s(B). Provided that the Feshbach resonance state <p les (r) is 
unit normalised, i.e. (0 re sl0res) = 1> the singular diagonal ma- 
trix element of G C \(B, E) is given by the following formula: 



<0 res |G cl (fi,£)|0 res > = [E-E m (B)Y 



(29) 



Based on Eq. (|29l , the single resonance approach takes advan- 
tage of the near degeneracy of the resonance energy with the 
entrance channel dissociation threshold illustrated in Fig. [TO] 
At typical cold collision energies, E = p 2 /m, the resonance 
detuning AE(B) = E — E reti (B) in Eq. ( 1291 is negligible com- 
pared to the spacings between the discrete energy levels of 
H c i(B). Consequently, the closed channel Green's function in 
Eq. (f27) is dominated by its virtually singular diagonal contri- 
bution associated with the Feshbach resonance level. These 
estimates thus lead to t he following approximation dChildl 
ll974tlG6rale7aZl 12004 : 



G ci (B,E)* |0 res ) 



1 



AE(B) 



(<Pn 



(30) 



The single resonance approach consists in replacing the 
bare Green's function in Eq. d2"7T i by the right hand side of 
Eq. (f30b . This replacement implies that the functional form 



of the closed channel component of the dressed continuum 
wave function p (r) is given by the Feshbach resonance state 

4>ies(r)- The associated overlap factor (0 res |0p ) is determined 
by the energy and magnetic field dependent amplitude: 



A{B,E) = (4^\W\<fg)/AE(B). 



(31) 



The dressed state component |0 p ) on the right hand side of 
Eq. ( 1261 can then also be eliminated in favour of the product 
of l^res) and A(B, E). These replacements in Eqs. (1261 and 
(f27) give the following explicit formulae for the dressed con- 
tinuum states in terms of the bare Feshbach resonance and 
background scattering states: 



I0p S > 



|0 p +) ) + GUE + /O) W|0 res >A(fl, E), 



\(/>t) = \cf> Ies )A(B,E). 



(32) 
(33) 



The amplitude A{B, E) of Eq. ( l3TT l may be expressed in terms 
of the same bare states by inserting Eq. (l32l into Eq. Oil . 
This yields: 



A(B,E) 



<0na|W|0 p +) > 



AE(B) - <0 res |WG bg (£ + iO)W\<l>re s ) 



(34) 



We note that the dressed continuum wave functions of 
Eqs. (l32l and ([33l depend on the magnetic field strength only 
through the detuning AE(B) of the resonance energy E les (B) 
in the denominator of Eq. d34l . Within the typical range of 
experimental magnetic field strengths E res (B) is a linear func- 
tion of B to a good approximation. This linear dependence 
may be represented in terms of an expansion of E KS (B) about 
the resonant field, Z? res , at which the Feshbach resonance level 
crosses the dissociation threshold of the entrance channel. In 
accordance with the choice of the zero of energy in Fig. [TO] 
B les is determined b y the relation E^B™ ) = 0. This implies 



the simple formula (Moerdiik et al, 



1995 



£res(S) - Pks(B - Bres)- 



(35) 



Here the slope of the linear curve, yu les , is the difference be- 
tween the magnetic moments of the Feshbach resonance state 
and a pair of asymptotically separated non-interacting atoms. 

The single resonance approach of Eqs. d32l , ( f33l > and ( f34b 
gives an exact representation of the dressed continuum states 
provided that the spatial configuration of a closed channel 
atom pair is restricted to the Feshbach resonance state <p KS (r). 
Such an assumption is associated with the following formal 
replacement of the closed channel Hamiltonian: 



H cl (B) -> |0 res >£ res (B)<0„ 



(36) 



This simplification of the complete two-channel Hamiltonian 
of Eq. ( fT6l ) yields low energy scattering amplitudes, f{§,p), 
sufficiently accurate to determine the magnetic field depen- 
dence of the s-wave scattering length a(B). 



3. Width and shift of a zero energy resonance 

The s-wave scattering length is determined by the long 
range asymptotic behaviour associated with the dressed con- 
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tinuum states in terms of the zero momentum limit of the scat- 
tering amplitude. This lim it is well repr esented by the follow- 
ing partial wave analysis (lTavloilll972h : 



1=0 



p->0 



(37) 



Here £ labels the quantum number associated with the orbital 
angular momentum and P((cos §) is a Legendre polynomial. 
The s-wave scattering amplitude is related to the phase shift of 
Eq. ( fT3l through the formula fo(p) = ft exp [/£(/?)] sin^(p)/p. 
We note that the limit p — > implies rotational symmetry of 
the entire wave function </>p 8 (r), i.e. independence of the scat- 
tering angle §, because the incident plane wave in Eq. (l22l is 
isotropic at zero momentum. The higher angular momentum 
components fdp) of the scattering amplitude, besides the s 
wave associated with t, — 0, are usually negligible in appli- 
cations to binary collisions in cold gases, due to their propor- 
tionality to p 2e . 

The behaviour of 0p 8 (r) at asymptotically large inter-atomic 
distances r can be inferred from Eq. (l28T i and from the bare 
amplitude, fbg(&,p), associated with the background scatter- 
ing wave function, p +) (r). In accordance with the explicit 
representation of the dressed continuum wave function of 
Eq. (l32b . the scattering amplitude is thus given by the follow- 
ing expression: 

m(27r^) 3 <^ ) |lV|0ies) 

m p) = /b g (#, p) - - — 4^2 A ( B < £ )- (38) 

In the zero momentum limit this expression recovers the res- 
onanc e enhanced s-wave s cattering length in terms of the for- 
mula (Moerdijk et al. I ll995h : 



a(B) = a bg 1 - 



AB 
B-B 



(39) 



Its parameters are the background scattering length, abg, the 
width, AB, and the position of the zero energy resonance, 
Bo. The background scattering length is associated with the 
bare scattering amplitude, /b g (#, p), via a relation analogous 
to Eq. (f37j- The zero momentum limit of Eq. (f38T > determines 
the width of the resonance to be: 



AB = 



m{2nti) i 



Ann 2 



UbgPres 



jw\<p ~ 



(+)\ 



(40) 



We note that Eq. ( 1391 predicts the scattering length a{B) to 
assume all values from -oo to +oo due to its singularity at Bo, 
the measurable magnetic field strength at which the zero en- 
ergy resonance occurs. The width, AB, characterises the dis- 
tance in magnetic fields between the position of the singular- 
ity, Bo, and the zero of the scattering length. The denominator 
of the amplitude A(B, E — 0) of Eq. (l34l determines the physi- 
cal resonance position, Bo, relative to the zero energy crossing 
point of the bare Feshbach resonance level, B res , to be: 



The absolute magnitude of Bo enters the two-channel ap- 
proach as an adjustable parameter because the magnetic field 
strength B, es is not directly measurable. The resonance shift, 
Bo - B res , however, determines characteristic properties of the 
highest excited dressed vibrational bound state, such as, e.g., 
the admixture of the closed channel spin configuration to its 
wave function. 



4. Dressed molecular bound states 

In analogy to Eqs. d26i i and (l27T i. the two-channel dressed 
molecular bound states are determined in terms of the bare 
Green's functions by coupled matrix equations of the follow- 
ing form: 



|^ g > = G hg (E b )W\<p c b l ), 
|<^> = G d (B,£ b )W|/ g >. 



(42) 
(43) 



The discrete bound state energies, Eb, are negative, i.e. below 
the dissociation threshold of the entrance channel. Both chan- 
nels are therefore closed and, consequently, the atom pairs are 
confined by the potential wells of Vb g (r) and V c i(r) in Fig.fTOl 

Their associated wave functions, <p^ s (r) and <p^(r), vanish ac- 
cordingly in the limit of large inter-atomic distances, simi- 
larly to the long range asymptotic behaviour of diatomic halo 
molecules of Eq. (fJJ. The components of the dressed two- 
channel bound states can therefore be interpreted in terms of 
probability amplitudes, subject to the normalisation condition: 



<^ E i^ E ) + <^Vb> 



1. 



(44) 



Bo = B res - <0 re s|WGb g (O)W|0 res )//Zn 



(41) 



Figure QT| illustrates the magnetic field dependence of the 
highest excited dressed vibrational bound state, referred to in 
the following as the Feshbach molecule, for the example of 
the 155 G zero energy resonance of 85 Rb. The wave func- 
tions 0^ s (r) and (p^(r) represent exact solutions of Eqs. ( l42l i 
and ( |43T > associated with the model potentials of Fig. [10] us- 
ing /-i res /h = -3.46MHz/G. The entrance channel compo- 
nents <p^ g (v) of Fig. QT| are extended over a wide range of 
inter-atomic distances. This range increases beyond all limits 
as the magnetic field strength B approaches the zero energy 
resonance position of Bo ~ 155 G. The bond length of the 
Feshbach molecules in Fig. QT| i.e. their mean inter-atomic 
distance, is given, for instance, by (r) =521 aBohr at 160G 
and 4255flBohr at 155.5 G. Figure QT| also suggests that the 
functional form of the closed channel compon ent, ^'(r), is 
virtual ly independent of B. Numerical studies dKohler et al.l 
2003b) reveal that (^'(r) is proportional to the bare reso- 
nance wave function, <^ res (?"), which indicates the applica- 
bility of Eq. (13 6t . The probability of detecting an atom 
pair in the closed channel spin configuration, i.e. (0jr |^') = 
j dr ^'(r)! 2 , decreases from about 4.7 % at B = 160 G to only 
0.1 % at B = 155.5 G. Given that $f(r) is spatially confined 
in the same manner as the resonance state, the suppression 
of the closed channel component is a consequence of the in- 
creasing range of ^ 8 (r) and the normalisation condition of 
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FIG. 1 1 (Colour in online edition) Entrance and closed channel com- 
ponents of the highest excited vibrational bound state associated with 
the 155 G zero energy resona nce of 85 Rb versus the inter-atomic dis- 
tance r (Kohlere/a/., 2003b). The wave functions were determined 
using the model illustrated in Fig.[l0]for magnetic field strengths of 
160 G (upper panel) and 155.5 G (lower panel), respectively. The 
nodes of the wave functions at short distances are associated with the 
vibrational states supported by the bare potentials. The long range 
behaviour of the entrance channel wave functions, beyond the van 
der Waals length of 82 OBohr (dot-dashed line), is large ly determined 
just by the binding energy of the Feshbach molecule ( Kohl eref all 
2003a). We note that the inter-atomic distance is given on a logarith- 
mic scale. 



N h = ^1 + (<pr e s\WG 2 bg (E h )W\<p ies ). (46) 

In the single resonance approach all bound state energies E b 
are determined by a constraint on Eq. (|43T > which can be de- 
rived by multiplying Eq. ( I4"2l i by {<f> KS \W from the left. This 
leads to the following formula: 

E b = firesiB - fi res ) + <<MWG bg (£ b )W|0 res >. (47) 

We note that Eq. (l47T i recovers Eq. (PTTT i in the limits E b — * 
and B — > Bo, provided that the magnetic field strength ap- 
proaches Bo from the side of positive scattering lengths. This 
confirms that the binding energy of the Feshbach molecule 
vanishes at the measurable resonance position Bo. Such a 
weak bond implies that the properties of the near resonant 
highest excited dressed vibrational state are determined solely 
by the scattering length, a, in analogy to the general findings 
with respect to halo dimers in Section HI] The range of mag- 
netic field strengths in which the Feshbach molecular state as 
well as its energy depend just on a is usually referred to as the 
universal regime. 



D. Universal properties of Feshbach molecules 

In accordance with the results of Section [Til the properties 
of halo dimers, such as, e.g., the large spatial extent of their 
wave functions in Fig.[TT] can all be inferred from Eq. ([1), the 
universal formula for their binding energy. 



1. Universal binding energy 

The formal deriv ation of Eq. ([T]) for Feshbach molecules 
dGoral et all (2004) relies upon an explicit determination of 
the matrix element involving Gb £ ,(£'b) on the right hand side of 
Eq. (|47| >. This may be performed on the basis of the following 
resolvent identity: 



Gb g (£b) = Gbg(0) - £ , bGb g (0)Gbg(£'b). 



(48) 



Eq. ( |441 . All these trends reflect general properties of Fesh- 
bach molecules which can be readily explained using the sin- 
gle resonance approach. 

In analogy to the treatment of the dressed continuum lev- 
els, the single resonance approach to the two-channel bound 
states is equivalent to the pole approximation to the closed 
channel Green's function of Eq. OOl ). This approximation ren- 
ders Eq. d43l ) into a practical form. Its analytic solution may 
be inserted into Eq. ( I4"2l i which yields the following unit nor- 
malised dressed molecular two-component state: 



(45) 



Here TVb is the associated normalisation constant whose ex- 
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' G bg (E h )W\<p KS ) ' 


{ K> J 




I0res) 



Multiplication of Eq. (|48j by G bg '(0) = -H bg from the left 
and by GzHEb) - E b - H bg from the right readily verifies 
this identity. The contribution to Eq. d47b from the first term, 
Gbg(0), on the right hand side of Eq. (l48b yields the energy 
shift jU res (Bies - Bq) due to Eq. (f4~TT) . To evaluate the contri- 
bution from the second term, -£ , bGb g (0)Gb g (£'b), it is instruc- 
tive to employ the spectral decomposition of the bare entrance 
channel Green's function: 



Gbg(z) 



dp 



Z - p 2 /m 



+ Gl(z). 



(49) 



Here the energy argument z is either zero or E b . The quan- 
tity G^ g (z) includes all contributions to Gb g (z) from the bare 
bound states whose energies, depicted in Fig. [TO] are usu- 
ally far detuned from typical energies of Feshbach molecules 
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in the universal regime. Thus neglecting G? (z) in the prod- 
uct -£ , bGb g (0)Gbg(£'b) determines the zero bound state energy 
limit of Eq. (03 to be: 



fi K ,(B-B Q )-m z E b dp 



I- 



!S IWp +) > 



0. (50) 



p 2 (p 2 + m\E b \) E b 



The corrections neglected on the left hand side of Eq. 
are all linear in E b , while the leading contribution involving 
the momentum integral is proportional to Vl-Ebl in the limit 
Eh — > 0. In particular, the asymptotic behaviour of the inte- 
gral can be determined explicitly via a change of variable in 
spherical coordinates from the modulus of the momentum, p, 
to the dimensionless quantity p/ -\fm\Eb\. This yields: 



dp 



„|W|^ +) >| 2;r 2 |<0 res IW o +) >| 



p 2 (p 2 + m\E b \) E b ^o (m|£ b |) 1/2 



(51) 



Equation d40b may be used to eliminate the matrix element in- 
volving the coupling, W(r), on the right hand side of Eq. ( Bil l 
in favour of the product a bg AB fi KS , The solution of Eq. 
with respect to the bound state energy thus reads: 



tr 



m[-a bg AB/(B - B )] 2 



(52) 



As anticipated (iDonley et al .1 l2002h for the general reasons 
outlined in Section[II] this formula exactly recovers Eq. ([TJ in 
the limit B — > Bo, as the singular contribution to a(B) on the 
right hand side of Eq. (1391 exceeds |flb g |. 



2. The Wigner threshold law 

The preceding derivation of the low binding energy be- 
haviour of the right hand side of Eq. d47l > may be extended 
also to positive continuum energies, E - h 2 k 2 /m, which deter- 
mines the width of the Feshbach resonance level dMies et all 
l2000t iMoerdijk et all 1 19951) . This width in energy is re- 
lated to the decay rate of the bare resonance state into a 
given continuum level associated with the wave number k 
via Fermi's golde n rule (Diirr et al, 2004b; Goral et al, 2004; 
lHaque and Stool l2005l iMukaivama et all 120031) . The ap- 
proach underlying these derivations is based on the spectral 
decomposition of Gb g (z) in Eq. (|49l , in addition to the follow- 
ing formula: 



p 2 /m 



= -md(E-p z /m)+?> 



1 



E — p 2 /m 



(53) 



Here "z — E + z'0" is the complex argument of Gb g (z) intro- 
duced in Eq. ( |26l , and P indicates the principal value of the 
momentum integral in Eq. d49j >. In accordance with Fermi's 
golden rule, the width of the resonance level is given by the 
modulus of the imaginary part of the following matrix element 
of the bare entrance channel Green's function: 



We note, however, that Eq. d53l l is not suitable for determining 
the left hand side of Eq. ( l54l in the case of negative bound 
state energies, z = E b , as the 6 function contribution vanishes 
identically. The derivation of Eq. d52l therefore relies upon 
Eq. g8). 

As the imaginary part of Eq. d54l is associated with the 
decay of the bare resonance state, the product a bg p le;i AB is 
always positive and characterises the strength of the inter- 
channel coupling. Equation (l54l also determines the low en- 
ergy dressed continuum wave functions of Eqs. ( 1321 and ( |33t 
via Eq. ( l34l . In particular, the s-wave scattering amplitude of 
Eqs. ( IT?! and (1381 is given by its well kn own general expan- 
sion about zero momentum (Tavlorl ll97i . 



foffik) 



k^0 



-a(\ - ika), 



(55) 



which is applicable provided that the condition ka <K 1 is ful- 
filled. The proportionality to k — '\fmE/H of the imaginary 
parts on the right hand sides of Eqs. d54b and d55l ) reflects 
a general predi ction by Wigner, known as the threshold law 
(Wigner, 1948). In the context of an analytic continuation to 
imaginary wave numbers, k = i y/m\E b \/ft, Eq. (154-b also yields 
the universal bound stat e energy of Eq. ( T52T ) using Eq. d47b 
(IDuine and Stoofl 120041) . The universal regime may there- 
fore be considered as the range of magnetic field strengths 
in which the Wigner threshold law applies to the energy of 
the Feshbach molecule in Eq. ( |47T >. Its extension about Bo 
depends on the interplay between the spin exchange (or dipo- 
lar) and entrance channel interactions as well as on the prop- 
erties of the Feshbach resonance state (Kohleref al, 2003b; 
iMarcelis et aZl.l2004HMies et all 120001 iNygaard ef »/.U20()6l). 

Figure [12] for instance, illustrates several theoretical ap- 
proaches to the bound state energy of the 85 Rb2 Feshbach 
molecule in the vicinity of the 155 G zero energy resonance. 
Although the universal formula of Eq. ([TJ determines the 
asymptotic behaviour of E b in the limit B — > Bo, it provides 
a reasonable approximation only in a small region between 
Bo ~ 155 G and about 158 G. Both the full coupled channels 
calculation (circles) and the two-channel predictions (solid 
curve) depend on specific properties of the particular 85 Rb 
zero energy resonance besides the scattering length. These 
approaches fully recover the measured Feshbach molecular 
energies (see the inset of Fig.[T2l over the entire experimental 
range of magnetic fiel d strengths from about 156 G to 162 G 
(IClaussen et all 12003b . The slope of the linear magnetic field 
dependence of E b in the limit of high fields in Fig. [T_2]is deter- 
mined by the bare resonance energy, E les (B) (dot-dashed line). 
This indicates an increasing admixture of the closed channel 
spin configuration to the dressed bound state. 



3. Closed channel admixture 

In accordance with Eq. d45l l and the unit normalisation of 
the bare Feshbach resonance state, the closed channel admix- 
ture to the Feshbach molecule is determined by the wave func- 
tion normalisation constant N b of Eq. (|46T i to be: 



<0 res |WGb g (z)W|0 res > ~ p KS (B KS ■ 

k — >0 



b„ 



r 2 dr\<f> KS (r)\ 2 = N h 



(56) 
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FIG. 12 (Colour in online edition) The bound state energy of the 
85 Rb 2 Feshbach molecule as a function of the magnetic field strength 
in the vicinity of the 155 G zero energy r esonance. The circl es in- 
dicate a full coupled channels calculation (Kokkelmans, 200 2]). Th e 
solid curve results from a two-channel approach dGoral et fl/.ll2005l) . 
For comparison, the dotted curve represents the universal estimate 
of Eb, while the dashed curve includes the leading correction to 
Eq. l[l) due t o the van der Waals tail of the background scatter- 
ing potential dGribakin a nd Flambaum. 1993: iKQhler et all l2003ah 
given by Eq. l |66t of Subsection |III.F| We note the pronounced shift 
of about 9G between the measurable resonance position, Bq, and 
the magnetic field strength S lcs at which the resonance energy (dot- 
dashed line) crosses the dissociation threshold of the entrance chan- 
nel. T he inset compares the theoretical approaches with measure- 
ments dClaussen et all 12003) indicated by squares. 



In the context of field theoretic approaches to the many- 
body physics of cold gases, the inverse square of the nor- 
malisation constant, N b 2 , is sometimes refe rred to as the 
wave f unction renormalisa t ion co n stant Z(B) dBraaten et ali 
20031: iBruun and Pethickl 12004 iDuine and Stoofl l2003at 
Goral et all 120041) . The quantity Z(B) can be expressed in 
terms of the difference in magnetic moments of the Feshbach 
molecular state and a pair of separated atoms in the entrance 
channel spin configuration, i.e. dE b /dB, via the following for- 
mula: 



(57) 



This exact result follows directly from Eqs. d46b and d47b us- 
ing the general relation 



G 2 bg (E b ) = -J;G bg (£ b ). 



(58) 



Within the universal regime of magnetic field strengths the 
closed channel admixture to the Feshbach molecule can there- 
fore be inferred from the derivative of Eq. (Q]) with respect to 
B. This yields: 



Z(B) 



2a h /(ma ) 

B-»B fi [es ABa bg ' 



(59) 



Here the limit B — ■> Bo is performed on the side of positive 
scattering lengths of the zero energy resonance. 

In accordance with Eq. j59l , the wave function normali- 
sation constant, N b , diverges as the magnetic field strength 
approaches Bo, i.e. in the limit a — > oo, due to the propor- 
tionality of its leading contribution to ^[a. This implies that 
the closed channel admixture to the Feshbach molecule of 
Eq. d56l i is negligible in the universal regime and vanishes at 
the measurable zero energy resonance position. Within this 
limited range of magnetic field strengths about Bo the Fesh- 
bach molecule can therefore be described in terms of just its 
entrance channel component, in analogy to Section|II] In par- 
ticular, the wave function b s (r) is determined by Eq. (O at 
inter-atomic distances large compared to the van der Waals 
length of Fig. QT| In accordance with Eq. ([3]), the mean dis- 
tance between the atomic constituents of universal Feshbach 
molecules is well estimated by its asymptotic value of one 
half of the scattering length in the limit a — > oo dKohler et all 
2003b). In the universal regime of magnetic field strengths the 
Feshbach molecules are therefore proper diatomic halo states 
with all the general properties described in SectionllTl 



E. Experimental signatures of universality 

The suppression of the closed channel admixture to the 
Feshbach molecular state has been directl y observed for 
the e xample of near resonant 5 Li2 dimers (Partridge et ai, 
120051) . Figure [13] shows a comparison of these measure- 
ments with the universal estimate of Eq. d59l as well as 
with Eq. (|68]» of Subsection IIII.FI The estimate of Eq. d68b 
is based on Eq. ( 1571 ) using the approximation for E b of 
Eq. d66l ), which is il lustrated in the inset of Fig. [12l The un- 
derlying approach (Grib akin and Flambauml 1 1993b accounts 
for the corrections to the universal bound state energy of 
Eq. ([1) due to the van der Waals interaction of the back- 
ground scattering potential in addition to the scattering length 
a(B). Its predictions agree well with exact coupled chan- 
nels calculations of th e closed channe l admi xture to the Fesh- 
bach molecular state (Partridg e et all 120051) . Several experi- 
ments to date have also directly confirmed the universal limit 
of the binding energy, in analogy to the inset of Fig. [T2l 
dClaussen et all 120031 : jPonlev et ali 120021). for a variety of 
alkali atomic species dBa rtenstein et all 120051: iMoritz et all. 
l2005blRegal et adl2003al) ~ 



1. Collisional relaxation 

Besides these demonstrations of universality, the long range 
nature of highly excited Feshbach molecules is manifest in 
their lifetimes with respect to deeply inelastic collisions. 
Such molecular loss may occur due to relaxation into tightly 
bound diat omic states upon collisions with surrounding atoms 
or dimers (Cubizolles et al .1 120031: iMukaivama et all 2003; 
iRegal et all 2004a ). The associated loss mechanism was dis- 
cussed first, in the context of cold gase s, for collisions be- 
tween H2 dimers and hydrogen atoms dBalakrishnan et ali 
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FIG. 13 (Colour in online edition) Occupation of the closed chan- 
nel resonance level relative to the number of atom pairs in a bal- 
anced mixture of two spin components of 6 Li versus the magnetic 
field strength in the v icinity of the 834 G zero energy resonance 
dPartridge et all [2005 ) . On the low field side of the resonance po- 
sition Bo (dot-dashed line) the gas was prepared as a Bose-Einstein 
condensate of 6 Li 2 Feshbach molecules. The experimental data (cir- 
cles) indicate the suppression of the closed channel admixture to the 
dressed molecular state, i.e. Z(B) = N^ 2 , over three orders of mag- 
nitude in agreement with Eq. j68t of Subsection IIII.FI (solid curve). 
The dashed curve shows the universal estimate of Eq. l |59t for com- 
parison. We note that the bond length of the Feshbach molecules of 
Eq. l|3} reaches the order of magnitude of the average inter-atomic 
spacing of the gas of 10 4 a% ohl at a magnetic field strength of about 
800 G. In this strongly interacting regime the purely two-body theo- 
ries are bound to break down. On the high field side of Bo the rem- 
nant small resonance state population indicates pairing phenomena 
in a cold, strongly correla ted two spin compone nt Fermi gas with a 
negative scattering length dPartridge et q/.Ll2005h . 



119971) . In these threshold-less reactions the energy lost 
through de-excitation is transferred to the relative motion of 
the products in accordance with momentum conservation. The 
associated relative velocities are sufficiently high for the scat- 
tered particles to leave an atom trap. The density «d of dimer 
molecules in a homogeneous gas is therefore depleted in ac- 
cordance with the following rate equation: 



«d/«d = -#ad«a - K iA H A . 



(60) 



Here n a is the density of atoms and K aA and denote the 
inelastic loss rate constants associated with atom-dimer and 
dimer-dimer collisions, respectively. 

Several theoretical studies of collisional relaxation of alkali 
metal systems have been performed in the limit of low vibra- 
tional e xcitation of the initial dimer states dCvitas et all 2002 , 



2005allbh as well as for Feshbach molecules dPetrov et al. 
20041 l2005bl) . In the case of collisions between alkali metal 
atoms and their dimers in excited states, whose bond lengths 
are smaller than / v dw, ab initio calculations suggest the in- 
elastic rate constant s , K a j, t o be on the order of 10~ 10 cm 3 /s 
dCvitas et all l2002i 12005 alibi) . For such species the colli- 
sional relaxation rates in cold gases would be too large for 



the observation of any phenomena relying upon the equilibra- 
tion of the molecular component. The expected short life- 
times of alkali dimers have been confirmed, for instance, via 
photo- association of 87 R b atoms in a Bose-Einstein conden- 
sate dWvnaref a/.ll2000l) . 



2. Lifetime of Feshbach molecules in Fermi gases 

Therefore it came as a surprise that cold Feshbach 
molecules produced from incoherent mixtures of two spin 
components of fermionic atoms could be stabilised for up 
to several seconds at densities of about 10 13 atoms/cm 3 
dCubizolles et all 120031: IJochim et all l2003al: iRegal et all 
l2004at Strecker et all 20031) . Under such conditions the ra- 
tio of elastic to inelastic collisions is sufficiently large to al- 
low for an efficient evaporative cooling of the dimers. As a 
consequence of their stability, even the Bose-Einstein con- 
densation of Feshbach molecules has been observed in the 
vicinity of the broad zero energy reson ances of 6 Li and 4Q K 
at about 830 G and 202 G, respective l y dGreiner et all 120031: 
IJochim et al. l l2003btlZwierlein et all 120031) . The associated 
resonance widths, AB, are on the order of 300 G for 6 Li 
(jBartenstein et all 120051) and about 8 G in the case of 40 K 
dGreiner et all 120031) . Such a broadness provides the oppor- 
tunity of performing measurements in the universal regime of 
magnetic field strengths. 

Both species of dimers in these experiments consist of pairs 
of unlike fermionic atoms in the lowest energetic Zeeman 
states. In the case of 4() K the associated atomic spin compo- 
nents are determined by the quantum numbers (/ = 9/2, w/ = 
-9/2) and (/ = 9/2, m/ = -7/2). Here / labels the angu- 
lar momentum quantum number of the hyperfine level with 
which the Zeeman state correlates adiabatically at zero mag- 
netic field, and mj indicates its spin orientation with respect 
to the field axis. Similarly, the pair of atomic Zeeman states 
associated with the 830 G zero energy resonance of 6 Li is de- 
scribed by the quantum numbers (/ = 1/2, my = 1/2) and 
(/ = 1/2, ntf = -1/2). We note that the Pauli exclusion prin- 
ciple allows such atom pairs to interact via s waves given that 
their spin wave functions are anti-symmetric. Consequently, 
all the results of Section [TT] on isotropic halo dimers apply to 
these Feshbach molecules consisting of unlike fermions. All 
experiments on molecular Bose-Einstein condensation have 
been performed at near resonant magnetic field strengths for 
which t he associated dimer bo nd lengths can be as large as 
1 00 nm dZwierlein et al. .1 120031 

The stability of such Feshbach molecules with respect to 
collisional relaxation relies upon the Pauli exclusion princi- 
ple in addition to the separati on of length scales associa ted 
with the initial and final states dPetrov et all |2004 l2005bb . In 
accordance with Eq. the spatial extent of the initial halo 
wave function is determined by the scattering length, a. The 
bond length of the final, deeply bound target level may be 
estimated by the van der Waals length, Z v dw, which is much 
smaller than a. Both dimer-dimer and atom-dimer relaxation 
therefore require at least three fermions to come together at 
short distances on the order of / v dw- One pair among these 
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atoms necessarily shares the same spin state and can inter- 
act at most via p waves. As the momentum scale associated 
with the initial Feshbach molecular state is determined by the 
wave number k ~ I /a, the inelastic loss rate constants are 
suppressed by powers of £/ vc iw ^ 1 ■ 
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FIG. 14 (Colour in online edition) Loss rates of diatomic Feshbach 
molecules (upper panel) and remnant atoms (lower panel) versus the 
magnetic field detuning, B — Bo, in the vicinity of the 202 G zero 
energy resonance in a cold gas of 4Q K with a peak density of about 
1.5 x 10 13 atoms/cm 3 (Ree al et all 12004*3) . The dimers were pro- 
duced using the technique of adiabatic sweeps of the magnetic field 
strength described in Section ITVl Their number is denoted by N A , 
while N a refers to the number of remnant unbound atoms. The circles 
in the upper panel indicate the measured dimer loss, whose general 
trend is well fit by an cT 23 power law as a function of the scatter- 
ing length (solid curve). Significant deviations from this scattering 
length dependence occur in a small region of magnetic field strengths 
where the bond length of the dimers of Eq. {3} is comparable to the 
average distance of the atoms in the gas. The lower panel shows the 
loss rate of remnant unbound atoms in gases with (circles) or without 
(squares) deliberate production of Feshbach molecules. The associ- 
ated trends suggest that the atomic loss is largely unaffected by the 
presence of a dimer component in the gas (Regal et ai, 20043). 

Based on the halo wave function of Eq. (O, the precise a 
dependences of the suppression factors for atom-dimer and 
dimer-dimer relaxati on have been predicted to be a -3 33 and 
cT 2 - 55 , respectively dPetrov et all |2004 l2005bb . While these 
predictions strictly apply just to the universal regime of mag- 
netic field strengths, Fig.[14]illustrates that their general trends 
agree with measurements on c old gases of 4Q K wi th a compo- 
nent of Feshbach molecules (Rega leToZl l2004ah . These ob- 



servations are consistent with the reduction of a predominant 
dimer-dimer relaxation in the limit of large scattering lengths, 
where the lifetimes reach about 100 ms. Conversely, as the 
magnetic field strength is tuned away from the zero energy 
resonance, the molecular lifetimes approach those small val- 
ues of less than a millisecond, typical for short ranged alkali 
dimers in cold gases. 

Besides their scaling properties with respect to the scat- 
tering length, the inelastic loss rate constants of Feshbach 
molecules depend sensitively on inter-atomic interactions at 
short distances below / v( jw- The experimental trends regard- 
ing the collisional relaxation of alkali dimers confined to such 
length scales are inconclusive. Remarkably long lifetimes 
on the order of seconds were reported, for instance, even in 
the case of short ranged, closed-channel dominated Feshbach 
molecules produced in a fermionic 6 Li gas in the vicinity 
of the narrow 543 G zero energy resonance dStrecker et all 
120031) . Such alkali dimers, however, are not described by the 
halo wave function of Eq. (|2]). Their observed stability, there- 
fore, suggests a mechanism for the suppression of collisional 
relaxation beyond the scaling of the associated loss rate con- 
stants with powers of the inverse scattering length. 



3. Lifetime of Feshbach molecules in Bose gases 

Feshbach molecules consisting of identical Bose atoms 
generally tend to be less stable than their fermionic coun- 
terparts. Large collisional relaxation rate constants on the 
typical order of 10~ 10 cm 3 /s have been reported for those 
dimers associated in cold gases in the vicinity of compara- 
tively narro w zero e n ergy resonances of 8 7 Rb, 133 Cs and 23 Na 



dDurr efq/.l l2004at iHerbig et all 20031; iMukaivama et all 
120031: lYurovskv et a/.Ul999bl l200ol). With a width of about 



1 1 G dClaussen et all 120031) the 155 G zero energy resonance 
of 85 Rb is by far the broadest among the bosonic species from 
which Feshbach molecules were produced. The atomic con- 
stituents of these dimers are prepared in the excited Zeeman 
state determined by the quantum numbers (/ = 2, mj = -2) 
which, in contrast to the electronic ground state, can be mag- 
netically trapped. Pairs of such atoms are subject to inelas- 
tic spin relaxation collisions involving the threshold-less tran- 
sition of the atomic energy level to a deeper Zeeman state 
dRobertsef aZ.l fcOQOb). The associated <i-wave exit channels 
are illustrated in Table [II] It is the spontaneous dissocia- 
tion via spin relaxation, rather than collisional relaxation of 
the molecular vibrational state, that predominantly limits the 
observed lifetimes of Feshbach molecules produced in dilute 
gases of 85 Rb. This decay mechanism has been unequivocally 
identified through comparisons with quantita tive predictions 
dKQhler et all 120051: iThompson et adl2005bl) . 

Figure [15] shows the monotonic increase of the measured 
lifetimes over three orders of magnitude as the magnetic field 
strength approaches the resonance position, Bq. The life- 
time of these Feshbach molecules is density independent and 
reaches several tens of milliseconds. This behaviour of the 
molecular decay can be explained in terms of a probabilistic 
estimate based on the average volume, *V = 4^(r 3 )/3, occu- 
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FIG. 15 (Colour in online edition) The lifetime of 85 Rb 2 dimers 
versus magnetic field strength in the vicinity of the 155 G zero en- 
ergy resonance. The circles indicate measurements in a cold thermal 
gas wi th a peak density of 6.6 x 10" atoms/cm 3 (Thompson et ai, 
2005b), while the squares refer to exact coupled channels calcu- 
lations accounting for the sponta neous dissoc i ation of Feshbach 
molecules due to spin relaxation ( Ko hler et al.l 12005). The solid 
curve represents an asymptotic estimate of the lifetime based on the 
universal halo wave function of Eq. ((2} and the spin relaxation loss 
rate constant, K 2 (B), associated with a non-degenerate Bose gas in 
the limit of zero collision energy. The two-body theory fails only 
close to resonance when the bond length of the molecules becomes 
comparable to the average inter-atomic distance of the gas. 



pied by an atom pair in the universal halo state of Eq. d2]i. In 
analogy to Eq. ([3]) this volume is determined in terms of the 
scattering length a by the formula *V = 7ta 3 . The event rate, 
y, associated with the spin relaxation of such a bound atom 
pair may be expressed in terms of the loss rate constant Kj for 
spin re laxation in a non-degenerate Bose gas ( Roberts et al.l 
2000b). This yields y = K 2 /(4'V) dStoof et al.l Il989l) . The 
lifetime, t = 1/y, is therefore proportional to the ratio of 
the magnetic field dependent quantities a 3 (B) and ^(B), i.e., 
t = 4na 3 (B)/K2{B), where, near resonance, KiiB) can be 
evaluated at zero collision energy. The same formula for 
the molecul ar lifetime also fol lows rigorously from Fermi's 
golden rule dKohler et al.l 120051) and exactly recovers the re- 
sults of the full coupled channels calculations of Fig. IT31 in 
the limit B — > Bq. The magnetic field dependence of K2(B) 
can be inferred from general prop erties of inelastic collisions 
dBohn and Juliennel 1 19971 1 19991) . Consequently, the stabil- 
ity of Feshbach molecules in both two spin component Fermi 
gases of 40 K and dilute vapours of 85 Rb directly probes the 
halo nature of these dimers near resonance, despite their dif- 
ferent decay mechanisms. 



F. Classification of zero energy resonances 

The singular behaviour of the scattering length as well as 
the associated existence of a universal regime of magnetic 
field strengths are common to all Feshbach resonance phe- 



nomena in cold gases. Two-body universality implies, in par- 
ticular, that the long range properties of the dressed diatomic 
bound states are well described solely in terms of an effective 
entrance channel interaction. The range of validity of such 
single channel approaches provides classification schemes for 
zero energy resonances. 



1 . Size of the universal regime 

The halo wave function of Eq. (O may be interpreted as 
the only bound state supported by a contact pseudo potential, 
gd(r), describing the low energy spectrum of an atom pair with 
a positive scattering length. Due to their simplicity, contact in- 
teracti ons are widely used i n theoretical descriptions of col d 
gases (iDalfovo et a/.lll999UDalib^d[|l999tlR"anderialll995l) . 
In these approaches the properties of inter-atomic collisions 
enter the many-body Hamiltonian in terms of a coupling con- 
stant, g - 4-nh 2 a/m, depending on the microscopic potential 
just through the scattering length, a. The associated binary s- 
wave scattering amplitude is given by the following formula: 



fo(M) = -a/(l + ika). 



(61) 



We note that Eq. ( l6"TT > also recovers the general low wave 
number expansion of Eq. d55l l in the Wigner threshold law 
regime. Consequently, the contact interaction approach pro- 
vides a minimal implementation of the universal low energy 
two-body physics, applicable to both positive and negative 
scattering lengths. In the latter case the contact pseudo po- 
tential does not support any bound state. 

The assumption of universality of the Feshbach molecu- 
lar state implies the closed channel admixture of Eq. ( TSoT ) 
to be negligible. For magnetic field strengths, B, far inside 
the width of the zero energy resonance, i.e. \B - Bq\ «: AB, 
Eqs. ( T39T > and d59l > therefore yield the following condition nec- 
essary for the applicability of the contact interaction approach: 



B-B 



AB 



2H 2 /(ma 2 bg y 



(62) 



The energy ratio on the right hand side of Eq. ( I62l i provides an 
upper estimate for the extent of the universal regime relative 
to the width of the zero energy resonance. A ratio small com- 
pared to unity indicates a weak coupling between the closed 
and entrance channels. Conversely, broad zero energy reso- 
nances with a large product |jU res AB\ tend to favour negligible 
admixtures of the bare resonance level to the dressed Fesh- 
bach molecular state over a substantial fraction of their width. 
This trend may be enhanced by a small energy H 2 /(ma^ s ) 
which indicates the presence of bound vibrational (ob g > 0) 
or virtual (ob g < 0) near resonant energy levels of the bare 
background scattering potential. 

A small closed channel admixture to the Feshbach molecu- 
lar state, however, is not necessarily identical to the universal- 
ity of the binding energy. While for the example of the com- 
parati vely weakly coupl ed 1007 G zero energy resonance of 
87 Rb dDurr et aZ.ll2004bl) the right hand side of Eq. ([62]) gives 
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0.1, it is as large as 81 in the case of the 155 G zero energy res- 
onance of 85 Rb. The estimate of Eq. d62| i for 85 Rb, however, 
is inaccurate with respect to the extension of the universal 
regime in units of the width, AB, given the results of Fig. [12] 
The inset of Fig.Q~2]reveals that this inaccuracy originates pre- 
dominantly from corrections to the universal binding energy 
due to the van der Waals tail of the bare 85 R b background 
scattering potential ( Gribakin and Flambaumi 119931) . In the 
case of such entrance-channel dominated zero energy reso- 
nances even the lowest order corrections to uni versality may 
be des cribed just by a single channel approach dKohler et al.l 
l2003al) . The minimal requirements on the associated effective 
potential, V{B, r), besides its long range asymptotic behaviour 
of Eq. are closely related to general properties of alkali 
dimer energy wave functions. 



2. Entrance-channel dominated resonances 

In accordance with Fig. QT| the vibrational bound state 
wave functions of 85 Rb2 Feshbach molecules consist of short 
and long range contributions. The characteristic scale for such 
a spatial separation is the van der Waals length of Eq. (0. Due 
to the deep wells of realistic background scattering potentials, 
the entrance channel wave functions of alkali dimers are well 
described by the semi-classical Wentzel-Kramers-Brillouin 
(WKB) approximation at short inter-atomic distances, r <§: 
/ v dw- Their behaviour at large separations, r » / v dw, is de- 
termined mainly by the van der Waals interaction of Eq. 
in addition to the energy of the state. In both spatial re- 
gions the functional forms of the associated asymptotic so- 
lutions ^ftheradial Schrodinger equation are known analyti- 
cally jGaoi Il998allbi: IGribakin and Flambaumi Il993l) . As the 
interactions between alkali atom pairs are dominated by large 
van der Waals coefficients, on the order of thousands of 
atomic units (the atomic unit of C(, is 9.5734xl0~ 26 J nm 6 ), the 
asymptotic wave functions can be matched. Such a matching 
procedure provides the basis of accurate semi-classical treat- 
ments of b ound as well as continu u m entrance channel wave 
funct ions (Flambau m et al.l 1 19991 IGribakin and Flambaumi 
119931) . 

Provided that the interaction of a pair of alkali atoms 
is well described by an effective entrance channel poten- 
tial, V(B, r), the semi-classical approach to the zero en- 
ergy wave function determines the scattering length to be 
( Grib akin and F l amb aum| 1 1993b : 



a = a [1 - tan((/) W KB - tt/8)] . 



(63) 



Here a is usually referred to as the mean scattering length and 
(fv/KB is the semi-classical phase shift. In accordance with 
the WKB approach, ^>wkb consists of the following integral 
between the zero energy classical turning point, ro, associated 
with the effective potential, V(B, r), and infinite distances: 



VWKB 



i r 



dr ^-mV{B, r). 



(64) 



coefficient a of Eq. (l63l depends just on the van der Waals 
length, Z vc j W . Its explicit expression in terms of / v( j W and Eu- 
ler's F function reads: 



a = 



/vdw r(3/4) 
V2 1X5/4) 



0.95598 x L 



(65) 



In accordance with Eq. d63l . the parameter a determines the 
characteristic scale of the scattering length. This average po- 
tential range is modulated by the poles of the tangent function 
provided that its argument, <£>wkb - tt/8, is close to an odd in- 
teger multiple of n/2. As i^wkb increases, each singularity of 
a indicates the emergence of an additional vibrational bound 
state in the potential well. 

The functional form of Eq. d63l using a realistic interaction 
is largely analogous to the formula for the scattering length of 
the simplified sq uare well plus hard core model of an inter- 
atomic potential (IGribakin and Flambaumlll993l) . In the case 
of this exactly solvable model the mean scattering length, a, 
recovers the finite outer radius of the well. The analogy be- 
tween these realistic and simplified interactions therefore im- 
plies that in the limit a » a the energy of the highest ex- 
cited vibrational state is well approximated by the following 
asymptotic formula: 



E\, * -H 2 / |m(a - a) 2 ] . 



(66) 



While (fiwKB is thus sensitive to the entire well of the interac- 
tion, from V(B, ro) = to its long range tail of Eq. 0, the 



An independent analysis based on effective range theory 
for realistic inter-atomic interactions dFlambaum et a/.lll999t 
lGaolll998al) confirms the magnitude of the range param eter a 
of Eq. ( l66l ) up to a constant factor on the order of unity 
l2004t) . The remnant uncertainties associated with Eq. d66l ) 
may be related to the divergence of the effective range ex- 
pansion for an y potential with a long range van der Waals tail 
(lTavloii[T972h . Figure [T2l illustrates the accuracy of Eq. ( 1661 ) 
within the experimentally relevant range of binding energies 
of near resonant 85 Rb2 Feshbach molecules. 

According to the semi-classical approach, at magnetic field 
strengths in the vicinity of entrance-channel dominated zero 
energy resonances the diatomic bound and continuum spec- 
tra are determined mainly by a(B) and a. In such a case, any 
potential, V(B, r), which at each magnetic field strength, B, 
accounts for the scattering length as well as the van der Waals 
tail of Eq. (0 provides a suitable description of the near res- 
onant binary physics. This conclusion is rather intuitive be- 
cause cold collisions are characterised by de Broglie wave 
lengths too large to resolve details of the effective entrance 
channel interaction besides its long range behaviour. The ad- 
justment of V(B, r) to recover the exact magnetic field depen- 
dence of the scattering length of Eq. ( f39b may be achieved, for 
instance, by varying the radius of its hard core. In the context 
of quantum defect theory such a description of alkali dimer 
spectra in terms of just the parameters a and Cf, has been rig- 
orously deri ved for a ran ge of energies much wider than the 
cold regime lGaolll998al) . 

As a is positive, the semi-classical estimate of the bound 
state energy of Eq. d66l ) is always below the universal predic- 
tion of Eq. (0 and becomes singular in the limit a(B) — > a, 
outside the range of validity of Eq. ( |66| |. This unphysical 
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behaviour is counterbalanced in Eq. d47l i by an increasing 
closed channel admixture to the Feshbach molecule, in ac- 
cordance with Eq. (l59l l. which tends to impose a linear slope 
on E\,(B). The principal question of the applicability of sin- 
gle channel approaches outside the universal regime of mag- 
netic field strengths is d ecided by wh i ch on e of these trends 
prevails near resonance (iKohler et al.l [2004). Consequently, 
the bound state energy, E\,(B), in the vicinity of an entrance- 
channel dominated zero energy resonance is subject to the fol- 
lowing inequality: 



E b (B) + 



m[a(B) — aY 



E b (B) + 



H 2 



ma 2 (B) 



(67) 



Within the range of validity of Eq. ( |66l ) the admixture of the 
closed channel resonance state to the Feshbach molecule is 
small compared to unity and can be determined from the bind- 
ing energy via Eq. (157) . This yields: 



Z(B) = 



2ah 2 /(ma 2 ) (1 - a hg /a) 2 
fi KS ABa bg (1 -a/a) 3 



(68) 



The accuracy of Eq. d68l ) in applications to 6 Li2 Feshbach 
molecules in the vicinity of the entrance-channel dominated 
834 G zero energy resonance is illustrated in Fig. [13] 

Based on Eq. d67| > a more practical criterion can be de- 
rived from a low energy expansion of the right hand side 
of Eq. (14Tb using a spe cific implemen tation of the general 
two-channel approach (Gor al et al. 1 120041) . This yields a di- 
mensionless parameter, 77, whose smallness indicates the va- 
lidity of Eq. d66|) beyond the un i versal regime of magnetic 
field strengths ( Stol l and Kohleii [2005). Consequently, an 
entrance-channel dominated zero energy resonance fulfils the 
condition 



a H /(ma ) 
77 = <sc 1. 



(69) 



Such a criterion also results from a n adiabatic description of 
Feshbach resonances dPetrovll2004 . 

In the opposite limit, 77 » 1, a zero energy resonance 
shall be referred to, in the following, as closed-channel domi- 
nated. Closed-channel dominated zero energy resonances are 
typically narrow and their universal regime of magnetic field 
strengths is experimentally largely inaccessible. The descrip- 
tion of their physical properties therefore crucially relies upon 
an explicit treatment of at least two scattering channels. 



G. Characteristic parameters of zero energy resonances 

For any two-channel approach to be sensible its implemen- 
tation should describe the two-body energy spectrum beyond 
the Wigner threshold law domain. Otherwise the same physics 
could be captured simply by the contact pseudo interaction of 
Subsection lHI.FI The approach should therefore recover both 
the scattering length of Eq. d39l and the binding energy of 
the Feshbach molecule beyond the universal regime. Even 
within such a comparatively wide range of energies a vari- 
ety of inter-atomic potentials are capable of describing the 



same two-body physics. It is the objective of this subsec- 
tion to provide a minimal set of physical parameters that ev- 
ery two-channel approach should account for, and to illustrate 
a practical implementation for a typical experimental situa- 
tion. The adjustment of the Hamiltonian of Eq. ( TToT ) in the 
single resonance approach of Eq. (|36] l will be performed on 
the basis of the general form of its energy spectrum derived in 
Subsection IIH. CI According to these derivations, the dressed 
two-channel energy states depend on the bare states associ- 
ated with the background scattering via G bg (z), on the reso- 
nance energy, E res (B), and on the product W\<p KS ) characteris- 
ing the inter-channel coupling. Just these quantities need to be 
adjusted. The specific form of the resonance wave function, 
<p res (r), does not affect the two-body spectrum in the single 
resonance approach. 



1 . The background scattering potential 

In accordance with Subsection lHI.FI all implementations of 
Vhg(r) which recover a bg in addition to Eq. (|5]l yield equivalent 
energy spectra beyond the cold regime dGao|]l998al) . provided 
that their number of levels is large compared to unity. A hard 
sphere in addition to the van der Waals tail of Eq. ©, for 
instance, provides such a background scattering potential with 
a minimal number of parameters. Its explicit expression reads: 



VUr) = 



-C 6 /r 6 



r < r 
r > r 



(70) 



The associated backgrou nd scattering length is given by the 
following exact formula dGribakin and Flambaumi[l993h : 



flbg = a[l - tan(<^wKB - 37r/8)] ■ 



(71) 



Here a is the mean scattering length of Eq. d65l ) and <^>wkb = 
2l 2 m /r^ is the semi-classical phase shift of Eq. d64| >. The dif- 
ference of 7r/4 between the arguments of the tangent functions 
in Eqs. d63l and ( TtTT > is due to the discontinuity of Eq. ( T70l at 
the core radius ro. Equation ( TTTb may be used to determine 
ro in such a way that V bg (r) of Eq. d70l ) exactly recovers the 
background scattering length. Similar procedures have been 
performed in two-chan nel approaches using continuous im- 
plementations of Vbs(r) (Kohler et al. , 2003b; Mar celis et all 
|2004 iMies et all l2000fc iNvpaard et all |2006() . An example 
of such an effective interaction is illustrated in Fig.[l0]in addi- 
tion to its energy levels. Table|IV|provides the parameters flbg 
and C(, characterising the background scattering potential for 
several experimentally relevant zero energy resonances. 

Equation (TTOb as well as those equivalent implementa- 
tions of Vbg(r) that explicitly include the long range asymp- 
totic van der Waals interaction of Eq. ((5]) are, in principle, 
suited to describe several bare energy levels. The number 
of bound states supported by Eq. (TTOt , for instance, can be 
arbitrarily increased by decreasing the core radius ro under 
the constraint of a fixed background scattering length. The 
description of such wide energy ranges, however, is usually 
beyond the scope of two-channel approaches. In addition, 
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TABLE IV Parameters characterising the background scattering potential associated with several experimentally relevant zero energy reso- 
nances. Those values of a\, g and E-i that are unreferenced refer to calculations performed for this review. The energy, £Li, associated with the 
highest excited vibrational state of the background scattering potential is given only in the cases of isolated resonances. The atomic unit of the 
van der Waals dispersion coefficient, Ce, is 9.5734 x 1CT 26 Jnm 6 , «Bohr = 0.052917 nm is the Bohr radius and 1 G=10~ 4 T. 



Species 



So [G] 



«ba [«Bohr] 



C 6 [a.u.] 



\E-i\/h [MHz] 



"Li 

23 Na 
40 K 

85 Rb 
87 Rb 
133 Cs 



543.25(5) (Streckeref al. . 2003) 
834.149 (Bartenstein et al. 2005) 
853 ( Stengeref q/.. 1999) 
907 (Stengeref al. 1999) 
202.10(7) (Regale? al. . 2004c) 
224.21(5) (Regal and Jin. 2003) 
155.0 (Thompson et al. 2005a) 
1007.40(4) (Volz etal. 2003) 
19.90(3) (Chin et al. 2004b) 
47.97(3) (Chin et al. 2004b) 
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-1405 ( Bartenstein et al. 2005) 
63.9 (Mies et al. 2000) 
62.8 (Mies et al. . 2000) 
174(7) (Loftus etal. 2002) 
174(7) (Loftus etal. . 2002) 
-443(3) (Claussen efaL 2003) 
100.5 (Volz etal. 2003) 
163 (Julienne et al. 2004) 
905 (Julienne et al. 2004) 



1393.39 (Ywetal. 1996) 
1393.39 (Ym etal.. 1996) 

1561 (Kharchenko efql.. 1997) 208 

1561 (KharchenkoefqL 1997) 218 

3897 (Derevianko etal. 1999) 8.6 

3897 (Derevianko etal. 1999) 8.6 

4703 (van Kempen et al. 2002) 218 

4703 (van Kempen et al. 2002) 24.0 

6890(35) (Leo etal. 2000) 

6890(35) (Leoef al. .2000) 0.045 



the explicit treatment of the van der Waals tail of the back- 
ground scattering potential is largely impractical in applica- 
tions to the many-body physics of dilute gases. Several im- 
plementations of two-channel approaches, therefore, use ef- 
fective low energy interactions to recover different aspects 
of the cold collision physics under the condit i ons of reso 
nance enhancement dBruun and Pethick, 2004; Chin, 12005 



Goral et al. 



Drummond and Kheruntsvan, 2004; Duine and Stoof, 2004 
"2004H Kokkelman s et al.U2002b . 



Figure [16] shows a typical range of energies relevant to 
the Stern-Gerlach separation of Feshbach molecu les from a 
87 Rb Bose-Einstein condensate dDurr et all l2004ah . This ex- 
perimental technique based on exposing a mixtu re of atoms 
and dimers to an inhomo geneous mag netic field (IChin et all 
120051: iDiirr et all l2004at Herbig et all 120031) is illustrated in 
Fig. [3] for the example of 133 Cs2- The relative force between 
the two components of the gas is proportional to the field 
gradient as well as to the difference in magnetic moments 
of a Feshbach molecule and a pair of separated atoms, i.e., 
dEb/dB. In the rubidium case the magnetic field dependence 
of E\,(B) in Fig. Q~6] is sensitive to the avoided crossing of 
dressed energy levels of Rb2 due to the highest excited bare 
vibrational level of Vb g {r). In accordance with the size of 
the positive back ground scattering length of about lOOflBohr 
(IVolz et al. I l2003h . this level is sufficiently close to the disso- 
ciation thres hold that its energy , E-\, is directly probed by the 
experiment 

An effective entrance channel interaction suitable for a de- 
scription of this experiment needs to recover both the precise 
value of flbg and the bare vibrational level with the energy E-\, 
beyond the Wigner threshold law regime. The associated min- 
imal bare Green's function, Gb g (z), is therefore required to re- 
produce the bound state pole at the energy argument z = E-\. 
Similarly to Eq. fl30l ), this requirement lead s to the following 
separable representation dGoral et q/.ll2004l) : 

Gbg(z) = Gofe) + G (z)Urbg>Tbg(z)<ttg|G (z). (72) 

Here the form factor, \xbg), needs to recover, via the relation 
Go(-E-i)i^bg) K \<P-i), the bare vibrational state satisfying the 



Schrodinger equation H\,o\<p-i) — Zi-il^-i). The resonance 
term, Tb g (z), may be represented, by the following ratio: 



Tbg(z) = 



1 - ^bg(^bg|Go(z)|^bg> 



(73) 



As Tbg(z) is required to reproduce the singularity of the bare 
Green's function in the limit z — > E~\, the amplitude, £bg, is 

given by l/<#bg|Go(£-i)ltf>g>. 

The separable representation of Gb g (z) of Eq. d72b provides 

the exact Green's function associated with an effective poten- 
tial, 



< = Ittg&gCttgl. 



(74) 



determined by the amplitude, £b g , and the form facto r, |yh g )- 
This follows directly from the resolvent identity dTavlorl 
119721) . i.e.. 



Gbg(z) = G (z) + G (z)V b e g ff G bg (z) 



(75) 



which is readily verified upon multiplication by G b g(z) = 

z - #bg from the right and by GZ (z) = z + S 2 V 2 jm from 
the left. Iterating Eq. d75l l yields the Born series which re- 
duces to a geometric series in the special case of the separable 
potential of Eq. d74l i. Its exact sum is given by Eq. d72l i with 
the resonance term of Eq. ( 1731 ). 

Such effective interactions ar e commonly employed in the 
theory of few-part icle systems dLovelacei fl964t Mitral [i962t 
Yamaguchi, 1 19541) as well as in condensed matter physics 
dSchriefferlll964 . In applications to cold gases the form fac- 
tor is not resolved because it is sensitive only to the physics on 
distance scales on the order of the van der Waals length. Its 
associated functional form is therefor e arbitrary and ma y be 
chosen, for instance, to be Gaussian dGoral et all 120041) . In 
the momentum space representation this yields: 



ipkbg) = xbg(p) = 



exp 



(-*) 



{2nhf' 2 



(76) 
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FIG. 16 (Colour in online edition) Avoided crossing of the highest 
excited vibrational levels of 87 Rb 2 (upper panel) and the magnetic 
moment of the Feshbach molecule (lower panel) versus the magnetic 
field strength in the vicinity of the 1007 G zero energy resonance. 
The solid curve of the upper panel indicates the bound state energy 
E b (B) of the Feshbach molecule, while the dashed curve refers to the 
next more tightly bound dressed vibrational level. The dotted and 
dot-dashed lines are associated with the energies E^/h = -24 MHz 
of the bare highest excited vibrational level of Vb g (r) and E ICS (B) 
of the closed channel resonance state, respectively. The crossing 
between the bare levels at 1001.7 G leads to the measured varia- 
tion of the magnetic moment of the Feshbach mo lecule indicated by 
the circles in the lower panel toiirr et all l20 04a). For comparison, 
the so lid and dashed curves refer to ex act coupled channels calcula- 
tions (Ivan Kempen and VerhaaJ. [2004) and a two-channel approach 
jG6ralefqZ.ll2004l) , respectively. 



Here <7b g accounts for the range of the interaction, and (r|p) = 
exp(/p ■ r/K)/(2nK) 3 ^ 2 denotes the plane wave associated with 
the relative momentum p. Given this choice of form factor, 
the parameters cr bg and £bg are determined by the requirement 
that the bare Green's function of Eq. ( 1721 exactly recovers a\, g 
and E-\. The zero energy limit of the scattering amplitude as- 
sociated with Eq. d72l . i.e., /b g (0) = -raTb g (0)/(47rft 2 ), yields 
the condition 



1 + x/ V^' 



(77) 



where x = m^b g /(47r/i 2 cr bg ) is a dimensionless variable. In 
addition, the bare energy level determined by the pole of the 



resonance term of Eq. d73l gives the following relation: 
1 - 4= f V^ye^erfcCy) - ll = 0. 



(78) 



9 poo 2 

Here erfc(y) = -^= I e~" du is the complementary error func- 
tion with the argument y = ■\/m\E-i\crb g /h. 

For the example of the 1007 G zero energy resonance of 
87 Rb the energy £_i can be determined using the potential 
of Eq. ( 1701 1. To this end, its parameter Tq should be cho- 
sen in such a way that the number of vibrational levels is 
large compared to unity and that the known quantities a\, g 
and C(, of Table ITVl are recovered. This yields E-i/h = 
-24 MHz in agreement with c oupled channels calculations 
(Ivan Kempen and Verhaaii 120041) . Based on the precise val- 
ues of £_i and flb e , Eqs. (|77| > and dTSJ, in turn, give the range 
parameter and the amplitude of the separable potential to be 
°"b g = 44a B ohr and m^ g /(4-nh 2 ) = -339«Bohr, respectively. 

We note that the above adjustment of the effective entrance 
channel interaction is restricted to zero energy resonances 
with a positive background scattering length, i.e., «bg > 0. In 
the opposite case, a\, g < 0, separable interactions do not sup- 
port any bound state, similarly to the universal contact pseudo 
potential of Subsection IIII.FI Among the experimentally rel- 
evant examples of Table HV1 only the broad, entrance-channel 
dominated zero energy resonances of 6 Li and 85 Rb have a neg- 
ative background scattering length. Associated two-channel 
approaches describing the dressed energy levels beyond the 
Wigner threshold law regime exist, at least, in applications to 
the broad resonance of 8S RbdG6ralefa/.U 2005 ; Kohler et a 



12004 iKokkelmans and Holland |2002|) . In the case of 6 Li 
a double resonance approach suitable fo r the description of 
many -body systems has been suggested ( Kokk eimans et al\ 
120021) . 

Figures[T2land[T3lreveal, however, that in both cases, 85 Rb 
and 6 Li, the Feshbach molecule is well described by an effec- 
tive single channel interaction, over a wide range of magnetic 
field strengths, in accordance with Subsection IIII.FI An ap- 
propriate magnetic field dependent single channel separable 
potential, suitable for a pplications to few-body bound states 
(Stoll and Kohler, 2005) as well as the dynamics of cold gases 
(IKohler et all l2003al) . may be constructed on the basis of 
Eqs. (1741 and d76l . The associated constant range parameter 
of the form factor is given by <x bg « The adjustment 

of £b g via Eq. (l77l to the magnetic field dependent scattering 
length a{B) of Eq. fl39l l instead of just «bg ensures, in turn, 
that the near resonant energy E\,{B) of the Feshbach molecule 
recovers Eq. 



2. The resonance energy 

In accordance with the single resonance approach of 
Eq. d36l >. the difference in energies of the resonance level and 
the entrance channel dissociation threshold, E res (B), charac- 
terises the closed channel part of the Hamiltonian. The as- 
sociated relative magnetic moment, yK res = dE res /dB, may be 
inferred from a Stern-Gerlach separation experiment, relying 
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upon the force experienced by an atom due to the inhomoge- 
neous magnetic field, B. Given any definite Zeeman state, 
this force is of the general form F a = -V£ a where E & is 
the Zeeman energy. The magnitude of F a = -(<9.E a /<9.B)V|B| 
depends on the orientation quantum number rtif of the to- 
tal angular momentum with respect to the field direction. 
This dependence is w ell described by the Breit-Rabi formula 
dBreit and Rabilll93ll) : 



F a = ±- 



2m f 



2(l + ^ X + X > 



1/2 



■^BohrV|B|. (79) 



Here i is the quantum number associated with the nuclear spin, 
//Bohr = 9.27400949 x 10~ 24 J/T denotes the Bohr magneton, 
and the Lande factor, gj * 2, refers to the electronic magnetic 
moment. The dimensionless variable, x = gjftBotxW/Ehi, 
depends on the field strength, |B|, in addition to the hyper- 
fine energy split t ing in the absence of magnetic fields, 
(lArimondo et all 1 19771) . For any given orientation quantum 
number, ntf, and atomic species, the sign of the force of 
Eq. (l79l is determined by the Zeeman multiplet. 
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FIG. 17 (Colour in online edition) The Zeeman multiplets of 87 Rb 
associated with the total angular momentum quantum numbers / = 
1 (solid curves) and / = 2 (dashed curves). The inset shows the 
magnetic moments, dE^/dB, determined by Eq. ([79} and indicates 
the orientation quantum numbers ntf of those Zeeman states that are 
relevant to Fig. 1161 

For the example illustrated in Fig. [17] the nuclear spin 



quantum number i = 3/2 of 87 Rb ( Ari mondo et all Fl977h in 
addition to the electronic spin give rise to two Zeeman mul- 
tiplets. These states adiabatically correlate with the hyperfine 
levels of total angular momentum quantum numbers / = 1 
and / = 2 in the limit of zero magnetic field. The s-wave en- 
trance channel spin configuration of any pair of the identical 
Bose atoms in the experiments of Fig. [16] is characterised by 
the quantum numbers (/i = I, rn.fi = +l;fi = l,*"/-, = +1), 
referring to their electronic ground states. In accordance with 
Eq. ( |79l and the zer o field hyperfin e struc ture splitting of 
Etf/h = 6834 MHz dArimondo et all 1 19771) . the total mag- 
netic moment associated with these ground state atoms at 
Bo = 1007.4 G is determined by 28EJ8B = -h x 2MHz/G. 
This derivative varies by less than 1 % over the range of mag- 
netic field strengths displayed in Fig. [TBI 

The predominant spin exchange interaction couples the 
pairs of 87 Rb ground state atoms to four s-wave scat- 
tering channels characterised by pairs of Zeeman levels, 
(fi, rn.fi; fz, Mfi), whose total angular momentum orienta- 
tion quantum number is conserved, i.e., tn.fi + rn.fi = +2. 
These closed channels are therefore described by (1, 1; 2, 1), 
(1,0; 2, 2), (2, 0; 2, 2) and (2, 1; 2, 1). In accordance with the 
inset of Fig. [17] their associated total magnetic moments 
at resonance amount to 0, ft X 0.9MHz/G, h x 1.9MHz/G, 
and h x 2MHz/G, respectively. The magnetic moment of 
the 87 Rb2 Feshbach molecule in the lower panel of Fig. [16] 
consists, in principle, of contributions from all five chan- 
nels wei ghted by t h eir a dmixtures. Both the experimen- 
tal data (jDurr et all 2004a) a nd the co upled channels pre- 
dictions ( Ivan Kempen and Verhaan, 120041) . however, are con- 
sistent with a Feshbach resonance state predominantly con- 
sisting of the pairs of Zeeman states just from the / = 2 
multiplet. The Breit-Rabi formula therefore yields an esti- 
mated difference in the magnetic moments of the resonance 
state and the entrance channel spin configuration of about 
ju res = h x 4 MHz/G in the vicinity of the zero energy reso- 
nance position. In an analogous manner the size of the param- 
eter ju res may be inferred from Eq. d79l ) for a variety of species. 
This procedure is particularly useful in the case of genuinely 
two-channel problems, such as, for instance, the spin config- 
urations relevant to the 20 2 G zero energy resonance of 40 K 
dBruun and Pethicfl [2004). The recommended values of ju les 
of Table [V] are all based on coupled channels predictions. 



3. Inter-channel coupling 

The inter-channel coupling leads to the decay width and 
energy shift of the bare resonance level due to its interaction 
with the background scattering continuum and entrance chan- 
nel bare vibrational states. In accordance with Eq. ( 1541 ). both 
quantities are determined by the parameters AB of Eq. ( f40b 
and Bo - B res of Eq. d4TT > in the Wigner threshold law domain. 
While the width in the magnetic field strength, AB, is mea- 
sured routinely, the associated shift is not directly observable. 
Its magnitude, however, may be inferred using general ideas 
dev eloped in the context of multichannel quantum defect the - 
ory dJulienne and Miesl 1 1989t [Mies and RaouIH l2000l |2004h . 
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TABLE V Parameters characterising the inter-channel coupling associated with several experimentally relevant zero energy resonances. Those 
values of AB and /j les that are unreferenced refer to calculations performed for this review. In accordance with Eq. l |69l >, the size of the parameter 
77 indicates whether a zero energy resonance is closed- or entrance-channel dominated. 



Species 


£>a \\>\ 




a ill riv/rtTv/m 


i] 


Li 


34j.2D(5) (Mrecker et al., 20U3) 


A 1 

u. 1 


2.5 


1213 




834.149 (Bartenstein et al, 2005) 


-300 (Bartenstein et al. 2005) 


2.8 


0.02 


23 Na 


853 (StenserefaL 1999) 


0.01 (Mies efa/..2000) 


5.24 (Mies et al. 2000) 


1090 




907 (SteneerefaL 1999) 


1.0 (Mies et al. 2000) 


5.24 (Mies efqZ..2000) 


11 


40 K 


202.10(7) (Reeal etal. 2004c) 


7.8(6) (Greineref al. 2003) 


2.35 (Nvsaarde?aL2006) 


0.46 




224.21(5) (Reeal and Jin. 2003) 


9.7(6) (Reeal and Jin. 2003) 


2.35 (Nveaardef al .2006) 


0.37 


85 Rb 


155.0 (Thompson et al. 2005a) 


10.71(2) (Claussen et al. 2003) 


-3.26 (Kokkelmans. 2002) 


0.04 


87 Rb 


1007.40(4) (Volz et al. 2003) 


0.21 (Diirr etal. 2004b) 


4.2 (Diirr ef al. 2004a) 


5.9 


133 Cs 


19.90(3) (Chin et al. .2004b) 


0.005 (Julienne et al. 2004) 


0.798 (Julienne et al. 2004) 


437 




47.97(3) (Chin et al. 2004b) 


0.15 (Julienne et al. 2004) 


2.09 (Julienne et al. 2004) 


0.99 



This yields the following approximate formula: 



a 



1 - flbg/fl 
1 + (l - abg/a) 



(80) 



Here a is the mean scattering length of Eq. ((65). Conse- 
quently, the size of the resonance shift depends just on the 
quantities abg, AB, and Cc which are all accessible to experi- 
mental studies. Entrance-channel dominated zero energy res- 
onances, such as the example of 85 Rb illustrated in Fig. [12J 
tend to have large shifts comparable to the size of AB. The 
predicted magnitude of Bq - B les = 9 G in Fig. fl2]is consis- 
tent w ith the coupled channels binding energies ( Kokk elmansl 
120021) . In the case of the closed-channel dominated 1007 G 
zero energy resonance of 87 Rb in the upper panel of Fig . [Tolthe 
value of Bq - B res - 0.07 G of Eq. d65t is significantly smaller 
than the width. For all the examples given in TableslTVland[Vl 
the corrections to Eq. ( f80b are negligible. 

The quantity W\(p KS ), characterising the inter-channel cou- 
pling in the single resonance approach, should be adjusted in 
such a way that the matrix element of the bare Green's func- 
tion in Eq. ((54) recovers both AB and Bo - B res . Such an 
adjustment ensures, in particular, that the bound state energy 
of Eq. d47b properly interpolates between the universal and 
asymptotic regimes of magnetic field strengths. An explicit 
minimal implementation of the two-channel Hamiltonian of 
Eq. ( fT6l > may therefore be based, for instance, on the separa- 
ble background scattering potential of Eq. ( f74b in addition to 
the following general expression: 



W\<t> KS ) = WH- 



(81) 



Here the amplitude £ determines the overall magnitude of the 
inter-channel coupling and \x) its functional form. As the off 
diagonal potential, W(r), is not resolved by the typically large 
de Broglie wave lengths associated with cold collisions, the 
form factor \x) may be chosen as a Gaussian fun ction in mo- 
mentum space, similarly to Eq. ((76). This yields dGoral etal\ . 



2004): 



exp 



<pUr> =x(p) 



(-#) 



(2ttH) 3 / 2 



(82) 



The associated range parameter, <r, and the amplitude, are 
completely determined by the requirement that the imaginary 
and real parts of Eq. ((54) recover the physical quantities AB 
and Bo - B res , respectively. For the purpose of this adjustment, 
it is convenient to introduce the average range parameter cr 



-^(<t 2 + crf^) 1 / 2 . The imaginary part of the matrix element 



V2 V 

of the bare Green's function on the left hand side of Eq. (1541 ) 
yields the first condition 



AB 



m\(\ 2 



4nh 2 a bg iJ. le 



(83) 



The associated widths of several experimentally relevant zero 
energy resonances are summarised in Table IVl Given the bare 
Green's function of Eq. d72) and Gaussian form factors, also 
the real part of Eq. ((54) can be evaluated analytically. This 
leads to the second condition, 



Bo - B rl , s — AB 



1 _ *a-(e:\ 2 

«bg 1 yfer W 
\ VW cr J 



~\Jno 



(84) 



whose value for the resonance shift on the left hand side is 
provided by Eq. ( f80b . Equations d83b and ((84), in turn, simul- 
taneously determine the parameters cr and £ characterising the 
inter-channel coupling. We note that the overall phase of the 
amplitude £ is irrelevant to the physics described by the asso- 
ciated Hamiltonian. 

For the example of the 1007G zero energy resonance 
of 87 Rb such an adjustment yields cr = 22«Bohr and 
m\(\ 2 l(Anfi 2 cr) — h x 10 MHz. Given the effective background 
scattering potential of Eq. ( 174) and yU res of TableM this proce- 
dure provides a complete implementation of the two-channel 
single resonance approach. Its predictions with respect to 
dressed binding energies are illustrated in the upper panel of 
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Fig. Q2] as well as in Fig. Q~2] The lower panel of Fig. [16] 
shows the magnetic moment of 87 Rb Feshbach molecules de- 
termined from the product /i les Z(B) via Eqs. (l46l l and (|57| | 
using the two-channel single resonance approach. The over- 
all offset of the data of about -2 MHz refers to the magnetic 
moment associated with the entrance channel spi n configura- 
tion. The comparisons with the measurements dDurr et all 
12004a) and with the results of coup led channels calcula- 
tions (van Kempen and Verhaarl [20041) indicate that the effec- 
tive two-channel approach can fully recover the microscopic 
physics within the experimental energy range. 



IV. ASSOCIATION OF FESHBACH MOLECULES 

The different experimental techniques for molecular asso- 
ciation in cold gases all depend in one way or another on 
the properties of the diatomic energy spectra. Several ap- 
proaches to the production of cold Feshbach molecules are 
based on t he relaxation of an a t omic gas into d i mers n ear 
resonance (Jochi m et al 1 l2003at IZwierlein et al. 1 120031) or 
on dynamical sweeps of the magnetic field s t rength across 



2003 



Regal et all |2003a; 



£n (ICubizolles et all 120031; [Purr et all l2004al: iHerbig et al. 



Strec ker et all I2003L IXu et all 



2003b . Both techniques take advantage of the degeneracy of 
the Feshbach molecular energy, Et,, and the threshold for dis- 
sociation into free atoms in the limit B — > Bo. Relaxation 
of an atomic gas into dimer molecules requires collisions of 
at least three atoms to balance the energies and is commonly 
employed, to date, just in two spin component mixtures of 6 Li 
Fermi gases. The conceptually simpler molecular association 
via magnetic field sweeps seems to be more generally appli- 
cable to both Bose and Fermi gases and will therefore be the 
main subject of this section. Figure [TBI illustrates its princi- 
ple which relies, for the example of the 1007 G zero energy 
resonance of 87 Rb, upon the adiabatic transition from the di- 
atomic zero energy continuum level to the bound state energy 
E\,(B) with decreasing B. Conversely, Fig. [12] suggests that 
in the case of the 155 G zero energy resonance of 85 Rb free 
atom pairs may be associated to Feshbach molecules by in- 
creasing B across Bo- The difference in energy of the colliding 
atoms and diatomic molecules is absorbed by the time depen- 
dent magnetic field. As, in general, the field is spatially ho- 
mogeneous, the association process does not affect the centre 
of mass momentum of the atom pairs. From this viewpoint, 
the Feshbach molecules produced by m agnetic field s weeps 
or other related dynami cal techniques (iDonlev et al.l 120021; 
[Thompson et all l2005al) are as cold as the atomic gas they 
originate from. 



A. Linear sweeps of the magnetic field strength 

In an idealised treatment of the molecular association the 
magnetic field strength may be assumed to vary linearly in 
time. This implies the relation 



where B is usually referred to as the ramp speed and f res is 
chosen to be the time at which the bare resonance energy, E res , 
crosses the dissociation threshold of the entrance channel. The 
field strength, Z? res , associated with f, es is indicated in Fig. [12] 
In accordance with Eq. Q51 >. the resonance energy is then also 
a linear function of time, 



^res(0 — ^res \t ^res)? 



(86) 



with a constant derivative, E res = iu KS B. This formula pre- 
supposes that the magnetic field sweep is sufficiently slow for 
the electronic degrees of freedom of the atoms to adiabati- 
cally adjust to the changes in the magnetic field strength. To 
achieve mol ecular association of an atom pair, E m% needs to 
be negative jMies et al.l [2000). Consequently, the Feshbach 
resonance level is swept downward in time. This requirement 
determines B through the sign of yu les . Conversely, upward 
sweeps of the Feshbach resonance level across Bo lead to heat- 
ing of the atomic cloud. This general principle may be readily 
verified on the basis of the magnetic field dependence of the 
discrete spectrum of dressed energy levels of a trapped atom 
pair. 



1 . Adiabatic association of Feshbach molecules 

Tight traps containing just a single pair of atoms may be 
reali sed, for instance, by the individual sites of an optical lat- 
tice (Jaksch et al, 2002; Tiesinga et al, 2000) or by micro- 



fabricated materia ls ( Folman et all I2002E iHinds et al.l 



Long et al. ,l2003t Muller et q/.[|l999HThvwissen et al.l 



2001 



1999; 



Weinstein and Libbrechtl 119951) . The two-body energy spec- 
tra associated with periodic potentials of optical lattices have 
been studied theoretica l ly in the Wigner thresho ld law do- 
main (lOrso et al. 1 120051: IWouters and Orsol 12 006) as well as 
observed experimentally dMoritz et a~. , 20051) . While in the 
limit of high excitations the tunnelling of atoms is significant, 
the deepest localised diatomic levels of tight lattice sites are 
well described by the harmonic oscillator approximation. As 
a consequence, the centre of mass and relative motions of an 
atom pair confined to a single site can be treated separately, 
similarly to the two-body problem in free space. 

A spherically symmetric harmonic confinement modifies 
the bare entrance channel Hamiltonian associated with the res- 
onance enhanced interaction as follows: 



h 2 , 

#b g = V 2 + V bg (r) + VtrapCr). 

m 



(87) 



Here V trd? (r) denotes the potential energy of the isotropic har- 
monic oscillator in the barycentric frame which is given in 
terms of the reduced mass, m/2, and the angular frequency, 
Who, to be: 



1 m 



2 „2 



VtrapM = 



(88) 



Similarly to Eq. ( [25] ). the associated bare vibrational energy 
levels are determined by the stationary Schrodinger equation, 



B(t) = B les + B(t-t les ) 



(85) 



Hbg<f>v( r ) = E v <p v (r). 



(89) 
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Here the index v = . . . , -2, -1, 0, 1,2,... labels the vibra- 
tional excitation in such a way that v = correlates adiabati- 
cally, in the limit o>h — > 0, with the dissociation threshold of 
the entrance channel. According to this counting scheme, neg- 
ative indices, v < 0, are associated with the bare vibrational 
levels E-i, E-2, etc., indicated in Fig.[l0] In the limit of low 
vibrational excitation, v > 0, the spatial extents of the bare 
states, including <po{r), are characterised by the trap length, 
Oho = V^/( nl<y ho)- This length scale usually very much ex- 
ceeds the modulus of the background scattering length, |ab g |. 
To first order in abg/«ho the energies of the e xcited levels 
are w ell approximated by the following formula dBusch et all 
1998): 



E v 



l +2v+ Jafv + 1/2W 
2 \ tt\ v )a ho 



(90) 



Here (' + ,'/ 2 ) is a combinatorial. The dressed energy levels may 
be determined via the two-channel Hamiltonian of Eg. dT6| > 
using the single resonance approach of Eq. (f36b (iMies et all 
|2000|) . or via a single channel energy dependent cont act in- 
teraction dBlume and Greenel 120021 Bolda e7aZll2002l) . Both 
methods yie ld spectra consistent with full coupled channels 
calculations dTiesinga et adl20"ool) . 
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FIG. 18 (Colour in online edition) Schematic illustration of the 
molecular association of a pair of ground state 87 Rb atoms via a 
downward magnetic field sweep in a spherical harmonic atom trap 
with an oscillator frequency of v ho = 39 kHz ( Thalhamm er et all 
2006). The bare vibrational levels (v = 0, . . . , 6) associated with 
the background scattering quasi continuum and the Feshbach reso- 
nance energy, £ rcs (B), are indicated by dotted and dashed lines, re- 
spectively. Solid curves refer to the magnetic field dependence of 
the dressed energy levels in the vicinity of the zero energy resonance 
position of B = 1007.4 G. 

Figure [18] shows the magnetic field dependence of the 
dressed two-body energy levels of ground state 87 Rb atoms 
predicted by the two-channel approach fo r a tight trap with a 
frequ ency of Vho = u>h /(2n) = 39 kHz dThalhammer et all 
2006). An atom pair from a rubidium Bose-Einstein conden- 
sate prepared on the high field side of the 1007 G zero energy 
resonance and loaded adiabatically into an optical lattice site 
is well described by the v = state. According to Fig. [TjO 



this level adiabatically correlates with the Feshbach molecule 
on the low field side of Bo- Consequently, a magnetic down- 
ward sweep across Bq associates the separated rubidium atoms 
to molecules with certainty in the limit of zero ramp speed, 
i.e. B — > 0. An excited trap level with the vibrational quantum 
number v is transferred to v - 1 . Conversely, an adiabatic up- 
ward magnetic field sweep across Bq dissociates 87 Rt>2 Fesh- 
bach molecules and leads to heating transitions from v to v + 1 
in the excited trap levels. In general, the ramp direction for 
such cooling or heating transitions is determined just by the 
sign of ju res or, equivalently, by the time variation of the reso- 
nance energy, E res . 



2. Exact time evolution of a single atom pair 

While the dressed energy levels of Fig.[l8]reveal the mech- 
anism of Feshbach molecular association in the limit of zero 
ramp speed, the dynamics of the diatomic wave function is 
described by the Schrodinger equation, 



(91) 



Here H2B is the two-channel Hamiltonian of Eq. ( [Tol l in the 
single resonance approach of Eq. d36b . Its time dependence 
is determined by the linear variation of the resonance energy 
of Eq. (f86b . The diatomic state, IT'(O), has components in the 
entrance and closed channels, whose wave functions may be 
expanded into bare states, in accordance with the formulae 



V 

¥ d (r, t) = 4> KS (r)C KS (t). 



(92) 
(93) 



Such a basis set expansion of Eq. d9"TT l leads to the following 
dynamical equations for the associated time dependent coeffi- 
cients 

mC v (t) = E r C r (t) + <0v W res >C res (0, (94) 
ihC KS (t) = £ res (r)C res (f) + 2vfeWv>C v (0. (95) 

V 

This configuration interaction approach is particularly useful 
for numerical treatments of Feshbach molecular association 
using any form o f time dependent magnetic field variation 
dMies et q/.ll2000l) . The special case of linear sweeps belongs 
to those quantum mechanical p roblems whose exact dynam- 
ics can be treated analyticall y dDemkov and Osherovl 1 19681: 
iMacek and Cavagnerolll998l) . 

To this end, it is useful to split Ho^it) into stationary and 
time dependent contributions: 



dyr 



HmW = fttat + ffci(r). 



(96) 



Here the stationary Hamiltonian, // stat = H2-R(t K& ), is associ- 
ated with the magnetic field strength B les . In the single reso- 
nance approach the dynamical contribution, H c \{t), is given by 
the following formula: 



H c i(f) = |0 res ,cl>£ res (f)<0 res ,cl|. 



(97) 
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Similarly to the effective potential of Eq. (l74l i. the time depen- 
dent interaction, H c \{t), involves only a single state associated 
with the closed channel, the product |</> res ,cl) = |^ les )|cl). The 
following derivations will show that it is the separable form of 
Eq. Wl\ in addition to the linear dependence of the resonance 
energy fires on t which allow for the analytic treatment of the 
time evolution. 

The complete dynamics of an atom pair exposed to a linear 
magnetic field sweep may be inferred from the time evolution 
operator determined by the Schrodinger equation, 



ih—U 2B (t,t') 
ot 



H 2B (t)U 2B (t,t'), 



(98) 



in addition to the boundary condition U2B(t,t) = 1. The di- 
atomic state at time t is thus given in terms of U 2B {t, ?') and 
the state at time f by 



mt)) = u 2 B(tj')mt')). 



(99) 



Similarly, the dynamics in the absence of the interaction of 
Eq. Wl\ is described by exp[-/// stat (f - t')/h]. Associated 
with these free and complete time evolution operators are the 
retarded Green's functions, 



(100) 



(101) 



where 0(t — f') is the step function which yields unity if t > t' 
and zero elsewhere. The Schrodinger equation ( f98l > may be 
represented in terms of Eqs. ( 1100b and ( 1101b via the following 
formula: 

G«(r,0=G«(r-*') 

+ J dTG ( £ t {t-T)H d {T)G 2 + B \T,t'). (102) 



Differentiation with respect to the variable t readily verifies 
this relation by recovering the time derivative of the complete 
retarded Green's function of Eq. ( 1101b which is directly de- 
termined by Eq. d98l ). The integral representation of Eq. d98l l 
chosen in Eq. ( 11021 ) is usually referred to as the post-form of 
the dynamical equation. The associated pre-form reads 

+ f dTG%(t,T)H d (T)G ( ^(T-t') (103) 



and can be verified similarly to the derivation of Eq. ( 11021 ). 

While, in general, the operator equation ( 1102b requires a 
complete basis set expansion for its numerical solution, the 
separable form of Eq. ( f97T > reduces this problem to the calcu- 
lation of just the pair of matrix elements 



- O = cl|G™(f - Olfcc., Cl> 

j^fcO = <0 res ,cl|G«(f,f')l0 res ,cl>. 



(104) 
(105) 



This follows from Eq. ( f9Tb via multiplying Eq. ( 1102b by 
(0 res , cl| from the left and by |0 res , cl) from the right which, 



in turn, determines Eq. ( 11051 ) through the following integral 
equation: 



g 2 % f)=g£ t (t-f) 



.(+)/ 



(106) 



Given the free retarded Green's function of Eq. ( 1100b , the 
complete time evolution operator can be inferred from the so- 
lution of Eq. ( 11061 ) by inserting Eq. ( 1103b into Eq. ( 1102b and 
performing the time integrations. 

Using the convolution theorem, a Fourier transform turns 
the integral on the right hand side of Eq. ( 1106b into a product 
of functions. In the case of a linear magnetic field sweep such 
a procedure allows Eq. ( 11061 ) to be solved analytically. To this 
end, it is instructive to introduce the energy dependent matrix 
elements, 



gstat(z) = J 



(107) 



(108) 



associated with the free and complete retarded Green's func- 
tions. Here the regularised argument "z — E + z'0" ensures 
the convergence of the time integrals in the limit t — > +oo 
by approaching the real energy E from the upper half of the 
complex plane. Given the resonance energy of Eq. d86*l ), a 
Fourier transform renders Eq. ( 11061 ) into the following inho- 
mogeneous first order linear differential equation: 

.dg 2B (z,t') 



iti- 



3E 



= hp'( Z ,t')g2B(z,t') + l/E n 



(109) 



Its dependence on the inter-atomic interaction is incorporated 
in the energy derivative <p'(z, f) = d(f(z, t')/dE of the complex 
phase 



<Piz,t') 



1 



dE' 



+ E(t' - t tes )l%, 



(110) 



hE Tes Jo gstat(-Z') 

where 'V = E' + iO" denotes a regularised integration variable 
with the same imaginary part as z. 

The imaginary part of <p(z, t') can be inferred from the 
Hamiltonian H stat using a spectral decomposition of g st m(z) 
analogous to Eq. d49l ). In particular, Eq. d53l determines the 
sign of Im<£>(z, ?') to be: 



sign[Im^(z, f')] = -sign(£ , res )sign(Imz). 



(HI) 



Consequently, in the case of a downward ramp of the Fesh- 
bach resonance level, i.e. E les < 0, the damped, retarded solu- 
tion to Eq. ( 11091 ) is given by 



g2B(Z, t') 



Je 



dE' 



-i[tp(z,f)-<p(z'f)] 



(112) 



This formula may be verified from Eq. ( 11091 ) by differentia- 
tion with respect to E. The matrix element of the complete re- 
tarded Green's function of Eq. ( 11051 ) is just the inverse Fourier 
transform of Eq. dl 12b which yields: 

M , r dE r°° , e -««')-^-'/)] 

^-Js-J. (113 » 
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The exact time evolution operator determined by Eq. ( II 13b is 
applicable to the association of Feshbach molecules in free 
space as well as to the case of a trapped atom pair illustrated 
inFig.[H 



3. Landau-Zener approach 

Already in 1932 Landau and Zener have independently de- 
rived the foundation for a simple estimate of the probabil- 
ity for molecu l e production in linear magnetic field sweeps 
<lLandadfl9l2tlZenerl [l932). Their generic approaches may 
be interpreted in terms of a coupled system of two channels, 
each of which containing just a single state. In applications to 
Feshbach molecular association such a treatment is equivalent 
to the single resonance approach of Eq. ( f36b in addition to the 
following replacement of the entrance channel Hamiltonian: 



H\> g -> |0o>£o<<N- 



(114) 



Here |0o) may be interpreted, for instance, in terms of the 
zeroth vibrational state of the relative motion of a trapped 
atom pair with the energy Eq. This reduction of the two- 
channel continuum to a two-level system gives rise to analytic 
solutions of the coupled set of stationary Schrodinger equa- 
tions ( T20t and d2TT l. Given the simple form of the entrance 
channel Green's function in the Landau-Zener approach, i.e. 



Gbg(z) = I0o> — ^TT<M 
z-E 



(115) 



Eq. ( |47l ) determines the two-level dressed energies to be: 



E + = 



Eq + E KS \Eq - E n 



1+4 KWo)P 

(E () - £ res ) 2 



The magnetic field dependent slopes of the levels E + and 
£_ indicate, quite generally, a crossing of the discrete bare 
entrance- and closed-channel energies, Eq and E le;i , respec- 
tively. Analytic representations of the associated dressed two- 
component stationary states, \<p+) and \<p-), can be inferred 
from Eq. (l45i 

In accordance with Eq. (II 16b . a downward sweep of E KS 
across Eq transfers a pair of atoms in the initial state |0o,bg) 
into the final state |0 res ,cl) in the limit of zero ramp speed. 
Their energies adiabatically follow the E- curve. Given the 
linear variation of E les of Eq. (186b , such a simplified sce- 
nario presupposes the sweep to start and end asymptotically 
far from the crossing point of the bare levels. This, in turn, im- 
plies the formal limits t\ — > -oo and ff — > oo of the initial and 
final times, respectively. The adiabatic energy variation of the 
Landau-Zener two-level approach is similar to the harmonic 
trap case illustrated in Fig. [18] except that the quasi contin- 
uum of excited levels with the indices v > is neglected. 

Finite ramp speeds allow an atom pair to end up in a super- 
position of the entrance and closed channel bare states |0q, bg) 
and |0 res , cl), respectively. In accordance with Eq. (|1051 l and 
the assumption of a two-level system, the probability for an 



asymptotic transition between the bare levels is given by 

|2 



fO.res = |<0 res ,cl|t/ 2 B(ff,fi)l<A),bg> 
= 1 -|*'%2B(ff,fi)| 2 



(117) 



in the limits t[ — > -oo and ff — > oo. These time limits can be 
determined analytically using the stationary phase condition 
for the energy integrals over the oscillatory functions on the 
right hand side of Eq. ( II 13b . This exact approach relies upon 
the observation that asymptotically only those regions close 
to zeros of the derivatives of the phases tp'(z, ?i) and <p'(z,tf) 
contribute to the integrals. All the remaining energy ranges in 
which the complex exponentials are rapidly oscillating yield 
negligible averages. 

The positions of the stationary phases can be readily found 
in the case of a two-level system. An explicit determination 
of the stationary Green's function, G sta t(z) = (z — # s tat) as- 
sociated with the Hamiltonian H stat of Eq. (196b yields 

1 /gstat(z) = z - (cf> res \WG hg (z)W\<p ies ). (118) 



This formula and Eq. ( IllOl i give the derivative of the phase to 
be 

„ ^ E - £ res (f) - {<f> KS \WG hg (z)W\cf> res ) 
<p(z,i) = . (119) 

Consequently, ip'iz, f) vanishes at the dressed energies E ± as- 
sociated with the magnetic field strength at time f, in accor- 
dance with Eq. (1471 1, 

As only those regions of parameters E and E' in the close 
proximity of E ± significantly contribute to Eq. ( II 131 ), the 
phases may be expanded to second order about their stationary 
points. This yields: 

tfz, t) * <p( z± , t) + \v"(E ± , t)(E - E ± f. (120) 

Here z± — E ± + iO denotes the regularised energy parameter 
associated with E ± at time f. The second derivative, 

<p"(E ± , t) = - JL [l + {<f> KS \WG 2 bg (E ± )W\<p KS )] , (121) 

is always positive in the case of a downward ramp of EresW 
and its inverse plays the role of a variance in the complex 
Gaussian integrals over E and E' . Performing the Gaussian in- 
tegration associated with f; on the right hand side of Eq. dl 131 ) 

gives 

r°° „ e WzL.ft)+»r/4] 

dE' e l{fiz •'' ) * V27r ) OjEl - E), (122) 

Je yJ<p"{Elti) 

Here z± and E\. refer to the dressed energies at the initial mag- 
netic field strength. Their asymptotic values in the case of a 
downward sweep of E KS are given by E\ — > oo and E[_ — > Eq 
in the limit f; — > -oo. This implies ip"{E\,t\) — > -l/(ftE les ) 
and <^"(£l,fi) — > oo. Consequently, only the term associ- 
ated with E\ significantly contributes to the sum of Eq. ( 11221 ). 
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The integral over the parameter E on the right hand side of 
Eq. dl 131 ) can be evaluated similarly to Eq. ( 1122b . This yields 



ihg^(t t , to 



e 



-i[<K£t,«f)-to(£U)] 



(123) 



in the limits t\ — > -oo and ff — > oo. Here E f _ refers to the 
dressed energy with the asymptotic behaviour E- — > -oo at 
the final time of the magnetic field sweep. The imaginary parts 
of the phases <p(E_, ff) and (p(E\, t\) relevant to the transition 
probability of Eq. ( II 17b may be obtained from Eqs. (Ill 01 > and 
This leads to the Landau-Zener formulae, 



PO.res 
P0,0 



1 -e 

-2jr<5 L z 



-2n5u 



(124) 
(125) 



Here pq q denotes the probability for detecting the atom pair in 
the bare state |0o,bg) at the end of the magnetic field sweep, 
and 



K0res|W|0 O )| 2 



(126) 



is the Landau-Zener parameter. 

For a pair of atoms in an isotropic harmonic trap the ma- 
trix element involving the inter-channel coupling W can be 
inferred from Eq. d40b via the general approximate relation 



K0res|W|0 l ,)| 2 ^2 m «- 



»\W\4>o 



(+)■ 



,2 8Ey 

I dv ' 



(127) 



Here E v is the energy of the bare oscillator level of Eq. (f90b 
in the limit of low excitations in the Wigner threshold law 
regime. Under the assumption that the single particle trap 
length very much exceeds the background scattering length, 
i.e. Iflbgl/flho ^ 1, the density of states is given by dE v /dv = 
2ftu>h . Consequently, Eqs. (l40l and (I1271 > determine the 
Landa u-Zener parameter a ssociated with a trapped atom pair 
to be (lJulienne et at .1 120041) : 



cho 

5 LZ = 



Vo7z 



'ho 



a h „AB 



B 



(128) 



Figure [19] illustrates the validity of the Landau-Zener ap- 
proach in applications to Feshbach molecular association and 
dissociatio n in an optical lattice wher e the sites are filled with 
two atoms (Thalha mmer et a/ll2006l) . 

On the basis of Eq. dl 1 3b , the preceding determination 
of the asymptotic behaviour of transition amplitudes may 
be extended to an arbitrary numb er of energy states of the 
entrance channel inte ra ction (Demkov and Osherov, 1 1968 ; 
Macek and Cavagneroi 19981: Yurovskv a nd Ben-R euven , 



19981: lYurovskv et all Il999al) . This shows that despite the 
fact that the Landau-Zener approach neglects the bare excited 
and molecular levels, it gives the exact probability for the 
loss of atom pairs from the v = mode of a harmonic trap. 
At the end of an asymptotic linear downward sweep of E res 
all population either remains in the initial v = state or is 
transferred either into energetically lower levels, v < 0, or 
into the resonance state. A similar statement applies in the 
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FIG. 19 (Colour in online edition) Efficiency of Feshbach molec- 
ular association (upper panel) and dissociation into the lowest en- 
ergy band of an optical lattice (lower panel) via linear magnetic field 
sweeps acro ss the 1007 G zero energy resonance of 87 Rb versus the 
ramp speed ( Thalhammer et a/.Ll2006T) . The solid curves refer to the 
Landau-Zener formula of Eq. d!24b and the coefficient for a harmonic 
oscillator of Eq. (128) using the measured frequency of v ho = 39 kHz. 
The dashed curve indicates the same Landau-Zener prediction scaled 
by a factor of 0.95 which accounts for possible deficits in the ideal 
diatomic filling, for instan ce, due to tunnelling between lattice sites 
( iThalhammer et allfxtixb . 



case of asymptotic upward sweeps which are also exactly 
described by Eqs. ( 1124b . J 125b and ( 1 1 28b . Consequently, 
in the limits t\ — > -oo and ff — > oo transitions between 
states occur only in the intuitive direction of the sweep. The 
intermediate dynamics, however, involves all levels, i.e., the 
amplitudes associated with unintuitive transitions interfere 
away only at asymptotically large times. 



B. Magnetic field sweeps in Bose-Einstein condensates 

Several pioneering studies on the p roperties of zero en- 
ergy resonances of alkali atom pairs dCornish et all l2000t 



i pa 

llnouve etall 1 19981: IStenger et a/ill 9991) have been per f ormed 
in dilute Bose-Einstein co n densates ( Anderson et all 1 19951: 
iBradlev et all Il997l Il995t iDavis et aZ.l Il995l) . Condensa- 
tion of Bose atoms dBosel 1 1 924t lEinsteinl 1 1 924 fl925 ) occurs 
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when the occupation number N c of a particular single parti- 
cle state becomes comparable to the total number of atoms, 
N, in such a way that N c /N remains finite in the thermody- 
namic limit. The signatures of this phenomenon in dilute al- 
kali gases have been identified, for instance, through a spe- 
cific, narrow momentum distribution of the atom s or a char- 
acteri stic spectrum of collective excitations (Dalfov o et all 
1999). While low collision momenta facilitate the theoret- 
ical description of the two-body physics of molecular asso- 
ciation, the complex many-particle nature of Bose-Einstein 
condensates subject to dynamically resonance enhanced inter- 
actions becomes particularly significant on long time scales. 
In addition, the Feshbach molecules of some species of al- 
kalis associated with closed-channel dominated resonances 
proved to be unstable in the environment of a gas dDiirr et all 
l2004at iHerbig et all 120031: IXu et all 120031) . For this rea- 
son, such experiments were performed, in part, under condi- 
tions of free expans ion dDurr et all 2004al iMark et all 2005 ; 
lYurovskv and Ben-Reuvenu 2004, 12005b . Comparisons be- 
tween predicted and measured molecule production via mag- 
netic field sweeps in Bose-Einstein condensates are, therefore, 
less conclusive than in the case of a single atom pair confined 
to an optical lattice site. 



depletion dGoral et aZ.ll2004l) : 



4nh 

N l0S s=27r(N-l)(Nm 

m 



a\, g AB 



B 



(130) 



Here Nj'V is the uniform density of the gas in the volume 
e V. In accordance with the local density approximation, the 
total number of atoms lost from a Bose-Einstein condensate is 
given by the spatial average over the densities of all boxes. 

We note that the Landau-Zener estimate of Eq. ( 11301 ) is 
applicable to both sweep directions of the Feshbach reso- 
nance level. In the case of downward sweeps the condensed 
atoms are partly transferred into diatomic Feshbach molecules 
whose final number, M, equals one half of the atom loss of 
Eq. JT3Qb . i.e., 



Nl = Moss/2. 



(131) 



Upward sweeps lead to the production of correlated pairs with 
a comparatively high relative velocity depending on the ramp 
speed. In some experime nts using nonlinear field variations 
(iDonley et al. l l200lL l2002h such atom pairs were detected as 
a trapped dilute cloud with an average spatial extent much 
larger than the size of the remnant Bose-Einstein condensate. 



1 . Limit of high ramp speeds 



2. Two-level mean field approach 



In the limit of high ramp speeds the transfer of con- 
densed atoms into Feshbach molecules may be estima t ed us- 
ing the two-body Landau-Zener approach dGoral et all l2004t 
iMies et all |2000|) . To this end, it is instructive to divide the 
dilute gas into regions of virtually constant density. Their vol- 
umes, *V, can be chosen sufficiently large for the thermody- 
namic limit to be applicable. The state of an arbitrary atom 
pair of a uniform Bose-Einstein condensate in each of these 
periodic boxes is well described by the lowest quasi contin- 
uum energy level. Accordingly, the associated Landau-Zener 
parameter can be inferred from Eq. (1126b using the following 
replacement of the initial wave function: 



\<h) -» l0o +) ) V(2^) 3 /^. 



(129) 



Given that the size of a dilute gas is on the order of several |/m, 
the volume *V is sufficiently large for the transition probability 
of Eq. ( 1124b to reduce just to its first order approximation, 
Po.res ~ 2^lz- A typical order of magnitude of /?o,res is 10~ 6 



dMies et artkoOOh . 

While the association of a particular pair of condensed 
atoms is a rare event, the fact that each atom has all the others 
to interact with grossly enhances the efficiency of molecule 
formation. Just in the limit of high ramp speeds the Bose- 
Einstein condensate may be treated as a reservoir whose total 
atom number is hardly affected by an asymptotic magnetic 
field sweep across Bq. This assumption implies that the small 
fraction of lost atoms is well described in terms of the pairwise 
average of the microscopic transition probabilities. The num- 
ber of pairs in a box with N atoms is N(N - l)/2 which in the 
limit \B\ — > 00 yields the following estimate for the condensate 



A description of atom loss from a Bose-Einstein condensate 
consistent with both Eq. ( 1130b and the dynamical depletion of 
pairs during a magnetic field sweep across Z?o may be based 
on Eqs. (|94| l and d95l l. Such an extended two-level configura- 
tion interaction approach to the probability amplitudes, Co(f) 
and C res (t), associated with the zero energy mode of an atom 
pair and its depletion, respectively, is g iven by the formulae 
dGoral et q/.L[2004llJulienne et all 120041) : 



ihC (t) = EaCo(f) + (Nm l/2 g; es q } (t)C ies (t), 
ihC ms (t) = £ res (f)C res (f) + (Nm l/2 g les C 2 (t). 



(132) 
(133) 



Here N is the total number of atoms of the homogeneous gas 
in the volume r V. The inter-channel coupling is determined 
by the matrix element 



g KS = (2;r/O 3/2 <0 res |W|C>, 



(+)\ 



(134) 



in accordance with Eqs. d94l ), d95l l and ( 1129b . Similarly to the 
two-body configuration interaction approach, Eqs. (11321 ) and 
(11331 ) lead to a constant of motion, |C (f)| 2 + |C res (f)l 2 = 1- 
This implies an interpretation of the quantities Af|Co(f)| 2 and 
N\C res (t)\ 2 in terms of the numbers of atoms associated with 
the remnant Bose-Einstein condensate and correlated pairs, 
respectively. The nonlinear nature of Eqs. ( 11321 ) and ( 11331 ) 
in terms of the Bose enhancement factor, (N/ r V) 1 ' 2 Co(f), en- 
sures consistency with Eq. ( 1130b and, therefore, accounts for 
the surrounding gas. 

The long time asymptotic populations may be estimated an- 
alytically based on a linearised version of Eqs. d 1 32b and ( 1133 
using a static Bose enhancement factor, (NfV)^ 2 (Mie s et al. 
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2000). This treatment of the inter-channel coupling neglects 
the time dependence of the depletion of pairs of condensed 
atoms and leads to dynamical equations formally equivalent 
to those of the two-body Landau-Zener approach. In accor- 
dance with Eq. ( 1130b . the associated Landau-Zener coefficient 
is given by 



cBEC 
5 LZ 



(N/<V) ■ 



Anh 



a\,„AB 



B 



N6 



LZ- 



(135) 



The asymptotic condensate depletion and its remnant uniform 
density can be inferred from Eqs. (11241 i and ( 11251 ). respec- 
tively, using (5^ c instead of £ LZ . 

An approach similar to Eqs. (11321 i and ( 1133b but applicable 
to trapped gases beyond a local density treat ment was derived 



on the basis of a mean field approx imation (Drummond et al. 



1998 



Timmermans et al. 1 119981 Il999al 1 Tommasini et al. 
19981) . The associated many-body model Hamiltonian 
was originally introduced in the context of superconductiv- 
ity ( jFriedberg and Leel 1 19891; [Ranninger and Robaszk iewiczi 



tions: 

M¥(x, t) =// G p¥(x, t) + g* KS ¥*(x, O^esCx, t), (136) 
a*„(R, t) =H m (ff¥ m QL, t) + g res Y 2 (R, t). (137) 

Here the background scattering is included in terms of 
the usual Gross-Pitaevskii mean field Hamilto nian, Hpp, in 
the contact pseu do interaction approximation dGrossl 1 196 it 
jPitaevskiH Il96ll) . Given a spherically symmetric harmonic 
atom trap, Hep therefore consists of the following contribu- 
tions: 



Hgp = — 



h 2 V 2 



2m 



Ajrtfi 

+ — abgimof 

m 



(138) 



Typical frequencies Vh = a>ho/(2x) associated with such com- 
paratively weakly confining traps are on the order of 100 Hz. 
The generalised resonance energy, H KS (t), contains the centre 
of mass kinetic energy of correlated pairs as well as a mag- 
netic field shift from B res to the measurable position of the sin- 
gularity of the scattering length. Its explicit expression reads: 



h 2 V 2 
Am 



+ mu 2 \R\ 2 + fi res [B(t) - B ]. (139) 



The mean fields ^(x, f) and T^e^R, f) refer to the amplitudes 
of the densities of atoms in the condensate at the position x 
and of correlated pairs with the centre of mass R, respectively. 
Similarly to the configuration interaction approach, Eqs. ( 1136b 
and ( 1137b give rise to a constant of motion consistent with the 
conservation of the total number of atoms: 

J" dx PF(x, t)\ 2 + JdR |¥ res (R, t)\ 2 = N. (140) 

Since its first applications in the context of Feshbach res- 
onances in the physics of cold gases, the concept of this 
two-level mean field approach has been continually extended 
to a variety of physical situations (Ivan Abeelen and Verhaarl 
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11985 ). This procedure leads to the following dynamical equa- .2 
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FIG. 20 (Colour in online edition) Loss of condensate atoms 
in upward sweeps of the Feshbach resonance level across the 
853 G and 907 G zero en ergy resonances of 23 Na versus the in- 
verse ramp speed, 1/|S| 1 Steng eref all Il999h . The experimen- 
tal dat a are compared to theoretical predi ctions using the Landau- 
Zener dGoral et a/.Ll2004lMies et a/.L i2000) and two-level mean field 
( Ivan Abeelen and Verhaaiill999l) approaches. The theoretical Fesh- 
bach resonance parameters employed refer to the values of a bg and 
AB given in Tables Hvl and IVl We note the differences in the ramp 
speeds of three orders of magnitude between the upper and lower 
panels which reflect the different widths of the 853 G and 907 G zero 
energy resonances. 



19991 iDuine and Stood l2003bt iGoralefaZi 



20011: 



Holland et all 1200 ill Yurovskv and Ben-Reuvenll2003blr 

Figure [20] shows that Eqs. (1136b and ( 1137b give reason- 
able agreement with the loss of condensed atoms observed 
in experiments invo lving zero energy resonances of 23 Na 
(Sten ger et all [l999). These measurements refer to asymp- 
totic upward sweeps of the Feshbach resonance level across 
Bo leading to the production of unbound correlated pairs 
of atoms (Ivan Abeelen and Verhaan Il999l) . We note that 
Eqs. (1136b and ( 1137b well account for the onset of atom 
loss over ranges of ramp speeds which differ by three or- 
ders of magnitude between the upper and lower panels of 
Fig. [20] A similarly fair agreement has b een reported for pre- 
dictions of several different approaches dKohler et all l2004t 
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iMackie et all l2002al lYurovskv and Ben-Reuverl l2003al) on 
the observed atom loss of a Bose-Einstein condensate ex- 
posed to magnetic field sweeps across the 1 55 G zero en- 
ergy resonance of 85 Rb (Cornis h et at ,, 2000). The dashed 
curves in Fig. [20] indicate the asymptotic Landau-Zener pre- 
dictions based on the coefficient of Eq. ( 1135b and the local 
density approximation. Both theoretical approaches recover 
the asymptotic behaviour of the loss of condensed atoms in 
the case of fast sweeps given by Eq. (1130b . Their functional 
forms differ i n the opposite, adiabatic limit o f low ramp speeds 
(IGoral et al.ll2004t Ishkh anvan et aUl2004l) . i.e., when the pa- 
rameter 1 l\B\ in Fig. [20] increases. We note that both the two- 
level mean field approach and its associated Landau-Zener es- 
timate of asymptotic populations do not distinguish between 
the sweep directions. Consequently, the molecule produc- 
tion in a downward sweep of the Feshbach resonance level 
is treated completely symmetrically to the heating of the gas 
in an upward sweep. 

While Eqs. (1132b and ( 11331 ) refer to a genuinely two-level 
system, the two-level mean field approach accounts for the 
background scattering continuum, in principle, through the 
parameter «bg of the mean field Hamiltonian of Eq. (11381 ). 
Such a contact pseudo interaction treatment presupposes 
a separation between the typical time scales associated 
with the evolution of the Bose -Einstein condensa t e and 
the diatomic collision al duration dKohler and Burnetii 120021: 
IProukakis et a/.[[l998l) . Accordingly, the two-level mean field 
approach can be derived formally in terms of the Markov 
limit of mi croscopic many-body theories of dilute gases 
(IGoral et al 1 12004) outlined in Subsection IV.CI The as- 
sumption of a separation of time scales between two- and 
many-body evolutions is violated during a magnetic field 
sweep across a singularity of the scattering length. This 
implies that Eqs. J 1 36b and ( 1 1 37b can describe the dynam- 
ics, at most, in the asymptotic regime where the dilute 
gas parameter, [N\a{B)\ 3 /^V] 1 ^ 2 , is small compared to unity. 
Similarly to the limitations of the two-body Landau-Zener 
model, the intermediate evolution of a Bose-Einstein con- 
densate is influenced by phenomena beyo nd the range of va- 
lidity of the mean field approximation (IGoral et all l2004t 
iHolland et all 1200 ll: iKohler et all |2004 . While most of the 
theoretical approaches agree in their predictions on the con- 
densate loss in the fast sweep limit of Eq. (1130b , the satura- 
tion of molecule production is a matter o f ongoing research 
dNaidon and Masnou-Seeuws[|2003[|2006l) . 



C. Molecule production in cold Bose and Fermi gases 

Cold Bose and two spin component Fermi gases are subject 
to a considerable momentum spread which generally tends to 
reduce the efficiency of molecule production via linear mag- 
netic field sweeps across a zero energy resonance. Its sig- 
nificance is particularly obvious in the case of slow asymp- 
totic sweeps, given the spectrum of dressed energies of an 
atom pair illustrated in Fig. [18] Only the v = mode adia- 
batically correlates with the Feshbach molecular level when 
the resonance energy is decreased, while all the excited states 



undergo cooling transitions. From this too simplistic view- 
point, for instance, two spin components of a dilute vapour 
of Fermi atoms distributed according to the Pauli exclusion 
principle would produce just a single molecule in the adia- 
batic limit of the ramp speed. Contrary to the case of Bose- 
Einstein condensates, binary physics alone is therefore not 
even sufficient to qualitatively explain t he observed s ubstan- 
tial molecule produc tion in such gases dRegal et all l2003at 
IStrecker et all 120031) . 

1 . Transitions from continuum to bound states 

A quantitative analysis of the problems associated with the 
theoretical description of molecule production in the presence 
of momentum spread may be based on the exact treatment of 
linear magnetic field sweeps of Subsection llV.AI Accordingly, 
the probability for transitions from an initial dressed contin- 
uum level of a pair of distinguishable atoms in a periodic box 
of volume *V to the Feshbach molecular state at the final mag- 
netic field strength is given by 

Pass(k)= ^^|<4|f/ 2 B(ff,fi)l^ k >| ■ ( 141 > 

Here the arguments of the two-body time evolution operator, f; 
and tf , refer to the initial and final times of the sweep, respec- 
tively, and Hk denotes the initial relative momentum of the 
atoms. Similarly to the derivation of Eq. ( 1123b . the asymp- 
totic transition probability including the background scatter- 
ing continuum can be determined analytically using the sta- 
tionary phase approach in Eq. ( 11411 ) in the limits t\ — > -<x> and 
ff — > oo. This yields: 

p ass (k) = |<0 res |W|<>|V 2 «> (142) 

Herezi = h 2 k 2 /m+iQ denotes the regularised energy argument 
of the initial phase. In accordance with Eqs. ( II 10b . (II 181 ) and 
(|53b . the exponent of Eq. ( 1142b is determined by the formula 

nH 2 r 2 

Im^fe,*i) = ^— *'0(fc-£O|<0 res W<£>| . (143) 

l-^resl *J 

Here the step function of the wave numbers, 9(k-k'), indicates 
that transitions occur just in the intui tive, downward direc- 
tion o f an asymptotic sweep across Bo dDemkov and Osherovl 
1968). According to Eq. d 143b . the exponential damping of 
Eq. (1142b increases in the limit \E res \ — > 0. This confirms the 
intuitive picture suggested by Fig.[18]that adiabatic sweeps in 
a diatomic system with a continuum of modes eventually lead 
to negligible molecule production. While the momentum de- 
pendence of the exponent of Eq. ( 11431 ) recovers experimental 
dissociation spectra (Mukaiva ma et all 120031) . in the context 
of molecular association in cold gases, Eq. (11421 ) gives rise to 
exact predictions just in the fast sweep limit. 

2. Fast sweep limit of molecule production 

The onset of molecule production in asymptotic magnetic 
field sweeps across zero energy resonances is sensitive to the 



36 



statistics associated with identical atoms. Subsection lV.BI pro- 
vides a strict approach to the determination of dimer popula- 
tions on the basis of the two-particle correlation function of 
the gas. Similarly to the prediction of the loss of condensed 
atoms of Eq. ( 1130b . however, low depletions can be inferred 
intuitively from the fast sweep limit of Eq. (1142b . treating the 
gas as a reservoir of atom pairs. In this limit, \B\ — > oo, the 
damping term of Eq. ( 1143b describing cooling transitions into 
continuum levels below the initial energy h 2 k 2 /m vanishes. 
In the context of cold collisions in the Wigner threshold law 
regime the matrix element involving the inter-channel cou- 
pling in Eq. (1142b may be evaluated at k = 0. Consequently, 
the transition probability of Eq. ( 1142b becomes momentum 
independent and recovers the result of the Landau-Zener ap- 
proach, /7 ass = 2ttSlZ- 

For the purpose of studies involving s-wave collisions, 
Fermi gases are usually prepared as incoherent mixtures of 
two different Zeeman states with occupation numbers N\ and 
N2- Accordingly, N = N\ + is the total number of atoms. 
Each one of the N\ atoms of the first component has N2 atoms 
of the second component to interact with via s-wave colli- 
sions. Classical probability theory and the incoherent nature 
of the initial state therefore lead to the following estimate for 
the number of di atomic Feshbach molecule s produced in the 
fast sweep limit (IChwedenczuk et at .1 12004 : 

N\ = 2jtNiN 2 6 lz . (144) 

We note that the two-body Landau-Zener coefficient, given 
explicitly by the right hand side of Eq. ( 1135b , is inversely pro- 
portional to the volume *V of the periodic box. The fraction 
of molecules, N^/N, is therefore proportional to the density, 
N yV, which allows for an extension of Eq. (1144b to inhomo- 
geneous gases via the local density approximation. 

Given a balanced mixture of Zeeman states, i.e. N\ = N2 — 
N/2, Eqs. (1144b and (1131b show that the onset of molecule 
production in a Fermi gas is half as large as in a Bose-Einstein 
condensate of identical density and resonance parameters. In 
the limit of zero temperature such a dilute, initially weakly in- 
teracting vapour of atoms may be approximately described by 
a pair of filled Fermi seas. As fast magnetic field sweeps im- 
ply small atomic depletions, it may be argued intuitively that 
two-body cooling transitions leading into occupied modes be- 
low the Fermi energy are suppressed by the Pauli exclusion 
principle. This suggests a description of the two spin com- 
ponent Fermi sea in terms of a single level over a significant 
range of ramp speeds . The associated Landau-Zener approach 
(IChwedenczuk et aZl 120041) consistent with the linear limit of 
Eq. ( 1144b is illustrated in Fig.|2T]in comparison to experiments 
on asymptotic magneti c field sweeps across t he 543 G zero en- 
ergy resonance of 6 Li dStrecker et al.l 12003). Accordingly, a 
suppression of cooling transitions gives an intuitive explana- 
tion for the onset of a substantial molecule production in two 
spin component Fermi gases. 

In thermal Bose gases each one of the N constituents can in- 
teract with all the other N—l atoms to form a dimer Feshbach 
molecule. As opposed to the zero energy mode of condensed 
pairs, diatomic wave functions in the presence of momentum 
spread need to be explicitly symmetrised which enhances the 
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FIG. 21 (Colour in online edition) The fraction of remnant atoms, 
1 - 2N { JN, at the end of asymptotic magnetic field sweeps across the 
closed-channel dominated 5 43 G zero energy res onance of 6 Li ver- 
sus the inverse ramp speed ( Stre cker et al.l I2003T) . The solid curve 
refers to a Landau-Zener es timate for the onset of mol ecule produc- 
tion in the fast sweep limit IChwedenczuk et ah, 2004) using the lo- 
cal density approximation. The inset shows the molecular conver- 
sion, 2NI/N, observed in a thermal Bose gas with a mean density 
of 1.3 x 10" atoms/cm 3 and a temperature of 40.6 nK using mag- 
netic field sweeps a cross the 155 G zero energy resonance of 85 Rb 
dHodbvefq/.[|2005h . For comparison, the dashed curve indicates the 
pairwise ensemble average over the transition probabilities, 2p ass (k), 
with respect to the Maxwell distribution of relative velocities. This 
semi-classical estimate for the onset of molecule production in Bose 
gases is based on Eq. ( 1142b and therefore consistent with the limit of 
Eq. ( 1145b . Both the solid and dashed curves are associated with the 
parameters a bg and AB of Tables|lV]and|V] respectively. 

association probabili ty of Eq. (|141b in the thermal average by 
a factor of two dStoof et oZI|1989|) 7As the number of interact- 
ing pairs is N(N - l)/2 « N 2 /2, the fast sweep limit of the 
number of Feshbach molecules produced in a thermal Bose 
gas is given by the following formula: 

N l d = 2ttN 2 8 lz . (145) 

This estimate is illustrated in the inset of Fig. |2T]in compar- 
ison to exper i ments associated with thermal clouds of 85 Rb 
(Hodb vef all l2Q05h . We note that a similar enhancement of 
dimer formation in thermal Bose gases as compared to con- 
densates was observed in the context of inelastic three-body 
recombination and reflec ts correlation properties (Bu rt et al.l 
ll997HKaganeFaZLll985h . 

According to Eqs. ( 1144b . ( 11311 ) and ( 11451 ). the onsets of 
molecule production increase by factors of two between bal- 
anced two spin component Fermi gases, Bose-Einstein con- 
densates, and thermal Bose gases of identical densities and 
resonance parameters. The validity of these statistical es- 
timates based on two-body physics depends on the signifi- 
cance of multiple collisions of each atom during a magnetic 
field sweep. This implies that the range of ramp speeds 
described by the fast sweep limits depends on the density 
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of the gas in addition to the specific nature of the Fesh- 
bach molecules associated with entrance- and closed-channel 
dominated zero energy resonances. Predictions accessing 
wider range s of ramp speeds require genuinely many-body 
approaches dJavanainen et ~al. , 2004 iPazv ef a/Il2005l l2004t 
Willia ms et aZ.U2004allbh . ~ 



3. Saturation of molecule production 

Magnetic field sweeps sufficiently slow to convert up to 
88 % of the atoms into Feshbach molecules were employed 
to produce Bose-Einstein condensates of dimers from bal- 
anced two sp i n com ponent mixtures of 40 K Fermi atoms 
dGreiner et a/.l l2003). Figure |4] shows a typical density profile 
of such a final molecular cloud (right image) in comparison 
to a thermal distribution (left image). Motivated by the theory 
of ideal Bose gases, the conditions for condensation in these 
experiments are likely to depend on the phase space density, 
n\AL, associated with the initial dilute atomic gas of density 
n\ and temperature T. Here A± = [2ntp- 1 (nik^T)] 1 1 2 is the 
thermal de Broglie wavelength. Given the experimental den- 
sities on the order of 10 13 atoms per cm 3 and the 202 G zero 
energy resonance parameters, mag netic field sweeps wit h low 
ramp speeds of typically 160G/s (Greiner et all [2003h may 
be considered to be adiabatic. Such processes smoothly alter 
diatomic wave functions but are not necessarily expected to 
change the occupation of states in phase space. It is there- 
fore plausible that not only the condition for molecular con- 
densation but also the saturated production efficiency, 2N f d /N, 
is determined just by the phase space density rather than the 
number of atoms or their temperature individually. 

Subsequent systematic studies of adiabatic magnetic field 
sweeps in both two spin component Fermi gases of 40 K 
and thermal Bose g ases of 85 Rb have supported this view 
dHodbv et all 12005). The measured molecule production ef- 
ficiencies of Fig. [22] were analysed in terms of a stochastic 
model relying upon the assumption that the probability for 
two atoms to form a dimer depends solely on their proximity 
in phase space. Based on the experimental data, the associ- 
ated proximity conditions turned out to be virtually identical 
for both atomic species, 40 K and 85 Rb. These observations 
can be understood from first principles using coupled Boltz- 
mann equations for the Wigner functions associated with den- 
sity matrices of sep arated atoms and pair s in the resonance 
state configuration (Will iams et all 120061) . Figure [22] illus- 
trates the acc uracy of this approach including quantum sta tis- 
tical effects dBlochl 1 19281; lUehling and Uhlenbeckl Il933l) as 
well as its classical gas limit. 

Similar experiments using adiabatic magnetic field sweeps 
across the 834 G zero energy resonance in cold two spin com- 
ponent Fermi ga ses of 6 Li yielded up to 8 0% Feshbach molec- 
ular conversion (ICubizolles et all 120031) . Such observations 
were analysed in terms of theoretical approaches based on 
the assumption of therm a l equilibrium throughout t he sweep 
dChin and Grimml 120041 iKokkelmans et all\200lk . In the 
magnetic field range of positive scattering lengths about Bo 
the Feshbach molecular level gives rise to a local minimum 
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FIG. 22 (Colour in online edition) The molecule production effi- 
ciency, 2/Vj/TV, versus the peak phase space density of cold Bose as 
well as two spin component Fermi gases. The circles refer to dimer 
populations measured at the end of adiabatic magnetic field sweeps 
across the zer o energy resonance s of 85 Rb at 155 G (inset) and of 
at 202 G dHodbv et all l2005h . The curves indicate predictions 
dWilliams et all l2006h based on coupled Boltzmann equations in- 
cluding quantum statistical effects (solid curve) as well as their asso- 
ciated classical gas limit (dashed curves). For comparison, the square 
and the diamond refer to the molecular production efficiencies ob- 
served for the lowest ramp speeds in earli er sweep experimen ts using 
zero energy resonances of 4 K at 224 G 1 Re gal et all Hp03a) and of 
6 Li at 543 G dStrecker et aZll2003l) . respectively. 



of the many-particle action dSzvmanska et all 120051) which is 
associated with a meta-stable state. As this state is energeti- 
cally favourable to a gas of separated atoms, it may be argued 
intuitively that production of dimer molecules should occur 
eventually also at a stationary magnetic field strength when 
the system equilibrates. This principle has been employed 
to convert balanced incoherent two spin component mix- 
tures of 6 Li atoms int o dilute vapours of Feshbach molecules 
dJochim etal, 2003a). While energy conservation inevitably 
leads to components of comparatively hot atoms and dimers, 
the stability of universal Feshbach molecules discussed in 
Subsection MI. El allows such gases to be cooled by evapora- 
tion. The associated increase in phase space density provides 
an alternative route to the Bose-Einstein cond ensation of 6 Li2 
dJochim et all l2003bl: IZwierlein et all |2003) besides the ap- 
proach of adiabatic magnetic field sweeps. 



D. Dissociation of Feshbach molecules 

Asymptotic upward sweeps of the resonance energy lead 
to the dissociation of Feshbach molecules which often serves 
as a precursor to their detection. To this end, a molecular 
component is usually spatially separated from the environ- 
ment of a remnant atomic gas using, for instance, the Stern- 
Gerlach technique illustrated in Section UlI.GI After this sep- 
aration the dissociation allows the fragments to be detected 
conventionally using probe lasers tuned to resonance with an 
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atomic spectral line. The energy provided by the time vary- 
ing homogeneous magnetic field during the sweep is trans- 
ferred to the relative motion of the constituents of a Fesh- 
bach molecule. Such correlated atom pairs with a relative 
velocity dependi ng on the ramp speed were detected in sev- 
eral experiments fourr et adl2004btlMukaiyama et a/.[l2003t 
IVolz et all 120051) . Their spectrum of kinetic energies of the 
relative motion is given by the following expression: 



s (n 2 k 2 /m) 



mfik 



\u 2 B(t t ,tiM 



(146) 



Here dQ. denotes the angular component of afk describing the 
direction of the momentum ftk, while \<p l h ) and |^^ k ) are bound 
and dressed continuum states associated with the initial and 
final magnetic field strengths, respectively. 

Similarly to the derivation of Eqs. (1142b and ( 11431 ). the en- 
ergy spectrum of Eq. j 146b can be determined analytically in 
the asymptotic limits t\ — > -00 and ff — > 00. As dissociation 
and association are related to each other just by time rever- 
sal, their transition probability densities are identical. Typical 
energies h 2 k 2 /m of atomic fragments are on the order of //K 
in units of the Boltzmann constant which is usually inside the 
Wigner threshold law regime, i.e. fc|flb g | <s 1- This implies 
that the matrix elements involving the inter-channel coupling 
in Eqs. ( 1142b and ( 1143b can be evaluated at k = and k' = 0, 
respectively. Consequently, the asymptotic dissociation spec- 
trum of Eq. ( 1146b is well approximated by the fo llowing for- 
mula dGoral et all 12004 IMukaivama et all 120031) : 



"di S s(£) = -^ ex P 
oE 



mE \a hg AB\E 
~h 2 h\B\ 
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We note that the integral of Eq. ( 1147b over all kinetic energies 
E of the relative motion gives unity. This implies that Fesh- 
bach molecules are dissociated with certainty in an asymptotic 
upward sweep of the resonance energy, provided that transi- 
tions to the highest excited entrance channel vibrational level 
are negligible. 

According to Eq. ( 1147b . the width of the dissociation spec- 
trum increases with increasing ramp speeds. The associated 
single particle kinetic energies are usually inferred from the 
velocities of the fragments which constitute a radially expand- 
ing cloud of atoms. Their average, Ediss> amounts to one half 
of the mean energy of the relative motion of all correlated 
pairs given in terms of Eq. (1146b by the expression 
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Accordingly, the Wigner threshold law approximation of 
Eq. J 147| > determines E& K to be 
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Figure [23] illustrates the accuracy of the predictions of 
Eq. ( 1149b on average dissociation energies detected at the end 
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FIG. 23 (Colour in online edition) The mean energies of the atomic 
fragments of 23 Na 2 Feshbach molecules dissociated by a linear up- 
ward sweep of the 907 G resonan ce level versus the ramp speed. Cir- 
cles indicate experimental data (Mukaivama et al while the 
solid curve refers to the predictions of Eq. i 149b using the parameters 
«b g and Ai? of Tables lIVI andlVl respectively. 



of asymptotic magnetic fiel d sweeps across the 907 G zero en- 
ergy resonance of sodium (IMukaivama et all 120031) . While 
the solid curve refers to independ ently determined parame- 
ters abg and AB dMies et all [2000), the good agreement with 
the measurements has motivated experim ents using molecu- 
lar di s sociation to charac terise resonances dBrouard and Platal 
l2005LlDurre7aI[l2004bl) . 

We note that the asymptotic energy spectrum of Eq. ( 11471 ) 
depends just on the product a\, g AB which, according to 
Eq. d39l , determines near resonant scattering lengths in the 
universal regime of magnetic field strengths. A similar state- 
ment applies to all Landau-Zener coefficients as well as the 
long time asymptotic occupations predic ted by the two-level 
mean field approach dGoral et al. , 2004). In the limits t\ — > 
-00 and ff — > 00 the number of correlated atom pairs pro- 
duced by linear magnetic field sweeps is, therefore, insensi- 
tive to the parameters /i les and Bq - B res which directly refer to 
properties of the resonance level. This suggests that such pop- 
ulations can be inferred from any one of the single channel 
approaches outlined in Subsections IIII.FI and IIII.GI which all 
give th e correct magnetic field dependence of the scattering 
length dGoral et ad 120041: iKohler et all |2004 . All deriva- 
tions of the exact diatomic dynamics associated with linear 
magnetic field sweeps of Subsection |IV.A| may be performed, 
for instance, using a single channel Hamiltonian with a time 
dependent separable potential. Despite the fact that such a 
two-body interaction is not suited to describe closed-channel 
dominated Feshbach molecular states, it exactly recovers the 
asymptotic dissociation spectrum of Eq. ( 1147b . 

In accordance with Eq. (1122b . the range of validity of the 
limits t\ — > -00 and ff — » 00 in applications to linear sweeps is 
determined by the accuracy of the stationary phase approach. 
Given a fixed ramp speed, this approximation is violated when 
the sweep starts or terminates too close to Bo for the variation 
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of the phases in Eq. ( 11 13b to produce sufficiently many os- 
cillations of the complex exponential functions. We note that 
according to Eq. ( II 19b the phase gradient with respect to the 
energy parameter is proportional to 1/B. This implies that 
slow sweeps tend to require smaller ranges of magnetic field 
strengths about Bo for the asymptotic limits to be applicable. 
Conversely, the faster a sweep the more it resolves details of 
both the intermediate dynamics and the initial and final states. 

In the idealised limit of a jump of the magnetic field 
strength across Bo the Feshbach molecular dissociation spec- 
trum of Eq. ( 1146b just probes the overlap of the dressed 
initial bound and final continuum states. Such a scenario 
has been realised using a narrow closed-channel dom inated 
zero energy resonance of 87 Rb dDiirr et all l2004bl) . In 
these experiments the initial magnetic field strength was cho- 
sen in such a way that the Feshbach molecule was virtu- 
ally identical to the resonance state, |0 res , cl). According to 
Eq. (T3"3"l l. this implies that the dissociation spectrum is de- 
termined simply by the modulus squared of the amplitude of 
Eq. (O, i.e. n diss (p 2 /m)2/(m/?) = An\A(B u p 2 /m)\ 2 . Its reso- 
nance denominator gives rise to a sharp max imum at E KS (Bf) 
dDiirr et all l2004ri lHaque and Stool 120051) slightly shifted 
by the real part of Eq. ( l54l . while the associated imaginary 
part yields the spectral width. Such dissociation jumps across 
Bo using a resonance state of cZ-wave symmetry can populate 
several outgoing partial waves of the atomic fragments. This 
approach was employed, for i nstance, fo r the spectroscopy 
of a of-wave shape resonance dDiirr et all 120051 : IVolz et all 
120051) . 



1 . Magnetic field pulse sequence 

A typical experimental pulse sequence dDonlev et all 
2002) is illustrated in Fig. [24] All the associated magnetic 
field strengths, B, are on the high field side of the 155 G 
zero energy resonance of 85 Rb where the scattering length 
is positive. Consequently, the inter-atomic potential supports 
a loosely bound dimer state whose energy as a function of 
B is shown in Fig. Q~2] The conversion of pairs of Bose- 
Einstein condensed atoms to Feshbach molecules in these ex- 
periments crucially relies upon intermediate order of magni- 
tude variations of the scattering length, a(B), and, accord- 
ingly, of the dimer bond length, (r) m a(B)/2. The pulse 
sequence of Fig. [24] starts at the evaporation field strength 
of about 162 G at which a dilute 85 Rb Bose-E instein conden- 
sate is usually produced dCornish et all l2000h with a scatter- 
ing length of about 200 flBohr- Given the experimental mean 
inter-atomic distance of about 12000 flBohr, the gas is there- 
fore initially weakly interacting. Each pulse tunes the scat- 
tering length to fl(B m i n ) m 9000 flBohr on time scales too short 
for the Bose-Einstein condensate to adjust its density. Conse- 
quently, at the magnetic field strength fi m j n closest to the zero 
energy resonance the gas is driven into the regime of strong 
inter-a tomic interactions. During the course of the experi- 
ments dClaussen et aZ.Ll2003l:lDonlev et aZ.ll2002l) . the station- 
ary field, Bevoive, associated with the evolution period separat- 
ing the pulses as well as the duration f evo i ve were varied. The 
variation of B e voive between 156 G and 162 G corresponds to 
the range of Feshbach molecular binding energies displayed in 
the inset of Fig.[T2]with associated scattering lengths between 

about 4000 flBohr and 200 flBohr- 



V. ATOM-MOLECULE COHERENCE 

Besides the production of Feshbach molecules via 
asymptotic linear sweeps across Bq, several app r oaches 
to date rely upon ramp seque nces dMark et all 120051; 
lYurovskv and Ben-Reuvenl 120051) or r esonant oscillating 
magnetic fields dThompson et all l2005at) . Such experimen- 
tal techniques are designed to improve the efficiency of dimer 
formation in Bose gases mainly by avoiding the region of large 
scattering lengths and accordingly strong inter-atomic interac- 
tions. The associated nonlinear magnetic field variations lead 
to dependences of the molecular population on properties of 
zero energy resonances beyond those included in the Landau- 
Zener parameter. 



A. Ramsey interferometry with atoms and molecules 

Pulses starting and ending on th e positive scattering le ngth 
side of a zero energy resonance dClaussen et all 120021) are 
of particular interest in this context. Their repeated appli- 
cation to a 85 Rb Bose-Einstein condensate gives rise to the 
coherent oscillations between the final comp onents of sepa- 
rated atoms and Feshbach molecules of Fig.[2](Clausse n et all 
l2003HDonlev et aZ.ll2002l) . 
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FIG. 24 (Colour in online edition) Scheme of a typical magnetic 
field pulse sequence (solid lines) employed in experiments produc- 
ing coherent oscillations between final ato mic and molecular c om- 
ponents in a 85 Rb Bose-Einstein condensate iDonlev et fl/.Ll2002l) . A 
pair of pulses with a minimum field strength, B m j„ = 155.5 G, is sep- 
arated by an evolution period of variable duration, f cvo | vc , and field 
strength, S cvo | vc . The zero energy resonance position, B ss 155 G, is 
indicated by the dashed line. 

All measurements probing the densities of atoms were per- 
formed after each pulse sequence had terminated and there- 
fore reflect the state of the weakly interacting gas at the final 
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magnetic field strength of about 162 G in Fig. [24] According 
to these observations, the fast magnetic field variation gives 
rise to the three components of the atomic cloud illustrated in 
Fig. [2] The occupation numbers associated with the remnant 
Bose-Einstein condensate, burst component and undetected 
atoms all oscillate with respect to each other as a function 
of the evolution time, f evo ive- It turned out that their common 
angular frequency, a> e , is accurately determined by the Fesh- 
bach molecular energy in the evolution period, EF voIve , via the 
relation a> e = \E™ lwe \/h. Donley et al. concluded that the 
undetected atoms were transferred into Feshbach molecules 
whose fast phase evolution as compared to the atomic compo- 
nents leads to interference fringes in the final occupations. 

In accordance with such an intuitive explanation, the first 
pulse provides the overlap between the dimer size and the av- 
erage distance between the condensed atoms which is crucial 
to the molecular association. As this fast magnetic field varia- 
tion is not resonant with the binding energy, it inevitably leads 
to the additional production of excited atom pairs which con- 
stitute a burst of atoms. During the evolution period all com- 
ponents of the weakly interacting gas evolve independently 
and coherently and therefore accumulate a phase difference 
of A<p = w e fevoive- Finally, the second pulse gets the Fesh- 
bach molecules and separated atoms to overlap once again 
and thereby probes Aip via an interference in their occupa- 
tions. This intuitive scenario is analogous to the principle of 
a Ramsey interfero meter using pairs of separated ato ms and 
Feshbach molecules dDonlev et q/.Ll2002l : IZolleil l2002) whose 
scheme is illustrated in Fig. [25] 



2. Dynamics of a single atom pair 
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FIG. 25 (Colour in online edition) Schematic illustration of the 
Ramsey interferometry with 85 Rb atoms and Feshbach molecules 
associated with the magnetic fi eld pulse sequence of Fig. l24l 
toonlev et flZl 120021 : IZoller. 2002). The first pulse drives the initial 
Bose-Einstein condensate into a coherent superposition of bound and 
unbound atom pairs. Between the pulses the individual orthogonal 
components evolve independently. The second pulse gets the atomic 
and molecular states to overlap and thereby probes their phase differ- 
ence. The inset shows a typical interference fringe pattern at the end 
of the pulse sequence as a function of the evolution time, t emlve . Cir- 
cles refer to measurements of the remnant Bose-Einstein condensate 



Claussen et ah, 2003), while squares indicate theoretical predictions 
Goral et all 2005). 



In analogy to the analysis of Subsection llV.CI the observa- 
tions of Fig. [2] can be qualitatively unders tood in terms of the 
dynamics of just a single pair of atoms dBorca et all 120031: 
iGoral et all 120051) in a periodic box of volume *V. To this 
end, the two-body time evolution operator of Eq. (f98T > associ- 
ated with the magnetic field variation of Fig. [24] may be split 
into its contributions of the first and second pulse and of the 
evolution period. This yields: 

U 2 B(k, fi) = U 2 (tf, f 2 )t/evolve(fevolve)£/l(fl, h)- (150) 

Here the first pulse starts at the initial time t\ and terminates at 
1 1 , while t2-t\+ fevoive indicates the beginning of the second 
pulse which ends at the final time of the sequence, ff . Accord- 
ingly, U\(ti, ti) and U 2 {tf, t 2 ) describe the two-body dynamics 
during the first and second pulse, respectively. In between the 
pulses the magnetic field strength is stationary, in accordance 
with Fig. [24] This implies that the associated time evolution 
operator, J/ e voive (/evolve), depends just on the duration ? eV oive- Its 
spectral decomposition in terms of dressed states of the rela- 
tive motion of an atom pair exposed to the magnetic field of 
strength B e voive is given by the following formula: 



ev.^evolve) = £ \ < f^ e ^ t ^{ 
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(151) 



Here the index v labels the vibrational quantum numbers of 
the box states in such a way that \(f> ev ° lve ) correlates adiabat- 
ically, in the limit of infinite volume r V, with the Feshbach 
molecular dressed state, |<^ volve ). Accordingly, the associated 
vibrational energy levels are denoted by £* volve which implies 
the asymptotic behaviour £'<-™ lve ~ ^evolve w j tn mcreasm g 
r V. The more tightly bound states associated with vibrational 
quantum numbers v < - 1 contribute negligibly to the fringe 
pattern of Fig.[25]and are therefore omitted in Eq. ( 11511 ). 

Similarly to the treatment of Subsection llV.CI the interfer- 
ence fringes in the population of the remnant condensate com- 
ponent shown in the inset of Fig.[25]can be described in terms 
of the following transition probability: 



Po,o 



(152) 



To this end, it is instructive to split the associated amplitude, in 
accordance with Eq. d 1 5 lb , into its contributions of the vibra- 
tional states with quantum numbers v > and of the Feshbach 
molecule. This yields: 



<4> f o\U2B(.t t ,t-M ) = D 2B +A 2B . 



(153) 



Due to the quasi continuum of levels, the first term, D 2 b, on 
the right hand side of Eq. ( 1153b tends to decay via phase dif- 
fusion as a function of ? eV oive- Its explicit expression reads: 



^2B=2^0 |t/2(ff ' f2)l ^ VOlVe)e 



/£™ olve fevoive /fi 



v=0 



(154) 
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The magnitude of the second term, A2B, on the right hand 
side of Eq. ( 1153b depends on the product of the amplitudes 
for Feshbach molecular association during the first pulse and 
dissociation into the zero mode due to the second pulse. Its 
evolution as a function of f CT oive stems just from the phase shift 
associated with the energy E^ ve , in accordance with the fol- 
lowing formula: 



Am =<4l^(ff,f2)ICi oh >^ 

x<0l? lve |f/i(*i,;i)l4>- 



^evolve 



(155) 



The modulus squared of the sum of interfering amplitudes, 
£>2b and A2B, therefore leads to the anticipated fringe pattern 
in its usual form, 

Po,o = |£>2b| 2 + |A 2B | 2 + 2|Z) 2B ||A 2B | sin( 

^evolve + <P), (156) 



in agreement with the illustration of Fig. [25] Here the angular 
frequency, <w e = |£ ,ev 1 olve |//i, is determined by the Feshbach 
molecular binding energy. The absolute phase, <p, is associated 
with both amplitudes Are, and £>2b and depends sensitively on 
the exact shape of each pulse. A similar statement applies to 
the efficiency of dimer production and dissociation which are 
both crucial to the fringe visibility, according to Eqs. ( 1155b 
and([B6}. 

The populations of Feshbach molecules as well as the burst 
spectrum composed of the excited dressed levels at the final 
magnetic field strength may be inferred from the associated 
transition probabilities, similarly to Eq. (11561 >. Their pair- 
wise averages in a volume *V in addition to a subsequent local 
density approximation give a good account of th e magnitudes 
of all components of the gas shown in Fig. [2] dGoral et al.l 
120051) . The range of validity of such a two-body estimate is 
set by the requirement that B em \ ve is sufficiently far from Bo 
that the mean distance between the atoms is much larger than 
the scattering length during the evolution period. As the 155 G 
zero energy resonance of 85 Rb is entrance-channel dominated, 
the two-body dynamics is well described in terms of both 
the effecti ve two- and single-channel approaches of Subsec- 
tion MM dGdraleLall [20051 ^^^eLM^mi). We note, 
however, that the contact pseudo interaction is insufficient. 
This is immediately apparent, for instance, from the inaccu- 
racy of the universal estimate of the binding energies in the 
inset of Fig. [T2]whose experimen tal data were de t ermin ed us- 
ing measured fringe frequencies dClaussen et al.l 120031) . The 
interference experiments are therefore sensitive to the van der 
Waals tail of the inter-atomic potential which is not accounted 
for by the Landau-Zener parameter associated with asymp- 
totic linear magnetic field sweeps. 



B. Number of dimers produced in Bose and Fermi gases 

Another significant difference in the descriptions of the 
linear sweep and magnetic field pulse techniques is asso- 
ciated with the treatment of the Feshbach molecular state. 
The results of Subsection IIV.AI suggest that the number of 
dimers produced via asymptotic magnetic field sweeps may 
be inferred from two-level approaches effectively identifying 



dressed Feshbach molecules with bare resonance states. For 
the Ramsey interfer ence experiments dClaussen et al I I2OOI 
iDonlev et al l l2002h operating just in the vicinity of the zero 
energy resonance a maximum final molecular conversion of 
about 16% was reported under the conditions referred to 
in Fig. [2] The closed channel admixture to the dressed 
bound state, however, is appro ximately 25 % at the final mag- 
netic field strength of 162 G (Ko hler et all 20041) and even 
smaller thro ughout t he pulse sequence ( Braaten et al.l 120031: 
Kohler et al. , 2003b). In order to predict the observed number 
of dimers, it is therefore necessary to include the long range 
nature of the Feshbach molecules produced in the description 
of their detection. 

An associated generic approach suitable for both magnetic 
field sweeps and pulse sequences may be based on the follow- 
ing observable for the populati on of any diatomic state, \<pd), 
in a gas consisting of N atoms dKohler et al.l l2003ah : 



1 N 



(157) 



Here the indices i and j label the atoms, \<pjj) refers to the 
two-body state, \<f>^), associated with the atoms i and j, and 
the factor of 1 /2 on the right hand side of Eq. ( 1157b prevents 
double counting of pairs. Given any state of the gas generally 
described by a density matrix, p^B , the number of atomic pairs 
in the state \(/>&) is consequently determined by the following 
expectation value: 



N d = <N d > = TrCNaPivB). 



(158) 



Here the symbol "Tr" refers to the trace over the degrees of 
freedom of all atoms constituting the gas. 

Many-body systems consisting of identical atoms are usu- 
ally described using the approach of second quantisation 
dFetter and Waleckal 197 lb . Accordingly, any N particle state 
is constructed by repeated application of creation operators to 
the vacuum, |Q), which refers to the absence of atoms. The 
associated field operators, if/ a (x) and i/^(y), annihilating and 
creating single atoms at the positions x and y, respectively, 
fulfil the following (anti-)commutation relations: 



<Mx)<Al(y) + (w¥<*W = <V<5(x - y)- 



(159) 



Here the minus sign is associated with bosons while the plus 
sign refers to fermions, and the Greek labels a and /3 indicate 
their Zeeman states. Given the diatomic wave function of the 
relative motion, (^(r) = (r; a,/3\(pd), the second quantised 
representation of the number operator of Eq. J 157b reads: 



N, 



afi.a'fi' J J J 

x^(x)^(y)^(y')^(x'). (160) 



Here the single-particle coordinates are given in terms of the 
centre of mass and relative positions by x = R + r/2, y = 
R - r/2, x' = R + r'/2, and y' = R - r'/2. 
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The approach of Eq. ( 1160b was introduced i n the context 
of sim ulations of the Ramsey fringes of Fig. [2] ( Kohler ef al., 
l2003allbh and subsequently applied to th e molecule produc- 
tion via magnetic field sweep s in Bo se dGoral et al .1 120041) 
and Fermi gases (Peral i et al. I 120051) . In accordance with 
Eq. (1160b . the expectation value of Eq. ( 1158b depends on the 
two-body correlation function, 

G 1W X '' y' ; x ' y } = <4w^(y)^(y')^(x')>. (161) 

Here the average refers to the final state of the gas, pNB(tf)- 
In the context of the dimer production via fast magnetic field 
sweeps of Subsection IIV.CI for instance, the dynamics of 
G® is determined simply by the two-body evolution opera- 
tor, U2s(tf,ti), in the limit of short interaction times. Given 
an ideal gas initial state, pNBih), associated with either Bose 
or Fermi atoms, the exact form of the initial two-body corre- 
lation function c an be inferred from Wick's theorem of statis- 
tical mechanics ( Fett er and Waleckal Il97ll iMatsubarai 1 19551: 
IWickLll950h . In the single resonance approach, the diatomic 
channel wave functions of the Feshbach molecular state oc- 
curring in Eq. (1160b are determined in terms of the entrance- 
and closed-channel components by the following formula: 



|bg)^ g (r) + |cl>^(r). 



(162) 



afi 



Based on these assumptions, the number of Feshbach 
molecules predicted by Eq. ( 1160b strictly confirms the fast 
sweep limits of Eqs. ( 1144b . ( 1131b and (|145b , in particular, 
their coefficients associated with the identical nature of the 
atoms. Besides these statistical estimates, the general observ- 
able of Eq. ( 1160b may be applied to a variety of dynamical 
magnetic field variations as well as physical quantities. Given 
the associated two-body correlation function, it recovers, for 
instance, not only the measured maximum conversion into 
bound dimers of about 16 % in Fig.[2]but also the populations 
of correlated p airs in continuum lev els which constitute the 
burst of atoms (Koh ler et ail l2003al) . It turns out that due to 
the nonlinear magnetic field variation of Fig. [24] the determi- 
nation of the time dependence of G (2) requires a description 
via techniques beyond the two-level mean field approach of 
Subsection|IVB] 



C. Dynamics of partially condensed Bose gases 

The precise many-body dynamics underlying all obser- 
vations of Fig. |2] can be derived from the time dependent 
Schrodinger equation dete rmined by the general Hamiltonian 
dFetter and Waleckal fT97lT) : 



a J 
+ Hm t . 



2m 



+ El + V^(x) 



i// a (x) 



(163) 



Here refers to an atomic Zeeman energy whose typical 
magnetic field dependence is illustrated in Fig.[l2]for the ex- 
ample of 87 Rb, and V^°(x) describes the harmonic confine- 
ment due to an atom trap. In the context of dilute gases, the 



inter-atomic interactions are predominantly pairwise which 
implies the following potential energy contribution: 

H i*t = ~ 2 f dx f dyV aM3 ,(x-y) 



xi/^( x )(/ay)iMy)<Aa'(x). 



(164) 



Here V a /3,a'j3'(x - y) denotes the two-body potential associ- 
ated with the incoming and outgoing spin channels \a',/3') 
and \a,f3), respectively. We note that Eqs. ( 1163b and ( 11641 ) 
are formulated sufficiently generally to treat the inter-atomic 
interactions on a microscopic level, similarly to the coupled 
channels theory of Subsection IIII.AI As a consequence, the 
atomic Zeeman levels, E^, appear separately in Eq. d 163b . 
while in the two-body Hamiltonian of Eq. (TTBT l the dissoci- 
ation threshold energy associated with the closed channel is 
included in the potential V C \(B, r). Accordingly, all interac- 
tions, V a p, a 'p'(x - y), vanish in the limit of infinite distances, 
r — |x - y| — > °o, similarly to the interaction matrix V m t(r) of 
Eq. (HQ. 



1 . Beyond mean field approaches 

Several techniques were employed in treatments of the 
many-body dynamics of dilute gases beyond the Gross- 
Pitaevskii and two-level mean field approaches of Subsec- 
tion IIV.BI These methods involve the Sch winger-Keldysh 
formalism (iKeldvshl fl965l ISchwingerl 1 1 96 ll) suitable for de- 
scriptions of phenomena associated with the evolution to- 
ward thermal equilibrium, such as the dimer production via 
adiaba tic magn etic field sweeps of Fig. l22l (Willia ms et all 
2004a b, 2006). Practical approaches describing equili- 
bration on even longer time scales, i.e. beyond the range 
of validity of the generalised Boltzmann equations derived 
from the Schwinger-Keldysh t heory, may be based on the 
two-particle irreducible action jBavml 119621: ICornwall et all 
Il974t iLuttinger and Ward Il960h . In the context of the dy- 
namics of cold gases, such techniques have been imple- 
mented, to date, using th e contact pseudo interaction in a 
single spatial dimension ( Gasen zer et al. , 2005I Hex et ail 
120051) . While the Ramsey fringes of Fig. [2] are sensitive to 
parameters of the inter-atomic potential besides the scattering 
length, the associated typical pulse sequence of Fig. [24] in- 
volves time scales sufficiently short for equilibration phenom- 
ena to be negligible. This implies that the dynamics of the 
gas is captured by extensions of mean field theory which ac- 
count for the two-body time evolution beyond the Markov ap- 
proximation dHolland et all 1200 it [Kohler and Burnettl 120021: 
iNaidon and Masnou-Seeuwa 12006: Pro ukakis et ali \1998). 

The usual quantities of interest in these short time quantum 
kinetic approaches may be expressed in terms of correlation 
functions, i.e. expectation values of normal ordered products 
of field operators, such as G (2) of Eq. ( 11611 ). Their general ex- 
pression reads (t/^(y) ■ ■ • if/ a (x)} t where all creation operators 
appear to the left of all annihilation operators, and the aver- 
age refers to the many-body state at time t. The number of 
field operators constituting the product shall be referred to in 
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the following as the order of the correlation function. Based 
on the Schrodinger equation, the associated dynamics is de- 
termined by 

Here // is the Hamiltonian of Eq. (1163b . and the symbol 
[A,B] = AB - BA indicates the commutator of the operators A 
and B. Due to the potential energy contribution of Eq. (1164b . 
the commutator on the right hand side of Eq. ( 1165b gives rise 
to products containing two more field operators than the ex- 
pectation value to the left. This implies that the dynamics 
of any one correlation function is determined by a coupled 
set of equations involving all the others. Approximate solu- 
tions to Eq. (1165b . therefore, often rely upon schemes for the 
truncation of the associated infinite hierarchy of dynamical 
equations. The range of applicability of such approaches is 
sensitive to the initial state of the many-body system. 

In the context of weakly interacting gases close to thermal 
equilibrium, an approximate form of all correlation functions 
may be determined using Wick's theorem of statistical me- 
chanics dFetter and Waleckal Il97ll iMatsubaral 1 19551 IWickl 
To this end, it is instructive to introduce the connected 
correlation functions (Weinberg] [Ull), sometimes referred to 
as cumulants. These quantities may be inferred recursively 
from a decomposition of each correlation function into a sum 
of all possible products of cumulants which preserve the or- 
der of appearance of the operators. Given a set of Bose field 
operators, A, B and C, for instance, the first three cumulants, 
denoted by (A) c , (BA) C and (CBA) C , are determined implicitly 
by the following relations: 

(A) =<A> C , (166) 

(BA) =<BA) C + {A) C {B) C , (167) 
(CBA) =(CBA) C + {BA) C {C) C + {CA) C {B) C 

+ (CB) C (A) C + (A} C (B) C (C) C . (168) 

The cumulants associated with Fermi field operators may be 
inferred from expressions similar to Eqs. (1166b . ( 1167b and 
(1168b . In accordance with the anti-commutation relation of 
Eq. (1159b . however, the number of commutations of field op- 
erators in the decompositions needs to be accounted for by an 
appropriate sign of each contribution. In both cases, Bose and 
Fermi atoms, Wick's theorem reduces to the statement that 
all cumulants containing more than two field operators van- 
ish provided that the gas is ideal and in thermal equilibrium. 
The magnitude of higher order cumulants therefore provides a 
measure of the deviations of the state from the grand canonical 
density matrix in the absence of inter-atomic interactions. 

This suggests that in the context of dilute, weakly inter- 
acting gases a reasonable quantum kinetic approach may be 
based on the transformation of Eq. ( 1165b into the associated 
set of dynamical equations for cumulants. In accordance with 
Wick's theorem, this coupled system allows for an approxi- 
mate truncation at any order of correlation usually dete rmined 
by ex ternal driving fields as well as the initial state dFrickel 
Il996l) . The dilute gas of 85 Rb atoms in the Ramsey interfer- 
ometry experiments illustrated in Fig. [25] was prepared as a 



Bose-Einstein condensate in the (/ = 2,mj = -2) Zeeman 
level. Such a coherent initial state gives rise to a mean field, 
*T„(x, t) = (iff a (x)),, describing the density associated with 
the macroscop i cally occupied mode via its modulus squared 
(iDalfovo et a/.l[l999l) . In accordance with Fig. [2] the conden- 
sate is depleted due to the near resonant magnetic field vari- 
ation of Fig. [24] This implies that, in addition to the mean 
field, a minimum set of cumulants describing this experiment 
is given, respectively, by the pair function and the one-body 
density matrix of the non-condensed component: 

4V(x, y, f) = <<My)<Mx)>, - 0Y/j(y, t), (169) 
T^Cx, y, t) = <^(y)^(x)>, - ¥ ff (x, r)^(y, t). (170) 

Accordingly, at any time, f, the density of atoms in the Zee- 
man state with the index a is given by IT^x, f)l 2 + r aa (x, x, f). 
Due to the weak interactions of the gas at the beginning of the 
magnetic field pulse sequence, the second order cumulants of 
Eqs. ( 1169b and ( 1170b are negligible initially, i.e. ^^(x, y, t{) = 
= r ff/3 (x, y, fj). 

To a first approximation, the dynamical equations ( 1165b 
may be transformed and truncated in such a way that they 
just include products of normal ordered cumulants contain - 
ing at most three field operators (Kohl erand Burnetii 120021) . 
This yields the following relation for the time derivative of the 
condensate mean field: 

iti%(x,t) = H™V a (x,t)+ ^ f dyVawrir) 

x^(y,«)(fr(y)^«> ( . (i7i) 

Here r = x - y refers to the relative coordinates of a pair 
of atoms at the positions x and y, and the one-body Hamil- 
tonian, H^ B , consists of the kinetic and Zeeman energies as 
well as the trap potential of Eq. ( 11631 ). The correlation func- 
tion, (iffp>(y)iff a '(x))u on me right hand side of Eq. ( 1171b is 
determined in terms of cumulants by Eq. ( 1169b . In addition 
to Eq. dl71| ), the dynamical equation associated with the pair 
function reads: 

m^ ap {x, y, t) = i H Ha'p'^'( x ' y. f ) 

a', /3' 

+ V M (r)WA^ tyP?(y, ?)]• (172) 

Here H^ a ,p, denotes the two-body Hamiltonian matrix asso- 
ciated with the incoming and outgoing spin channels \a',B') 
and \a,B), respectively, describing both the centre of mass and 
relative motions of an atom pair. 

Given this first order truncation scheme, Eqs. dl71| ) and 
dl72| ) uniquely determine the condensate mean field as well 
as the pair function. It turns out that the dynamical equa- 
tion associated with the one-body density matrix of the non- 
condensed component may be solved implicitly in terms of 
the pair function. This yields: 

r ap (x,y,t) = Y J I dzO„ r (x,z,f)[<I>e y (y,z,f)r. (173) 
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In accordance with the general expression for the associated 
observable, N = J dx ifr a (x)ifr a (x), the expectation value of 
the number of atoms, N, is strictly conserved by Eqs. (11711 ). 
( TT721 and ( fT73l at all times, i.e. 



x,oi 2 + r fl 



,(x,x,f)] 



N. 



(174) 



In addition, Eqs. ( 1 170b and ( 1173b preserve the positivity of the 
one-body density matrix (i^i(y)t/f a (x))t, a consequence of the 
unitary time evolution, which is not necessarily recovered by 
higher order approximation schemes. 

In the context of dimer formation in partially con- 
densed Bose gases, similar extensions of mean field 
theory are based on the H artree-Fock Bogoliubov ap- 
proach (IHolland et al.l 1200 lb as well as the reduced pair 
wave function approximation (IChernv and Shanenkol l2000t 
iNaidon and Masnou-Seeuwsl, [2003. 2006|). As these methods 
all originate from Eq. (1165b . their short time asymptotic limits 
agree with the exact result given by the perturbation expan- 
sion of the dynamical many-body Schrodinger equation. In 
addition, the corr elation functions predicted by any one of 
these approaches (IHolland et al.l 1200 ll IKohler a nd Burnett, 
l2002t INaidon and Masnou-Seeuwsl 120031 120061) are free 
of secular long time asymptotic behaviour, a common 
artifact of perturbation theory associated with a spuri- 
ous polynomial time de pendence. Both the cumulant 
( Kohl er and Burnet! 20021) and the reduced p a ir wav e func- 
tion approach (Naidon and Masnou-Seeuwsl 120031 120061) 
were formulated in such a way that they are compatible with 
the use of microscopic potentials beyond contact pseudo in- 
teractions. 

Most implementations of quantum kinetic approaches to 
the dimer production via Feshbach resonances to date are 
based on single-channel or two-channel single resonance bi- 
nary interactions illustrated in Subsection llll.BI This implies 
that the Bose-Einstein condensate mode is described by a sin- 
gle mean field, T^x, f), associated with the Zeeman level in 
which the gas is prepared. In addition, the entrance- and 
closed-channel components of the pair function may be in- 
ferred, similarly to Eq. ( 11621 i, from the following formula: 

J] |a,/J><Mx, y, t) = |bg>d) bg (x, y, f) + |cl>4>ci(x, y, t). (175) 



The resonance mean field of Subsection IIV.BI is determined 
by the relation O c i(x,y, t) = ^f KS (fi, i)0res(/)- Here R = 
(x + y)/2 refers to the centre of mass coordinates of a pair 
of atoms and r = |x — y| denotes their relative distance. Based 
on a two-channel implementation of Eqs. ()171b and (1172b . 
the two-level mean field approach (tPrummond et a 



Timme rmans et al ., 1998, 1999a; Tommasini et al. 



may be recovered by formally solving the dynamical equa- 
tion associated with the entrance-channel component of the 
pair function. A subsequent elimination of <&be(x, y, t) from 
the coupled equations (1171b and ( 1172b yields the functional 
form of the right hand sides of Eqs. (1136b and dl37b . The 
coupling constants of Eq. ( 1134b and of the entrance channel 



contact pseu do interaction are d etermined by the Markov ap- 
proximation (iGoral et a/.ll2004l) . 

We note that truncation schemes of higher order than 
Eqs. ( 1171b and ( 1172b can involve cross coupling terms that 
give rise to pair functions beyond the two-channel decom- 
position of Eq. ( 1175b . Such a scenario occurs in quan- 
tum kinetic approaches associated with both Bose and Fermi 
gases provided that one and the same atomic Zeeman state 
is shared between the two-body entran ce and closed channels 
(Bruu n et al .1 120051: iParish et al. I l2005h . This, in turn, restricts 
the applicability of common l y employed model Hamiltoni- 



ans (Friedb ere and Leelll989tlRanninger and Robaszk iewiczl 



11985 ) which separate out the resonance state in the description 
of Feshbach molecule production. 



2. The remnant Bose-Einstein condensate 

The dynamics of the condensate mean fi eld associated 
with the Ramsey interferometry experiments (iDonlev et al.l 
120021) w as describe d using the Hartee-Fock Bogoliubov 
method jKokkelman s and Holland! |2002|) and related tech- 
niques dMackie et a/.ll2002bl) . as well a s the first order cumu - 
lant approach of Eqs. (fTTTb and ( fT72l dKohleref all l2003al) . 
Such studies involving beyond mean field theories, can of- 
ten be simplified by eliminating the pair function from the set 
of dynamical o r associated eigenvalue equations, respectively 
(iBurnettl 1 1999b . In the context of Eq s. dTTTb and (fT72l , this 
procedure yields (IKohler and Burnera 12002): 



iW¥(x, f) ^//jb^x, t) 

Jr»oo q 
dt ,y ¥ 2 (x,t')—h(t,t'). (176) 
* of 



Here H\% is the one-body Hamiltonian associated with the ini- 
tial Zeeman state of the condensed atoms. The coupling func- 
tion on the right hand side of Eq. ( 11761 ). 

h(t,t') = 6(t - f'X27r/0 3 <0,bg|V£/ 2 B(f,f')|0,bg>, (177) 

is determined by the two-body time evolution operator of 
Eq. ( 1150b . Here (r|0) = l/{2nh) 3 ^ 2 denotes the zero momen- 
tum plane wave of the relative motion of an atom pair and V 
is the microscopic potential matrix of Eq. ( 1164b . 

The representation of Eqs. ( 1176b and dl77b shows that all 
findings of Subsection IV.AI about the two-body time evolu- 
tion and physical origin of the Ramsey fringes are included 
in the quantum kinetic approach. Its degree of accuracy in 
comparison to experimental data of the condensate compo- 
nent remaining at the end of a magnetic field pulse sequence 
dClaussen et al is illustrated in the inset of Fig. [25] The 

associated theory curve is based on an implementation of the 
two-channel singe-resonance approach of Subsection lHI.GI us- 
ing a separable background scattering potential dGoral et al.l 
2005). 
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3. Feshbach molecule and burst components 

In accordance with Eq. dl74l ). the total number of atoms can 
be decomposed into a mean field contribution as well as a non- 
condensed component, N nc (i), described by the density matrix 
of Eq. ( 1170l i. The physical signi ficance of N nc (t) in t he Ram- 
sey interferometry experiments dDonlev et al. , 2002) may be 
inferred from Eq. ( 1173b . To this end, it is instructive to replace 
the spatial average of Eq. ( 11731 ) using the completeness of the 
set of dressed bound and continuum energy states associated 
with the magnetic field strength at time t. This leads to the 
following relation: 

N us (t) = J"</r[ J dp |<R,0 P |<5(O>| 2 

+ |(R,0 b |(D(O>| 2 ]. (178) 

Here R may be interpreted in terms of the centre of mass po- 
sition of an atom pair, and |0(f)) refers to the pair function 
whose channel components in the spatial representation are 
given by Eq. J 169b . For simplicity, the spectral decomposition 
of Eq. ( 11781 ) includes just the dressed continuum states, \<p p ), 
and the Feshbach molecular state, |</> b ). More deeply bound 
levels are neglected. 

Equation ( 11781 ) suggests that the continuum contribution 
may be interpreted in terms of correlated atom pairs associated 
with relative momenta, p, while the bound state part yields the 
number of atoms converted into dimers. An analy sis based on 
the ob servable of Eq. j 160b confirms this view dKohler et ali 
l2003al) . Its practical implementation based on Eqs. dl711 l and 
( 11721 ) involves a cumulant expansion of the two-body corre- 
lation function in accordance with the first order truncation 
scheme. This yields: 



,(2) 

T a'/3',a/3 



(x',y';x,y;r)=<iA!(x)^(y)> f 

x <<My')<Mx')> ; . 



(179) 



Here the correlation functions on the right hand side can be ex- 
pressed in terms of pair functions and condensate mean fields 
via Eq. ( 11691 ). Using Eqs. ( 1162b and ( 1179b . it turns out that the 
diatomic number operator of Eq. ( 1160b gives rise to a dimer 
mean field, 



¥ b (R,t)=^r{<R,0bl*(O> 



V2 

+ jdr [^ E (r)]>(x,f)>P(y,?)}, (180) 

which determines the number of Feshbach molecules pro- 
duced via the relation N d (t) = JdR |¥ b (R, Ol 2 - 

The first term on the right hand side of Eq. ( 1180b recovers 
the bound state contribution of Eq. ( 1178b . In accordance with 
the long range of the 85 Rb2 entrance-channel wave functions 
of Fig. QT| the second term may be interpreted in terms of 
the overlap between the dimers produced and the surround- 
ing gas. Its magnitude is determined by the dilute gas param- 
eter squared, n c (f)a 3 , where « c (f) is the average condensate 
density and a refers to the scattering length associated with 



the magnetic field strength at time t dKohler et ali [2003b). 
Consequently, the existence of this overlap term reflects the 
breakdown of the concept of diatomic molecules in the envi- 
ronment of a strongly interacting gas. This breakdown occurs 
when the bond length (r) = a/2 of the universal dimer wave 
function of Eq. (O is comparable to the mean inter-atomic dis- 

— 1/3 

tance n c of the condensate. Similar phenomena are visible 
inFigures[T3[Hand[T5] 

At the end of the magnetic field pulse sequence of Fig. [24] 
the gas is weakly interacting and the second, overlap term 
on the right hand side of Eq. ( 1180b is negligible. Based on 
Eqs. ( 1176b and ( 1180b . the predicted number of atoms con- 
verted into dimers recovers the 16% maximum fraction of 
missing at oms of Fig.[2 ] as well as its modulation as a function 
of fevoive dKohler et ali |2003a). In accordance with Fig. [15] 
the fact that these Feshbach molecules were not detected re- 
flects their lifetime with respect to spin relaxation on the or- 
der of only 100 jus at the final magnetic field streng th of about 
162 G dKohler et a/.ll2005l:lThompson et adl2005bl) . An anal- 
ysis based on Eqs. ( 1160b and ( 1179b also shows that the contin- 
uum contribution on the right hand side of Eq. ( 1178b is associ- 
ated with correlated atom pairs which constitute the measured 
burst component of Fig . [2] ( K ohler et ali l2003al) . 



4. Three-component Ramsey fringes 



The 



phys i cal origin of the Ramsey fringes 
dDonlev et ali 12002) was t he su bj ect of sever a l theo- 
retical works dGoral ef ad 120051: iKohlere/ ali 12003d: 



IKokkelmans and Holland!, 120021: IMackie et ali I2002M) using 
different beyond mean field appro aches. Figure [26] shows 
such a prediction dGoral ef^Zll2005l ) referring to the magnetic 
field pulse sequence and number of atoms reported for the 
measurements of Fig. [2] The overall picture indicates that 
the quantum kinetic approaches support the interpretation of 
Subsection I V. Al provided that the gas is weakly interacting 
during the evolution period. 

Some measurements were performed also in the regime of 
near resonant evolution fields, B e voive < 157 G. Under these 
conditions, the scattering length as well as the spatial extent 
of the Feshbach molecule are comparable to the mean inter- 
atomic distance of the gas. As the atomic and dimer compo- 
nents are no longer orthogonal during the evolution period, 
the classic concept of a Ramsey interfero meter is expected t o 
break down when B e voive approaches Bo dGoral et ali 120051) . 
Such a phenomenon was observed in terms of a pronounced 
damping of the fringe pattern accompanied by an upward shift 
of its freq uency with respect to the two-body prediction of 
\El mlve \/h dClaussen et aZ.[|2003l) . This shift was subsequently 
attributed to genuinely many-particl e corrections to the intu- 
itive viewpoint il lustrated in Fig.[25l(Duine a nd Stoofl,l2003bt 
G6raleraZ.Ll2005h . 
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FIG. 26 (Colour in online edition) Predicted R amsey fringes be - 
tween three components of a dilute gas of 85 Rb dGoral et all 120051) 
at the end of the magnetic field pulse sequence of Fig. [24] versus the 
evolution time, /evolve- The circles refer to the remnant Bose-Einstein 
condensate, while the squares indicate the final population of cor- 
related atom pairs constituting the burst component. The remnant 
condensate and burst add up to the number of detectable atoms (di- 
amonds). Its difference to the t otal number of N = 17100 atoms 
(dotted line) in the experiments (Donley et al., 2002) determines the 
population of Feshbach molecules. The magnetic field strength in 
the evolution period of B e voive = 159.84 G associated with these pre- 
dictions refers to the conditions of the measured fringes of Fig. [2] 



VI. CONCLUSIONS AND OUTLOOK 

This article has given an overview of a variety of con- 
cepts associated with the description of properties of Fesh- 
bach molecules as well as the techniques for their produc- 
tion in the environment of a cold atomic gas. Its conclu- 
sions can be summarised as follows: Coupled channels the- 
ory provides the foundation for an accurate understanding of 
the diatomic molecular and collision physics. Its predictions 
can be recovered, in an experimentally relevant energy range 
about the dissociation threshold, by two-channel approaches 
characterised in terms of just a few measurable quantities. 
These involve the resonance position and width, the back- 
ground scattering length and the van der Waals dispersion co- 
efficient, as well as the difference in the magnetic moments as- 
sociated with the resonance state and a pair of asymptotically 
separated atoms. Universal properties of Feshbach molecules 
and the low energy collision physics are directly related to 
Wigner's threshold law. Such concepts associated with two- 
body physics treat identical bosons and fermions in different 
Zeeman states as well as distinguishable atoms in essentially 
the same manner. The fundamental transition amplitudes for 
the Feshbach molecular association of an atom pair via linear 
magnetic field sweeps can be represented analytically, in the 
context of two-channel single-resonance approaches. Such an 
exact treatment provides the foundation of the Landau-Zener 
approach to the dimer production in tight traps, such as op- 
tical lattice sites, as well as in cold gases in the fast sweep 
limit. While these limits including their quantum statistical 



phenomena can be understood largely in terms of classical 
probability theory, the saturation of molecule production in- 
volves the many-body dynamics of thermalisation. Given suf- 
ficiently short time scales, beyond mean field approaches can 
provide quantitative descriptions of Feshbach molecule pro- 
duction in Bose-Einstein condensates also in the context of 
non-linear magnetic field variations. 

Closely related to the enhancement of the s-wave scatter- 
ing cross sections discussed in this review are similar mag- 
netically tunable Feshbach resonance phenomena involving 
finite angular momentum quantum numbers, £ > 0. Associ- 
ated scattering amplitudes, fe(p), are usually negligible com- 
pared to the i-wave contribution in cold collision physics due 
to their p 2e scaling with the relative momentum of an atom 
pair. In gases of identical Fermi atoms involving just a single 
Zeeman state, however, only the odd partial waves contribute 
to diatomic scattering, in accordance with the Pauli exclusion 
principle. In such cases, p-wave (I — 1) resonances h ave been 
observed in gases of both 40 K (iRegal et all l2003bl) and 6 Li 



dSchunck et all 120051: IZhang et fl/.ll2004lh providing possibil 
ities of magnetically tuning interactions as well as associating 
dimers. 

Optical Feshbach resona nces dBohn and Julienne! 1 19971 
1999; Fedichev et all 1 19961) provide an alternative approach 
to tune i-wave scattering cross sections in cold gases. This 
technique relies upon coupling the internal states of a pair 
of separated atoms to a closed channel using laser light in- 
stead of a homogeneous magnetic field . According to ex- 
perimental studies in cold gases of 23 Na (Fatemi et al, 2000) 
and 87 Rb Bose-Einstein condensates dTheis et all 12004]) . sig- 
nificant changes of the scattering length tend to cause sub- 
stantial atom loss due to spontaneous photon emission. Re- 
cent predictions indicate, however, that such loss may be ef- 
ficiently suppressed in applications to cold gases of alkaline 
earth metal atoms, such a s bosonic calcium, as well as ytter- 
bium ( Ciurvlo et all 120051) . As these atoms lack hyperfine 
structure due to the absence of nuclear spin, magnetic Fesh- 
bach resonances are ruled out. Their optical counterpart may 
therefore provide the only possibility of tuning interactions in 
such gases. 

Inter-species magnetic Feshbach resonances were observed 
in bos on-fermion mix tures consisting of cold ga ses of 6 Li and 
23 Na dStan et a/.ll2004 as well as 40 K and 87 Rb dlnouve et all 
2004). Their binary physics can be described using the ap- 
proaches outlined in this review, whereas the collisional sta- 
bility of i nter-species Feshba ch molecules requires special 
attention dPetrov et all [2005a). The production of such di- 
atomic molecules is of particul ar interest because of a possibly 
polar character of their bonds (Do yle et q/.l l2b04). Cold polar 
molecules composed of Bose atoms have been produced vi a 
photo-association dKerman et al. , 20041 Wang et all 120041) . 



2002J), the long 



According to predictions (Bara nov et al. 
range, anisotropy and magnitude of the dipole-dipole inter- 
actions give rise to a host of new features in the associated 
fermionic and bosonic superfuids. Arrays of polar molecules 
in optical lattices may provide an opportunity for studies of 
super-solid phases (Gora l et a/.Ll2002l) . It is hoped that heavy 
polar molecules will enable us to improve tests of fundamental 
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physical symme tries, including measurements of the electron 
dipole moment dHudson et al. 1 120021: ISandarsl 1 1 975h . Stabili- 
sation with respect to collisional relaxation of excited molec- 
ular states may be achieved via transfer to their vibrational 
ground states, an ap proach also demo nstrated in the context 
of photo-association (Sag e ef Q/.Ll200l . Such a de-excitation, 
in turn, tends to enhance the polar character of their bonds. 

Few -body scattering phenomena, such as dimer- dimer reso- 
nance s recently observed in cold gases of 133 Cs2 dChin et all 
120051) . may provide possibilities of extending existing tech- 
niques of Feshbach molecule production to more complex 
species. One of the long standing goals in this context is 
an experimental confirmation of predictions associated with 
Efimov's effect in the three-body energy spectrum of identi- 
cal Bose particles (Efimq^, Il970i 1 1 97 lb . The Efimov spec- 
trum consists of an infinite sequence of isotropic three-body 
bound states accumulating at the dissociation threshold, which 
occurs in the limit of infinite s-wave scattering length of 
each two-body subsystem. This scenario is therefore di- 
rectly related to the long range nature o f weakly bound two- 
particle halo states dJensen et a l.. 2004) and may be realised 
via magnetic Feshbach resonances in cold gases. Signa- 
tures of the emergence of such Efimov states are predicted 
to occur in terms of a modulation of three-body recombi- 
nation loss rates as a function of the near resonant mag- 
netic field strength ( Braaten and Hammed l200ll lEsrv et all 
Il999fc iNielsen and Macekl 1 1 9991) . Recently, the first exper- 
imental evidence for this phenomenon has been reported 
( Kra emer et all l2006t) . Such inelastic scattering resonances 
may be exploited to associate meta-stable three-body Efimov 
molecules from cold Bose gases loaded into optical lattices, 
using the magnetic field swee p technique discussed in this re- 
view dStoll and K ohler, 20051) . 
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